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I. INTRODUCTION 



Weak decays of hadrons provide an important source of information on the parameters of the standard model, the 
structure of weak currents, and internal structure of hadrons. Therefore for many years weak decays have been in the 
focus of experimental and theoretical investigations. 

In the last decade, the main emphasis has been laid on weak B decays which allow to measure the unknown 
parametrs of the Cabibbo-Kobayashi-Maskawa (CKM) matrix describing the mixing of heavy c, b and t quarks, and 
the CP-violation: semileptonic and nonleptonic B decays induced by the charged-current quark transition h u,c 
provide a direct access to the CKM matrix elements Vub and Vcb and the weak CP-violating phase. 

Rare semileptonic B decays induced by the flavour-changing neutral current transitions b s and b d measure 
Vts and represent an important test of the standard model and of the physics beyond it. Rare decays are forbidden at 
the tree level in the standard model and occur through loop diagrams. Thus they provide the possibility to probe at 
the relatively low energies of B decays the structure of the electroweak sector at large mass scales from contributions 
of virtual particles in the loop diagrams. Interesting information about the structure of the theory is contained in the 
Wilson coefhcents in the effective Hamiltonian which describes the b ^ s,d transition at low energies. These Wilson 
coefficients take different values in different theories with testable consequences in rare B decays. 

However, along with the interesting parameters of the standard model, the decay rates and differential distributions 
measurable in the decay process contain also quantites related to the presence of hadrons in the decay process. 
Therefore the extraction of the interesting standard model parameters from the decay experiments requires reliable 
information on hadron structure and hadronic amplitudes of the weak quark currents. 

Theoretical description of hadronic amplitudes of the quark currents is one of the key problems of particle physics 
as such amplitudes provide a bridge between QCD formulated in the language of quarks and gluons and observable 
phenomena which deal with hadrons. The difficulty in such calculations lie in the fact that hadron formation occurs 
at large distances where perturbative QCD methods arc not applicable and therefore a nonperturbative consideration 
is necessary. 

The presence of heavy quarks with the masses much bigger than the confinement scale of QCD provide important 
constraints upon the long distance effects due to a new spin-flavour symmetry which emerges in this limit This 
symmetry leads to a heavy quark effective theory (HQET) |^ for QCD with heavy quarks. HQET restricts the 
structure of the expansion of the hadronic transition amplitudes in inverse powers of the heavy quark mass. 

In inclusive B decays the combination of the heavy-quark and the operator product expansions allow one to connect 
the decay rate of the quark bound in a heavy hadron to the decay rate of the free quark. An important consequence 
of the operator product expansion is the appearance of the corrections to the free quark decay rate only at the second 
power of the l/mg |^,||. Therefore, these corrections are numerically small. Providing quite reliable predictions for 
the integrated rates, the theoretical approach based on the OPE turns out to be less efficient for the description of 
the differential distributions. In this case a more precise information on the details of the b quark motion inside the 
B meson is necessary. 

For the exclusive b ^ c decays, HQET provides constraints on the structure of the expansion of the meson transition 
form factors in the inverse powers of the b and c quark masses and leads to the appearance of universal process- 
independent form factors at each order of the l/mg expansion ^,0]. Moreover, heavy quark symmetry provides the 
absolute value of the leading order universal form factor - the Isgur-Wise function at zero recoil (maximum momentum 
transfer) thus allowing for a relible description of the 6 — *■ c decay in this kinematical region. HQET however cannot 
calculate the universal form factors as functions of the momentum transfer. Moreover, the l/rric corrections for the 
b c transiton form factors turn out to be not small because the c quark is not sufficiently heavy. 

For the exclusive 6 — s- w, s transitions, when the final quark is light, relations between the form factors describing 
meson transition induced by different currents emerge in the region near zero recoil Q . In the opposite region of large 
recoils, where the final quark is fast, one can construct another effective theory, so called large energy effective theory 
(LEET) which allows the double expansion of the transition form factors in inverse powers of m?, and the energy E 
of the light quark produced in the weak decay LEET predicts the appearance of several universal form factors 
at the leading order of the 1/E and l/m^ expansion, but does not calculate these form factors, and also does not 
constrain the structure of higher order corrections. 

Thus, heavy-quark symmetry provides important constraints for the form factors but does not allow to calculate 
them in full detail. This task requires a detailed treatment of the nonperturbative effects. 

Theoretical approaches for calculating transition form factors are quark models QCD sum rules [E8H39|, 

and lattice QCD [^-Q. Approaches combining different methods are also extensively used p9|-^. 

In spite of the considerable progress made in the recent years theoretical uncertainties are still uncomfortably large, 
around 10-15%. The main difficulty in obtaining the full picture of the form factors for various decays and for all 

lies in the fact that methods more directly related to QCD, such as lattice QCD and QCD sum rules have only a 



4 



limited range of applicability, whereas the results from quark models are sensitive to the details of the model and the 
parameters used. 

QCD sum rules are suitable for describing the low region of the form factors. Sum rule calculations require 
various inputs such as the condensates or the distribution amplitudes of the light mesons produced in the heavy 
meson decay. The higher region is hard to get and higher order calculations are not likely to give real progress 
because of the appearance of many new parameters. The accuracy of the method cannot be arbitrarily improved 
because of the necessity to isolate the contribution of the states of interest from others. 

Lattice QCD gives good results for the high q^ region. But because of the many numerical extrapolations involved 
this method does not provide for a full picture of the form factors and for the relations between the various decay 
channels. 

Quark models do provide such relations connecting different processes through the meson soft wave functions, and 
give the form factors in the full range. However, quark models are not closely related to the QCD Lagrangian 
(or at least this relationship is not well understood yet) and therefore have input parameters which are not directly 
measurable and may not be of fundamental significance. 

Since quark models are not directly deduced from QCD it is important to match the results for the transition form 
factors obtained within quark models to rigorous QCD predictions for these form factors in the specific limit of heavy 
quark masses. 

The application of various versions of the constituent quark picture to decay processes has a long history. First 
models were based on a relativistic or nonrelativistic [pT|jl^ ] considerations. They could not fully incorporate the 
quark dynamics and spin structure and therefore could not satisfy rigorous relations for the transition form factors 
based of the spin-flavour symmetry in the heavy quark limit of QCD. 

A self-consistent relativistic treatment of the quark spins can be performed within the light-front quark model 
[ ^ , p8| . The model allows for a calculation of the reference- frame dependent partonic contribution to the form 
factor. However, the non-partonic contribution cannot in general be calculated. At spacelike momentum transfers the 
nonpartonic part can be killed by an appropriate choice of the reference frame. Therefore, the partonic contribution 
calculated in this specific reference frame gives the full form factor. At timelike momentum transfers the nonpartonic 
contribution does not vanish for any accessible choice of the reference frame, and the knowledge of the frame-dependent 
partonic contribution is not sufficient to determine the form factor. 

A relativistic dispersion approach to the decay processes based on the constituent quark picture overcomes the 



above difficulties. It has been formulated in |24-27^] and applied to the study of the long-distance effects in many 
pocesses involving the B meson: such as exclusive semileptonic p^- p7| , p9| - |6ll a nd rare [ |62|j63| B decays, the B — > ^Iv 
form factor and the weak annihilation in the rare radiative B ^ pj decay [|64 65 1, nonfactorizable corrections for the 
qo — B^ mixing |6^, calculation of the differential distributions in inclusive B X^lv decays p7| . 

The approach is based on a consistent treatment of the two-particle singularities of the Feynman diagrams describing 
meson interactions. Amplitudes of these processes are given by the relativistic spectral representations over the mass 
variables in terms of the wave functions of the participating mesons and spectral densities of the corresponding 
Feynman diagrams. In particular, meson transition form factors both at spacelike and timelike momentum transfers 
are given by the double spectral representations. 

We discuss in this paper applications of the dispersion approach to various processes involving heavy mesons, laying 
the main emphasis on the calculation of the form factors for heavy meson decays. 

Let us highlight the main features of our dispersion formulation of the constituent quark model: 

1. The physical picture 

The constituent quark picture is based on the following phenomena expected from QCD: 

• the chiral symmetry breaking in the low-energy region which provides for the masses of the constituent quarks; 

• a strong peaking of the nonpcrturbative meson wave functions in terms of the quark momenta with a width of 
the order of the confinement scale; 

• SL qq composition of mesons in terms of constituent quarks. 

As is well known, the qq component of the meson wave function in terms of current quarks dominates the exclusive 
form factors in the deep inelastic region, i.e. for large spacelike momentum transfers. A successesfull description 
of the mass spectrum of mesons as qq states in terms of the constituent quarks obtained in |Q shows that the qq 
approximation works well also in the soft region, where one however has to take into account the transition of the 
current quarks to constituent quarks. The qq meson component in terms of the constituent nuarks leads to a good 
description of the elastic form factors at small and intermediate momentum transfers ||69| , |70|| . Therefore, one can 
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expect the two-particle approximation to be quite reliable for the description of the range of momentum transfers 
relevant for heavy meson decays. 

An important shortcoming of previous quark model predictions was a strong dependence of the results on the special 
form of the quark model and on the parameter values. We demonstrate that once 

(a) a proper relativistic formalism is used for the description of the transition form factors and 

(b) the numerical parameters of the model are chosen properly (we discuss criteria for such a proper choice below), 

the quark model yields results in full agreement with the existing experimental data and in accord with the predictions 
of more fundamental theoretical approaches. In addition, our approach allows to predict many other form factors 
which have not yet been measured. 

2. The formalism 

For the description of the transition form factors in their full g^-range and for various initial and final mesons, a 
fully relativistic treatment is necessary. The dispersion formulation of the quark model provides for such a relativisitc 
treatment and guarantees the correct spectral and analytical properties of the obtained form factors. 

The form factors are given by the double spectral representations over the variables si and S2, the squares of 
the invariant masses of the initial and final Qq pairs, respectively. The integrations in si and S2 run along the 
two-particle cuts in the complex si and S2 planes. The spectral functions of these spectral representations involve 
the wave functions of the participating mesons and the double discontinuities of the corresponding triangle Feynman 
diagrams. 

We start our analysis of the transition form factors from the spacelike region < 0, where the double discontinuities 
can be calculated by means of the Landau-Cutkosky rules. 

The form factors in the decay region q^ > are obtained by the analytical continuation in q^. A specific feature 
of the timelike decay region of g^, is the appearance of the anomalous cuts in the complex si and S2 planes and, 
respectively, of the anomalous contributions to the form factors. Let us point out that the anomalous contribution 
as well as the normal one is completely determined by the wave functions of the participating mesons in the physical 
region. The anomalous contribution is small for small positive q^ but becomes increasingly important as q^ rises. 

The form factors obtained by this procedure obey all rigorous constraints from QCD on the structure of the 
long-distance corrections in the heavy quark limit: they develop the correct heavy-quark expansion at leading and 
next-to-leading orders in l/mg in accordance with QCD in the case when both quarks participating in the weak 
transition are heavy, i.e. have masses much larger than the confinement scale of QCD. 

For the heavy-to-light transition, i. e. when only the initial quark is treated as heavy, the transition form factors 
of the dispersion approach satisfy the relations between the form factors of vector, axial-vector, and tensor currents 
valid at small recoil. In the limit of the heavy-to-light transitions at small the form factors obey the lowest order 
1 / mg and 1/E relations of the Large Energy Effective Theory. 

We want to emphasise once more that form factors for meson decays in the physical decay region q^ > are directly 
calculated through the meson wave functions. 

3. Parameters of the model 

In previous applications of quark models the transition form factors turned out to be sensitive to the numerical 
parameters, such as the quark masses and the slopes of the meson soft wave functions. 

A possible way to control quark masses and the meson soft wave functions is to use the lattice results for the B ^ p 
form factors at large as 'experimental' inputs. The h and u constituent quark masses and slope parameters of the 
B, TT, and p wave functions assuming for them a simple Gaussian form are obtained through this procedure ]59| . The 
Gaussian wave functions of the charm and strange mesons and the effective masses rUc and mg are fixed by fitting 
the measured rates for the decays D [K,K*)lv |6l|. 

With these few inputs, numerous predictions for the form factors for the and decays into light mesons 
are obtained ]6l| ] which nicely agree with the experimental results at places where data are available. The calculated 
transition form factors are also found to be in good agreement with the results from lattice QCD and from sum rules 
in their regions of validity. 

Thus, in spite of the rather different masses and properties of mesons involved in weak transitions, all existing data 
on the form factors can be understood in the quark picture, i.e. all form factors can be described by the few degrees 
of freedom of constituent quarks. Details of the soft wave functions are not crucial; only the spatial extention of these 
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wave functions of order of the confinement scale is important. In other words, only the meson radii are essential for 
the decay processes. 

The paper is organised as follows: 

In Chapter II we present details of the composite system description using spectral representations over mass 
variables. The start with the amplitude of the constituent interaction in the low-energy region and its analytical 
properties. We then the consider the interaction of the bound state with the external electromagnetic field and 
construct the gauge and relativistic invariant amplitude of the bound state interaction. The elastic electromagnetic 
form factor is discussed and the relativistic wave function is introduced. The normalization condition for the wave 
function is obtained which corresponds to the electric charge conservation. 

Properties of pseudoscalar mesons are studied. A single spectral representation for the weak leptonic decay constant 
and double spectral representations for the elastic electromagnetic and the weak transition form factors in terms of 
the wave function are obtained at < 0. For comparison with the light-cone quark model results we rewrite our 
explicitly invariant spectral representations in terms of the light-cone variables and demonstarte the equivalence of 
the form factors obtained within our dispersion approach and the light-cone quark model for < 0. 

The weak transition form factor at q'^ > is obtained by the analytical continuation which discussed in full detail. 

We then consider the case when one of the quarks in the pseudoscalar meson is heavy, mq — » oo and analyse the 
size of the l/niQ effects in the form factors for quark masses in the range of the b and c quarks. 

Chapter III contains a detailed discussion of the weak transitions of the pseudoscalar mesons to pseudoscalar and 
vector mesons. The double spectral representation for both cases are constructed starting from q^ < and going 
to > by the analytical continuation. A procedure of fixing the subtraction prescription in the double spectral 

representations is discussed. 

We then perform the l/mg expansion of the dispersion form factors for the heavy-to- heavy transition to next-to- 
leading order accuracy in l/mg. For the heavy- to-light transition the expansion within the leading order in l/mg is 
developed. Full consistency with the structure of the heavy-quark expansion in QCD for both of these cases to the 
orders considered is verified. 

Spectral representation for the Isgur-Wise function and subleading universal form factors is obtained in terms of 
the wave function of an infinitely heavy meson. Numerical estimates for the universal form factors are given. 

In Chapter IV the choice of the numerical paremeters of the model is discussed and form factors for many D{Ds) 
and B[Bs) to light mesons are calculated. Convenient parametrizations for the calculated form factors are given. 
Strong coupling constants of heavy mesons are estimated by analysing the residues of the form factors in the poles 
located beyond the physical region of meson decay. The strong coupling constant obtaibed by this procedure are 
found in agreement with the results of the direct calculation within our dispersion approach. 

A detailed comparison with the experimental data and results from other approaches is presented. In all cases a 
good agreement with the available experimental data on form factors and strong coupling constants is observed. 

In Chapter V the B ^Iv form factors and the weak annihilation in rare B ^ pj decay is analysed within the 
factorization approximation. A detailed analysis of the contact terms in the weak annihilation amplitude is presented. 
A new contribution missed in the previous analyses of the weak annihilation is reported. Parameter-free numerical 
estimates of the relevant form factors are given. 

Chapter VI contains the analysis of the nonfactorizable effects in the B° — B° mixing due to soft gluon exchanges. 
Assuming the dominance of the local gluon condensate, the correction to factorization can be expressed in terms of 
the specific B meson transition form factors at zero momentum transfers. It is shown that the correction is strictly 
negative independent of the values of the form factors. The form factors are calculated within the dispersion approach 
and numerical estimates for them are obtained. 

In Chapter VII the application of the dispersion approach to inclusive B XJv decays is presented. Spectral 
representation for the integrated semileptonic rate in terms of the B meson wave function is constructed. The 

subtraction prescription is discussed and the absence of the 1 / mg correction in the ratio of the bound to free quark 
decay rates is verified. Differential distributions are calculated in terms of the B meson wave function. 

Conclusion summarises the main results. 
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II. SPECTRAL REPRESENTATION FOR BOUND STATE TRANSITION FORM FACTORS 



This Chapter presents a formaUsm for the relativistic description of hadron form factors within the constituent 
quark picture ||,|9l^. 

Our approach is based on representing the amplitudes of hadron interactions in the form of the dispersion integrals 
over the hadron mass in terms of the hadron soft wave function. This procedure corresponds to a consistent relativistic 
treatment of the leading two-particle singularities of the scattering amplitude and the bound state form factors at 
spacelike momentum transfer q. The form factors at the timelike momentum transfers corresponding to the decay 
process Mini — > M final are obtained by performing the analytical continuation in the variable from its negative 
to positive values. As a result, the weak decay form factors in the kinematical region < {Mini ~ Mfinaif' can be 
expressed through the bound state wave function. 

The application of spectral representations to the description of composite systems has a long history [f76|-p0[. 
Usually, spectral representations in q^ are considered, q being the momentum transfer. In this case anomalous 
singularities in q^ appear in the explicit form as separate contributions to spectral representations. 

We present here an approach to the bound state description based on spectral representations in mass variables. 
Within this approach anomalous contributions turn out to be included in the usual dispersion integrals for the form 
factors at spacelike momentum transfers, relevant for the scattering problems. Form factors at timelike momen- 
tum transfers, corresponding to the decay processes, are obtained by performing the analytic continuation from the 
spacelike region. In this case, anomalous cuts give separate contributions to spectral representations for the form 
factors. 

An important advantage of spectral representations in mass variables is the possibility to introduce in a consistent 
way a relativistic-invariant function which describes the motion of the constituents inside the bound state and which 
can be interpreted as the bound state wave function. We show how this function emerges when keeping in a relativistic 
and gauge-invariant way only two-particle singularities of the Feynman graphs. 

In Section |ll A| we give some details of describing relativistic bound state using spectral representations. 

In Section Tib we present all technical details of the description of a pseudoscalar meson within the dispersion 
approach (leptonic weak decay, two-photon decay, elastic electromagnetic form factor) and demonstrate the equivalence 
of the dispers ion approach and the light-cone constituent quark model | |57t . 

In Section |ll C| we study form factors describing weak transitins between pseudoscalar mesons starting with the 
region of spacelike momentum transfers. As the next step, we perform the analytic continuation to the region > 
and show how the anomalous contribution to the form factor in this region emerges. The origin of this anomalous 
contribution is co nnected with the non-Landau- type singularities of the Feynman triangle diagrams. 

In Section \l Y\ we analyse electroweak properties of pseudoscalar mesons using a simplified parameterization of the 
meson wave function based on the heavy quark symmetry. This allows us to study the dependence of the axial-vector 
decay constant fp and the heavy- meson elastic form factor on the heavy quark mass mg. In particular, we study 
the transition to the heavy quark limit mq oo and discuss the size of the subleading l/mg-corrections for the 
heavy quark mass in the region of b and c quark masses. We calculate weak decay form factors at timelike momentum 
transfers and compare our results with other theoretical analyses and the experimental data. 



A. Bound state description within dispersion relations 

In this section we present some aspects of the dispersion approach to the relativistic description of the bound states. 
For the sake of argument we consider the case of two spinless constituents with the masses mi and m2 interacting via 
exchanges of a meson with the mass /i. We start with the scattering amplitude of the real constituents 

A{s,t)^{k[,k'^\S\ki,k2), S={ki+k2)\ t={ki-k[f, 

kl = kf = ml, kl = k'i = ml. (2.1) 

The amplitude as a function of s has the threshold singularities in the complex s-plane connected with elastic rescat- 
terings of the constituents and production of new mesons at 

s = (mi -I- m2)^, {mi + m2 + ^if , (mi + m2 + 2/^)^ . . . (2.2) 

We assume that an S'-wave bound state with the mass M < mi +m2 exists, then the partial wave amplitude Ao{s) has 
a pole at s = M^. The amplitude A{s, t) has also t-channel singularities with thresholds at t — (rt/i)^; n = 1,2,3... 
connected with meson exchanges. If one needs to construct the amplitude in the low-energy region s > (mi + TO2)^ 
the dispersion N/D representation turns out to be convenient. Consider the S'-wave partial amplitude 
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Aois) 



dz A{s, t(s, z)), 



(2.3) 



where t{z) — —(1 — z)A(s, m^, m2)/2s, z = coa0 in the c.m.s. The Ao{s) as a function of complex s has the right- 
hand singularities related to s-channel singularities of A{s,t). In addition, it has left-hand singularities located at 



(mi -I- 7712) — (nfiy 



1,2,3.... They come from t-channel singularities of A{s, t). The unitarity condition 



in the region s ~ (mi -|- m2) reads 



ImAo(s) =p(s) |Ao(s)|^ p(s) 



A(s,m^,m|) 
IGtts 



(2.4) 



with p{s) the two-particle phase space. The N/D method represents the partial amplitude as Aq{s) = N{s)/ D{s), 
where the function N has only left-hand singularities and D has only right-hand ones. The unitarity condition yields 



D{s) = 1 



dS p{S)N{S) 

TT S — S 



(2.5) 



Assuming the function N to be positive we introduce G(s) — \J N{s) 

. G -xG G/ 




Fig. 1. One of the terms in the expansion of Ao{s) 



Then the partial amplitude takes the form 

Ao{s) = G{.s) [1 + B{.s) + B\s) + B%s) + . . .] G{.s) 

This expression can be interpreted as a series of loop diagrams of Figj^ with the basic loop diagram 

00 

Bis)= I ^m^. (2.7) 
J n s — s 

(mi+ni2)^ 

The bound state with the mass M corresponds to a pole both in the total and partial amplitudes at s = so 
B{IvP) — 1. Near the pole one has for the total amplitude 



A= (fci,fc2lP>T72- — 2(?'l^i'^2) + regular terms 



^xl{k[X)j^Xp{kiM) + --- (2.8) 

where Xp(^ii^2) is the amputated Bethe-Salpeter amplitude of the bound state. The dispersion amplitude near the 
pole reads 

A = N/ D + regular terms related to other partial waves 



(M2 - s)B'{M^) 

^ GUM') 
~ NP-s 

where is a vertex of the bound state transition to the constituents and 



(2.9) 
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J -JT^M^y- (2.10) 

(mi+m2)^ 



The singular terms in Eqs ( p.8[ ) and (2^) correspond to each other and hence 

X,(fci,fe)-G„(.)^^|^. (2.11) 

Underhne that among right-hand singularities only the two-particle cut is taken into account in the constructed 
dispersion amplitude. 

Let us turn to the interaction of the two-constituent system with an external electromagnetic field. The amplitude 
of this process — {k[, fc2| J^|/ci, fc2) in the case of a bound state takes the form 

= k'2) ^^, {p + p'),F{q^)^^Xp{ki, fc2) + . . . (2.12) 

where the bound state form factor is defined according to the relation 

{p'\J^\p)^{p'+p)^F{q^). (2.13) 

The dispersion amplitude with only two-particle singularities in the p^- and p'^-channels taken into account is 
given K% by the series of graphs in Fig.||. 



. G — ^ G -xGG 

P)0(i B(s) ## ## B(s') ^i^p' 

Fig. 2. One of the terms in the series for 



These graphs are obtained from the dispersion scattering amplitude series by inserting a photon line into constituent 
lines. The amplitude has the form 



T^(p',p, q) = 2p^(g)r(s', s, g^) + 



(2.14) 



P ^ s, p ^ s , q = p - p, p^{q) ^ (p q)^ 



The dispersion method allows one to determine the transverse part T{s, s' ,q'^) of the amplitude. Summation of the 
series of dispersion graphs in Fig.0 gives 



T{s\s,q^) = 



^(^) F(.',.,,^)- ^(^') 



1 - Bis) 



1 - B{s')' 



(2.15) 



Here 



T{s',s,q^) = 



dSG(s) dS'G(S') . 2n 
—^—^^A{s',s,q^), 
TT[S — s) Tr(s — S ) 



and A(s',s, 5^) is the double spectral density of the three-point Feynman graph with a pointlike vertex of the con- 
stituent interaction. 

The longitudinal part G is determined by the Ward identity 



C ^ . . {Bis') - B{s)} ■ ^("'^ 



1 - B{s) 



1 - B{s')' 



(2.16) 
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In the region s ~ s' 



M^, Tn develops both s and s' poles, so we write 



G,{M 



/f2\ 



AP - i 



less singular terms 



where 



(mi+m2)^ 



ds Gy jS) d~s'Gy{s') 



A(S',S,g2). 



(2.17) 



(2.18) 



is the bound-state form factor (see ( 2.11 ) and ( 2.12| )). Thus, the quantity {p'\J^\p) corresponds to the three-point 
dispersion graph with the vertices Gv One can derive the following important relation 



A(S',S,0) ^Ti5(s' -s)p{s). 



(2.19) 



This is a consequence of the Ward identi ty w hich r elates the three-point graph at zero momentum transfer to the 
loop graph. Making use of the relations (2.19) and ( 2.10| ), one obtains the normalization of the form factor at zero 
momentum 



F{Q) = 1, 



(2.20) 



which is just the charge conservation condition. The expression ( 2.18 ) gives the form factor in terms of the A^- function 
of the constituent scattering amplitude and double spectral density of the Feynman graph. 

If the constituent is a nonpoint particle, the expression ( 2.18 ) should be multiplied by form factor of an on-shell 
constituent. 

In general, for constructing the spectral representation of the amplitude describing the bound state interaction 
within a two-particle approximation, the following prescription is valid: the spectral density of the amplitude is just 
the spectral density of the corresponding Feynman graph multiplied by Gy . The amplitude obtained through this pro- 
cedure takes into account in a consitent relativistic-invariant way only two-particle singularities of the corresponding 
amplitude. 



B. Quark structure of pseudoscalar mesons 



The pseudoscalar meson P with the mass AI is considered to be a bound state of the constituent quark with the 
mass mi and the antiquark with the mass TO2. Therefore in order to derive the meson interaction amplitude for e.g. 
leptonic decay P Iv, two-photon decay of the neutral pseudoscalar meson 77, and the elastic electromagnetic 

interaction {P{p')\J ^\P{p)) , we start with the corresponding amplitudes of the constituent quark interactions {lv\QQ)^ 
{ll\QQ): a-iid {QQ\J fi\QQ) and single out poles corresponding to the pseudoscalar meson. The amplitude of the QQ 
interaction turns out to be the basic quantity for describing the bound state properties. Near the pole corresponding 
to the bound state with the quantum numbers — 0~, the amplitude is dominated by the 5- wave partial amplitude. 
The latter can be expressed in the two-particle approximation through the dispersion loop graph Bps with the vertex 



)°(fci, mi)z75Q°(-fc2, m2) 



G{s) 



(2.21) 



with a a color index, Nc — ?> the number of quark colors, k\ — , fc| — rn^, and (/ci -f fc2)^ = s 7^ AP. For on-shell 
constituents, the expression ( ^.21 ) is the only independent spinorial structure. 

The dispersion loop graph Fig. ^, which is connected with the meson vertex normalization, reads 




Fig. 3. Meson dispersion loop graph Bps{P^). 
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ds G2(s) 



tt[s — 



(mi+m2)^ 

with pps (s) the spectral density of the Feynman loop graph 

1 



Ppsis), BpAM^) = 1, 



(2.22) 



Pps{s,mi,m2) 



87r2 



dkidk2S{ki — mi)6{k2 — m2)S{p ~ ki — ^2) 



X (-1) Sp \^{ki + mi)ij5{m2 ~ fc2)i75 



Sits 



-{s - (mi - 1112) ) 9{s - {nil + TO2) ), 



(2.23) 



where s = p^ and 



A(s, m^, TO2) = (s + — 777,2) ^ 4sr77i 



Taking into account constituent-quark rescatterings leads to the renornialization of G (Section [I A) and the soft 
constituent-quark structure of the pion is given by the vertex 



)°(fcl,777l)z75Q°(-fc2,nT,2) 



G,(s) 



where G'„(s) = G{s)/B'{AP), such that 



Gl{s)pps{s,mi,m2)ds _ 
7r(s-M2)2 



(2.24) 



(2.25) 



Once the soft vertex (2.24) is fixed, we can proceed with calculating meson interaction amplitudes. 



1. Leptonic decay constant 



Let us consider the decay P Iv. The corresponding amplitude reads 

{Q\A^mp(p)) = ippfp. 



(2.26) 



In this expression = ^"(0)7^ 759^(0) is the axial-vector current where summation over quark colours is implied and 
Jp is the meson axial- vector decay constant. 

We must first take into account that the pion structure is described in terms of the constituent quarks whereas the 
current A^(0) = 97^^759 is given in terms of the current quarks. So let us consider the quantity 



(2.27) 



The bare matrix element which emerges before wc take into account the soft interctions of the constituent quarks has 
the structure 



Al = {<d\A^mQ{k^)Q{k2))bare 

= Q{-k2) [l^.l59W) +P^l5hl{p') + (fci - k2)p^5h°_{p')] Qih) 
If current quarks were identical to constituent ones we would have had 

g0(p2)^l, hl{p')^0, h'L{p')^0. 
It is reasonable to assume that the form factors and /ij. are not far from these values | |8l| ,| 



(2.28) 
(2.29) 



Soft rescatterings of the constituent quarks lead to the series of the dispersion graphs of Fig. 

Bp,(s) 



4_ 




Fig. 4. The series of dispersion graphs for {Q(5|^p(0)|0). 
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The loop diagram Bps is already known, so let us discuss B^{p) The latter contains the pseudoscalar vertex and 
the bare matrix element A^^. The spectral density of the corresponding Feynman graph reads (p^ = s) 



— — ^ / dkidk25{k'l - 'm\)5{kl ~ m\)5{p - fci - /C2) 

xSp (^{t^TsSaCs) + -Pm75^+(s) + (fci - ^2)^75^° (s)} (fci + mi) Z75(m2 - h) 
The trace is equal to 

- 4i(fci^m2 + k2f_tmi)gj^{s) + 4i{ki + fc2)^(fcifc2 + mim2)/i+(s) 

+ 4i(fci - fc2)^(fcifc2 + mim2)h°_{s) 



(2.30) 



(2.31) 



The expression for the loop graph takes the form 



ds G{s) A^/^(s, , mj) s — (mi — m2)^ 



7r(s — p2) IGtts 2s 

(mi+m2)^ 

x[(mi +m2)g!i(s) - s /i^(s) - (m? - m^)/i° (s)]. (2.32) 
The amplitude with the quark rescatterings taken into account has the same spinorial structure as the bare amplitude 



(0|A^(0)|g(fci)Q(fc2)) = g(-fc2) [7m75.9a(p') +P^l5h+{p^) + (fci - fc2),,75/i-(p')] Qih) 



(2.33) 



with 



9a{p') ^ 9"a{p') 
h_(p^)^h^_{p^) 

h+{p')^hl{p')- 



G{p^) f dsG{s) Ai/2(s,mf,m2) s- (mi -7712)2 



l-Bpsip^)] nis-p^) IGtts 2s 
x4[(mi + m2)5^(s) - s /7']^(s) - (m^ - ■ml)h°_{s)] 

The form factor /i+ develops the pole at p^ — A'P since Bps{A'P) = 1. Near p^ = A'P the pole dominates the amplitude 

1 



(0|A,JQQ) = (0|A^|P(p))- 



:{P{p)\QQ) + regular terms. 



M2 - p2 

Comparing the pole terms in ( 2.33| ) and ( 2.34 ) and making use of the relation 

{P{p)\QQ) - ^^G. 



/Nr., 



we find that 



with 



fp = A^jN, 



{Q\A^\P{p))^ip,.fp 
ds Gy{s) A-^/^(s,mf,m|) s - (mi - m2)^ 



(2.34) 



(2.35) 



(2.36) 



7r(s - Af2) 167rs 2s 
x[(mi +m2)g\{s) ~ s h°_^{s) - (m\ ~ Ta\)h^_{s)\. 



Assuming that in reality g\ and are not far from the limit (2.2£), we come to the relation 

ds Gy{s) A^/^(s, ml, ml) (s — (mi — m2)^) 



ff 



^(mi+m2)g!l(Af2) 



(mi+mg)^ 



7r(s - AP) 



Stts 



(2.37) 



(2.38) 
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2. Two-photon decay of a neutral pseudoscalar meson 



Let us consider the two-photon decay of the neutral pseudoscalar meson Pq. Its constituent quark structure is 
described by the vertex 

^'^'^G.. (2.39) 



The rate of the decay Pq (p) 2j can be written as 

T^^a^M^gl^^, gp^^ = Gp^^iM^ ,0,0), (2.40) 
where the form factor Gp^-y is connected with the amplitude 

(oi j^r (92) j^r (93)i^b)> = 26„,,30,^3Qf^g3^^Gp,,(A^^ qi qi). (2.41) 

The electromagnetic current J^™(0) — q{0)jfj_q{0) is defined through current quarks, whereas the meson structure is 
described in terms of the constituent quarks. So, for calculating the meson amplitude the constituent quark amplitude 
of the electromagnetic current is necessary. The latter is assumed to have the following structure 

(Q(fc')k"(0)7^g(0)|Q(fc)) = g(fc')7pQ(fc)/c(9'), q = k'-k. (2.42) 

The constituent charge form factor fc{q^) is normalized such that /c(0) — e^, t he co nstituent charge. The anoma- 



lous magnetic moment of the constituent quark is neglected in the expression (2.42), but it can be included into 
consideration straightforwardly. 

Hereafter we skip the intermediate steps and consider only the residue of the constituent interaction amplitude 
which determines the bound-state amplitude. From now on we shall use the following notations: 
f) is the off-shell momentum = s which is used for the calculation of the imaginary parts of the Feynman diagrams, 
p is the on-shell bound state momentum, p^ = IvP . 

The single dispersion representation for the form factor Gp^~f reads 

Gp,,{M\ql,ql) = fM)fM) J ^j^^Ap,,{s,qlql). (2.43) 
Here Ap^^ is the spectral of the triangle Feynman graph of FigjH and can be obtained from the following relation 



Fig. 5. The graph describing the decay P^' —> 77. 
' dkidk2dk36{p - fca - ^3)^(^2 - fci - 93)- 



xS'p (^i75(m - A;3)7q2(to -I- fci)7a3(m -f ^2)) , (2.44) 
where p^ — s. The trace reads 

Sp (ijsim - fc3)7a2(™ + fcl)7a3(™ + ^2)) = 'imea2a3l320392^l3^^ (2-45) 
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and we find 



(2.46) 



Substituting (2.46D into (2.43), one obtains Gp-y^ for off'-shell photons. For real photons one finds 



An 



'el 



irn? 



dsGy{s) 1 / 1 + ^l-im^/s 



(2.47) 



3. Elastic electromagnetic form factor 



The elastic electromagnetic form factor of a pseudoscalar meson with the mass M is given by the following matrix 
element 

{P{p')\JrmPip)) = ip'+P), F^\q'), (2.48) 

where ~ p'^ — AP, q — p~ p' and < 0. The form factor F'^^{q^) describes the amplitude of the photon emission 
from the bound state which depends on the three independent Lorentz invariants. It is convenient to choose such 
invariants as the squares of the three external momenta p^ , p'^, and q^ . Therefore the form factor F'^^{q^) is in fact 
the function of the three variables F'^^ [p^ , p'"^ , q^) , but considered for fixed values of the two invariants p^ — M^, 

Assuming again the structure of the constitent quark matrix element of the the electromagnetic as given by Eq. 
(2.42), the meson elastic charge form factor can be written in the form 

F'\q^) = h{q')H{q\ mlml) + f2{q^)H{q\ml ml) (2.49) 

in terms of the form factors H. The quantity H{q^ ,m\^ m^) describes the subprocess when the constituent 1 interacts 
with the photon, while the constituent 2 remains spectator (FigJ|). 



q 




Fig. 6. Triangle diagram describing //(g^, mf , m|). 

Clearly, the form factor H{q^,m\,m\) also implicitly depends on the invariant variables p^ and p'^, such that 
H{q^,ml,ml) = H{q^,p'^ = M'^^p''^ = M"^ ,m\,m\). As well known Q the form factor H{q'^,p'^,p''^,ml,ml) is an 
analytic function of the external mass variables p"^ and p'^ and at < satisifes the double spectral representation 

Hiq^mlmj) = J ""^1;^) ^^(-^ ^' q'\rm,m„m,). (2.50) 

Here Ay is the double spectral density of H over the variables and P'^. According to the Landau-Cutkosky rules. 
Ay can be calculated through the following procedure: one should place all internal particles on their mass shell, 
ki = ml, k'^ = , k\ = rr3^ but go off the mass shell for the variables p^ and p'^, i.e. instead of the on-shell external 
momenta p and p' consider off-shell external momenta p and p', such that p^ = s, p'"^ = s' , and {p' — p)^ = q^, but 
p' — p = q ^ q- The double spectral density A is then determined from the following relation 
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— / dkidk[dk26{kl - ml)S{k'^ ~ ml)S{k^ - ml)S{p - fci - k2)S{p' ~ k[ ~ fca) 



X{-1) Sp [{k[ +mi)7p(fcl +7711)^75(777,2 - fc2)«75) = 2pp((j)Ay(s',S,9^|777l, 7771, 7772) 



with 



The trace reads 



-Sp [{k[ + mi)-ff,{ki + 7771)75(7712 - fc2)75 



2fciu(s - (™i - "^2)^) + 2fci^(s' - (7711 - 7712)^) + 2fc2^g^ 



Multiplying both sides of (2.51) by and using (2.52) one obtains for q <Q 



Ay(s',s,g^|777i, 7771,7712) 



(s's + (s' + S — 9^)7712(7771 — 7772) — (77I1 + 77l2)(77ll — m2)^) 



4A3/2(.s',s,g2) 

X6I(S - (77I1 + 7772)^)e'(s' - (77I1 + 7772)^) 

x^^ {-q^{s' + s-q^ + 2{ml - 7772))^ + A(s', s, q^){q^ - Ami)) 



with A(s', s, g2) = (s' + s - q^Y - 4s's- 
At = one finds 



ls.v{s' ,s,q^ = 0| 777i, mi, 7M2) = 7rpps(s, 7711, 7772) (5(s' - s). 



and 



F^'(O) = (61+62) 



7r(s - APy 



■Pps{s,mi,m2) 



ei + 62. 



(2.51) 



(2.52) 



(2.53) 



(2.54) 



(2.55) 



This relation is the direct consequence of the Ward identity and corresponds to the electric charge conservation. 



4-. Dispersion approach in terms of the light-cone variable 



For some applications and for comparison with the light-cone technique, it is convenient to rewrite our explicitly 
relativistic-invariant spectral representations in terms of the light-cone variables 



fc+ = ;^(fco + fcz); 



(2.56) 



Most easily this can be done by in troducing the ligh t-con e variables directly into the integral representation for the 
form factor spectral density ( 2.51 ). The variables ( 2.56 ) should be connected with some specific reference frame, 
which can be specified by fixing components of the physical momenta p and q. It is convenient to choose the reference 
frame in which 



P± =0, q+ ^ 0, qj_ = -q . 



(2.57) 



Notice that this choice is only possible in the re gion q^ < 0. 

The choice of the reference frame in the form (2.57) allows us to choose components of the momenta in a convenient 
way. For the physical on-shell momenta|^ 



^We use the notation = (a+,a_,aj 
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M2 

P= (P+^o '0) 

, Mj + qj 

p = — -,-q±) 

2p+ 



(2.58) 



For the dispersion ofF-shell momenta: 



P = 71 — .0) 



p' = b-t 



q = (0, 



2p^ 



.-qi-) 



/ 9 

s -qi 



2p^ 



q±). 



(2.59) 



A specific and good feature of this choice is that the (+) and (_L) components of the on-sheU vectors and the 
corresponding off-sheU vectors are equal to each other, such that all off-shell effects in the description of a bound state 
are completely shifted t o the (— ) component of the momenta. 

Now, let us consider ( 2.51 ) and set /i = + in both sides of this equation. Using the result for the trace ( 2.52 ) for 
/i = +, and performing the k- integration, we obtain 



Ay(s', s, g^|TOi, mi, m2) 



1 



dx(P k I 



s' 



1-x 



IGtt J x{l — x) 

m\ _ ikj_+xq±f 
X x{l — x) 



5 s 



1-x 



x{l — x) 



(s' + s — 2(mi — TO2)^ 



1-x 



(2.60) 



Here x = k2+ /p+ and kj_ — k2±_ . 

Substituting ( |2.60 ) into (2.50) and performing s and s' integrations, one derives 

Hiql,m,,m2) = I dxd^k^i^ix, kM^, k, + xqM^, fc„ ,,) 

where the radial light-cone wave function of a pseudoscalar meson is introduced 



(2.61) 



tp{x,k±) 



Gv{s)^/s - (mi - TO2)^ 
7r3/2V8(s-M2)^a;(l-a;7' 
mJ m2 kj^ 



1-x 



x{l — x) 



and 



f3- 



s - (mi - m2)^ + kj_q^/ (1 - x) 



y/s- (mi - ra2Y\Js' - {mi - m2)^ 



Piqi_ - 0) - 1 



(2.62) 



(2.63) 



The function /3 accounts for the contribution of spins. It is different fro m un ity at g_L ^ because both the spin-nonflip 
and spin- flip amplitudes of the interacting quark contribute. The Eq.( ^.55| ) is the normalization condition for the soft 
radial wave function 



/ 



dx d'^k± \ip{x, A:^)p 



1. 



In terms of this wave function, the pseudoscalar meson axial-vector decay constant fp is represented as 



fp = 9 A 



^/27r3/2 



dx d'^k± tp{x, k±) 



m2(l — a;) -|- mix 
\/s- (mi - m2)2 



(2.64) 



(2.65) 



This expression can be easily deduced by introducing the light-cone variables into the dispersion representation ( ^.30 ), 
making use of (|3^) and examining the fJ- = + component of the axial current. 
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Similarly, introducing the light-cone variables into (2.44) yields the following expression for gp^^ 



The same expressions for the speudoscalar meson elastic form factor, leptonic decay constant, and the two-photon 
decay constant as (| 2.6i p-( ^.66 ) were derived within the light-cone quark model in refs | p^ , p3[ , with our Gy{s) just 



equal to ho{P) of rcf. ||14|. As we see later, the light-cone quark model and the dispersion approach also lead to the 
same expression for the form factor describing weak transitions between pseudoscalar mesons at < 0. This is 
not at all strange as both of these approaches are based on the assumption of dominance of the qq components in the 
description of the meson properties. 

However, the dispersion approach has several advantages compared to the light-cone quark model. In particular, 
the light-cone quark model faces at least two difficulties in applications to inelastic transitions: The first problem 
emerges already in the spacelike region of the momentum transfers in considering transitions induced by the current 
of higher spin. It is related to the proper choice of the current components to be used for the extraction of the form 
factors (so-called 'good' and 'bad' components and the problem to satisfy the angular condition, see for discussion 
and references). 

The second problem emerges for the description of hadron transition processes at > 0. The light-cone quark 
model can consistently take into account the spectator contribution, whereas contributions of the so-called Z-graphs 
(also called non-partonic, or pair-creation subprocesses) are in general beyond the scope of the light-cone treatment. 
The Z-graph can be calculated only for some exceptional simple forms of the light-cone wave functions. Unfortunately, 
in the light-cone treatment such contribution is inevitably present for processes in the region of timelike momentum 
transfers where it cannot be suppressed by an appropriate choice of the reference frame. 

As we shall see, both of these problems find their solution in the dispersion approach. 

C. Form factors of meson transitions 

In this section we examine the electroweak transitions of pseudoscalar mesons. First, we derive the transition form 
factors at < in the form of the relativistic dispersion representations. Second, these dispersion representations 
allow us to perform the analytic continuation in q^ and derive the form factors of semileptonic decays of pseudoscalar 
mesons at > 0. 

1. The pseudoscalar meson transition form factor at < 

The amplitude of the weak transition of pseudoscalar mesons AIi A'h (Fig.0) is determined by the two form 
factors and F_ 

{p2.M2\V^\puMi) = {pi + p2)f,F+{q^) + {pi ~ p2)f,F^{q'') 
{p2,M2\A^\pi,Mi) =0, 

pI = MI pI^MI p^-p2^q. (2.67) 

The weak currents are defined in terms of the current quarks 

V, = (0)7^92(0), 

- 9? (0)7^7592" (0). (2.68) 

The structure of the mesons is described in terms of the constituent quarks by the vertices 

Q2ik2)ij5Q3{-h) ^ f^2^ 

Ml . Gyl(p^), 

M2 : --j= Gv2(P2) (2.69) 



and Giii 2 are normalized according to Eq. ( |2.25| ) 
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For calculating the transition amplitude (2.67) we again need the constituent quark matrix element of the weak 
current which is taken in the form 



(Qi(fci)|gi(0)7M'?2(0)|Q2(fc2)> = QiikihMk2)f2i{q^) 



(2.70) 



As in the case of the elastic form factor, the transition form factors F± can be treated as the functions of the three 
invariants Pi, and q^, taken at p\ = and p\ = Af|. The form factors are analytic functions of the external 
mass variables p\ and p\ and can be represented by the following double spectral representations 



J 7r(si - Mf ) TT 



dsiG„i(si) d s2Gy2{s2) 

(,S2 - M|) 



A±(si,S2,g^|mi,m2,TO3) 



(2.71) 



Here A± are the double spectral densities of the triangle Feynman graph corresponding to Fig.|^ in the p\ and 
P2-channels 




Fig. 7. The dispersion graph for the decay {-P2|V^''|A) 



To calculate Aj- at < 0, we put internal particles on their mass shell, but consider external off-shell momenta 
pi ^ P2 + q, Pi = si, P2 = S2, = q^ ■ Then the double spectral densities can be obtained from the following equation 



1 

8^ 



dkidk2dk^5{k\ — m\)5{k1 — m^)5{k2 — in'^)5{pi 



h)S{p2 



x(--l) Sp (ki + mi)7p(fc2 + m2)i75(TO3 - fc3)«75 



= (pi -t-p2)^jA+(si,S2,9^|mi,m2,m3) + {pi - p2)^A_(si, S2, TO2, ma) 

The trace has the following form 

-Sp (^{ki + mi)jf_,{k2 + m2)z75(TO3 - ^3)«75 



(2.72) 



2fci^(si - (m2 - ms)^) + 2A:2^(s2 - (wa - toi)^) + 2/c3^(q^ - (mi - 7712)^). 

S±(si,S2,g2) 



Explicit calculations give 

A± (si , S2 , I nil , m2 , TO3) 



A(si,S2,g2) 



-A{si,S2,q |TOi,m2,m3) 



(2.73) 



(2.74) 



where B± are the following polynomials 



S+(si,S2,9^) = &+(si,S2,g^) {a{si, 7122,1713) + a{s2,m3,mi) - a((j^, mi, TO2)} 

+a(q^, mi, m2)A(si, S2, 9^), 
i?_(si,S2,g^) = 5_(si,S2,g^) {a(si,m2,m3) + a(s2, m-3, "?,i) - a(g^, mi, m2)} , 

+ (a(s2,m3,mi) - a(si,m2,m3)) A(si,.S2,g^) 
6+(si,S2,9^) = -g^(si + S2 - + mi + m2 - 2m3) - (mi - m2)(si - S2), 
6_(si, S2, 9^) = (mi - m2)(2si + 2s2 - g^) - (si - S2)(si + S2 - + mi + mj - 2m3), 
a(s, /ii, ^2) = s - (mi - A*2)^- 



(2.75) 



19 



In the above relations we have introduced A, the double spectral density in and p2-channels of the Feynman triangle 
graph r{pi,p2, q^) with the scalar constituents defined according to the relation 

1 f dkidk2dk35{Pi - - fc3)(5(P2 - fcg - ki) 



(27r)4 i J {mj-kf- iO){ml - fcf - iO){ml - fc| - iO) 

/dSl (lS2 A / 2l \ /f, rT/,N 

— 7 2\^('Si''^2,g mi,m2,TO3). (2.76) 

At < 0, A is given by the following integral 

A{si,S2,q^\mi, 1712,1713) = — J dkidk2dk36{pi - k2 - k3)S{p2 - fcs - ki) 

x6{kl - ml)6{kl - ml)6{kl - mj). (2.77) 
Introducing the light-cone variables we derive a useful expression 



Af 2, \ ^ f dXldX2dX3 ^2u Kt \ 

A(si,S2,g mi, 7712, ma) = 77— / d k3^d{xi - X2)d{l - xi - X3) 

IGtt J X1X2X3 



x^L-!^-!!i-^Vf«2-^-^-^^^^^tW\ (2.78) 

V X2 X3 X2X3J \ Xi X3 X2X3 j 

Hereafter we denote Xi = A;j+/p+,pi-)- = p2+ P+, — = <7^; a; = 2:3, k^ = fcaj,. For < one obtains 

2, ^ 6»(62_(si,s2,g^) - A(si,s2,g2)A(g2,m?,mi)) 
A(si,S2,g^ mi,m2,m3) = ^. (2.79) 

The solution of this 0-function gives the following allowed intervals for the integration variables s\ and S2 

(mi + m3)^ < S2, 

sr(s2, g^) < si < (s2, g^); (2.80) 



where 



±1 2A S2(m? + mi-g2) + g2(m?+mi)-(m?-mi)(m?-mi) 
Si l'S2, g j — 



2m\ 

AV^(.2,mi,mf)AV^(g^mf,mi) 
^ "h^l ^^-^^^ 

The final dispersion representation for the form factors at g^ < takes the form 

00 s^(s2,g^) 

P / 2N . . 2n /■ C?g2g^2(52) A C?glg^l(si) -B±(si, 52, g'^) 

^^(^^=^-^'^) y 4^2 - Ml) y .(si-Mf)A3/2(si,S2,g^) ^'-''^ 

(mi+m3)2 s^(s2,9=^) 

This representation will be the starting point for the consideration of the meson decays in the next section. 
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The light-cone representation 



For comparison with the hght-cone quark model results it is convenient to represent the spectral representation 
( 2.82 ) as the integral over the light-cone variables. Following the s ame l ines for the elastic form factor, we turn back to 
the equation (^.72| ) and again make use of the light-cone variables ( 2.56| ), choo sing t he reference frame 5+ — 0,pi± = 
(see |69| , [73| , |75[ for details). Setting /i — + and making use of the first line of ( 2.73| ) gives for 



A+{si,S2,q \mi,m2,m.3) = -— / 

iDTT J Xl 



dxidx2dx3 2, xf \x^^ \ 

d k3^d[xi - X2)o{l - Xl - X3) 



X2X3 



xS si 



X2 



X3 



X2X3 



6 S2 



ml ik3±+X3q±f 



Xl 



3 

X3 
2 



X2X3 

„2 



{xi{si ~ {m2 ~ TOa)^) + 2:2(52 - (™i - ^3)^) + - (™i - "^2)^)} (2.83) 



Substituting ( 2.83 ) into ( 2.71 ) yields the following expression for the form factor F+ 

dxd'^k± Gi,i(si) Gv2{s2) 



F+{qi) ^ f2l{qi) 



167r3a;(l - x) n{si - Mf) n{s2 - M|) 



si + S2 — {mi — ms)^ — (m2 — ms)^ + 



1 - X 



{-ql - {m-i - ms)^) 



(2.84) 



Introducing the radial light-cone wave funcion according to ( 2.62 ) leads to the familiar light-cone expression (cf. [g^ ) 

(2.85) 



F+iq'i) = f2iiq'i) J dxd'^k±tpi{x,k±^)ip2{x,kj_ + xq±_)f3+{x,k^,qj_), 

with /3+ given by the expression 

_ si + S2 - {mi - TO3)2 - (m2 - ms)^ + ^(-q2 _ (^j^ _ 7712)^) 
P+ ^ . — , 

2^ySl - (to2 - TO3)^VS2 - (TO3 - TOl)^ 

_ (miX -I- 7713(1 — x)){m2X + 7713(1 - x)) + k^{kj_ + xqj_) 
x{l - x)^/si - (777,2 - f"3)^Vs2 - [ms - 7ni)2 

The same expression for J^+(g^) was obtained in within the light-cone quark model. The (_L) component in Eq. 
(2.72) allows to determine the form factor F^{q'^). The (— ) component of Eq. (2.72) cannot be used for isolating the 
form factor beca use o f the difference between the (— ) components of the on-shell and the off-shell vectors as given by 
Eqs. (I^H) and ^M)- 



2. Transition form factors at cf > Q 

For the description of decay processes the form factors in the regio n < g'^ < (Mi — 1/2)^ are necessary. For 
deriving the form factors at > the dispersion representation (2.82) turns out to be a convenient starting point. 
In general, this representation has the following form 



rfsiG'^i(si) ds2G^,2{s2) B{si,S2,q'^ ^ 
7t{si - Mf) 7r(s2 - M|) A(si, S2, g2) 



El/ 2\ f I 2\ I ""'i'-f-uiV-'i; ",''2^u2V''27 -t^V'^ii ''2, y 2| N |c^ or\ 

F{q ) = /2i(g ) / — TTT:— 2V^(*i'^2,9 |mi, 7773, 7773) (2.86) 



where A is the double spectral density of the Feynman graph F with scalar constituents Eq. ( ^.76 ). This double 



dispersion representation defines the analytic function of q^ both at negative and positive values. It is important to 
point out that the functions Gy{s) have no singularities in the right-hand side of the complex s-plane |72[| , and B 
and A are polynomials. So the details of the dispersion integration at > are determined by the behavior of the 
quantity A. 

A detailed discussion of the double spectral representation and the anomalous singularities can be found in [pO| in 
connection with the deuteron elastic form factor. Anomalous singularities in decay processes were considered inT2q] 
for the case mi = 7173 — 0. We perform the analysis for arbitrary nonzero masses. 

Let us start with the single dispersion represenation for F in the variable p| . A standard calculation yields 



21 



ds2 



7r(.S2 - pI) 



o-2(p?,s2,g^), 



(2.87) 



where 



0-2(^1,32,9 ) = a+{p^,S2,q ) - (Pi,S2,g ), 
1 

CT^(si,S2,g^) 



167rAi/2(si,s2,g2) 



X log (-S2(si + - S2 +ml + ml- 2ml) " (^i " " "^s) ± X^^'^{s2,ml, ml)X^/'^{si,S2,q'^)j ■ 

Hereafter we assume TO2 > mi. The single dispersion representation reproduces the exact value of the Feynman 
integral ( |2.76| ). 

Let us consider the function (J2{Pi, S2, 9^) as the analytic function of the variable si = p\ at fixed S2 and > 0. 
As S2 < S2 such that 



Q (/^ + m\ — 

2 = 



2vr 



g-^ + — m^ 



+ (m§-m^), < (m2 - mi)^ 



(2.88) 



both of the functions cr^ and have square-root branch points on the physical sheet at — ^ \f(p)'^ 

sf = (-\/s2 + 'v/?)^' connected by the cut, see FigJ|. 







s7 




s\ 






St 


s\ 




s7 





Fig. 8. The location of the singularities of (T2 in the complex si plane as a function of S2 for > 0: a). S2 < S2; b). S2 > s^. 
Solid lines cuts located on the physical sheet, dashed lines - on the second sheet. 



The function has in addition a logarithmic cut from to on the physical sheet. Here given by ( ^.8C ) 
are the zeros of the argument of the logarithm in a2 ■ The function cr^ also has a logarithmic cut from to but 
on the second unphysical sheet of the Riemann surface of the square- root (dashed line in Fig.||a). 

The square-root cuts cancel in a2 = cr^ — a2 , and the logarithmic cut is the only singularity of (T2 on the physical 
sheet. Notice however that 172 has the two cuts from to and from sj~ to of the second sheet. 

The situation changes at S2 — S2 which is determined by the condition s{*(s2) = Si{^2)- For this value of §2 the 
logarithmic and square-root branch points coincide. As S2 increases, S2 > S2, ^^e branch point of cr^ moves up 
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through the square-root cut onto the physical sheet. At the same time the branch point sj~ of cr^ moves onto the 
second sheet (Fig.||b). Hence, the function cr^ acquires the logarithmic cut from to on the physical sheet, and 
(T^ still has the logarithmic cut which now lasts from sf" to . Both of the functions have also square-root cuts from 
Si to s^. In the difference (T2 = cr^ — cr^ the square-root cuts cancel each other, but the logarithmic cuts add. The 
resulting expression for the double spectral density takes the form 



A(si,S2,g^|mi,m2,TO3) 



0{s2 - (mi + m:iY)e{s^ < si < s'l) 

16Al/2(si,52,g2) 

29{q^)e{s2-sl)e{s^<s^<s-) 
16Ai/2(si,s2,g2) 



(2.89) 



One can check the double dispersion representation ( 2.76 ) with the spectral density A given by ( 2.89| ) to reproduce 
correctly the Feynman expression. The first term in ( 2.8E ) relates to the Landau- type contribution emerging when 
all intermediate pa rticle s go on mass shell, while the second term describes the anomalous contribution. 

The expression dH) for A is derived for TO2 > mi implying the 'external' S2 integration, and the 'internal' 
si-integration. The location of the integration region for this case is shown in Fig. ^ 



(m 1-1-013) 



S2 



S°(q2) 



(mi-i-mj) 



S°(q2) 



S2 



(m2+m3)" 



Si 



SI 



(m2H-m3) 



sf 



sV 



Si 

3- 



st 



Fig. 9. The integration region for m2 > mi and > for the order of integration J ds2 J dsi: a. S2 < S2 b. 
lines - cuts on the physical sheet, dashed lines - cuts on the second sheet. 



S2 > S§. 



Solid 



In addition to the quantity A, the spectral density of the representation ( 2.86 ) involves the factor 1/A(si, S2, q^) which 
is singular at sf' , the lower integration limit in the anomalous term: because 

As it has been discussed in p8|] , in this case an accurate application of the Cauchy theorem yields the subtracion term 
in the non-Landau contribution. Representing 02 as a contour integral, we must take into account the nonvanishing 
contribution of the small circle around the point s^. Underline once more that the presence of the factor Gi,i(si) 
does not change the argumentation as the function Gr,{s) has no singularities at si > (m2 -I- ma)". The final properly 
regularized representation for the form factors at < (7^ < (m2 — mi)^ takes the form (omitting the constituent 
transition form factor f2i{(f')) 



ds2Gy2{s2) f dsiGyi{si) B{si,S2,q'^) 



(mi+ma)^ 



7r(s2 - Mi) J 7r(si - M[) 16A(si, S2, g^) 



2eiq') 



dS2Gv2{s2) f dsi 

7r(s2 - Mi) J 167r(si - sf )3/2 



ai(si)B(si,S2,g') ai(sf)B(.sf,S2,g2 



(si-sf)3/2( 



Si 



Ml) (sf -sf)3/2(sf -Mf) 



(2.90) 



It should be pointed out, that although the representations (2.86) and (2.90) were derived for the case of pseudoscalar 
mesons, transition form factors of any hadrons have the similar structure. A particular choise of the initial and final 
hadrons yields a specific form of the function B. In the next chapter we shall calculate the double spectral densities 
for the form factors describing the pseudoscalar to vector meson transitions. 
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Bound state vertex and hound state wave function 



It is convenient to introduce the bound state wave function ^p{s) related to the bound state vertex as follows 

ip{s) = G^{s)/{s-M^). (2.91) 
In the nonrelativistic limit one finds </?(s) ~ 5'Arjj(fc), where the variable k is connected with s 




In the nonrelativistic limit the elastic electromagnetic form factor of a pseudoscalar meson takes a well-known form 

F{f ) = J ^NR{k)^NR{k + q}dk, cf = -q\ (2.92) 
The r elationship between the bound state wave function ip{s) and the light-cone wave function is given by the relation 

(Hi)- 

It should be pointed out that the analytic properties of if{s) in the region of s near the two-particle threshold 
are different for a truly bound state with a negative binding energy (like a deuteron as a two-nucleon system) and a 
bound state in a confined potential with a positive binding energy (like a meson as a qq bound state) : in the deuteron 
case Gv{s) is a smooth regular function near the threshold and the pole at s = is located only slightly below the 
two-particle threshold at s = 4M^. Therefore (p{s) is strongly peaked near the threshold and it is more convenient 
to analyse the deuteron form factors in terms of the vertex functions Gy{s). 

For the confined potential the situation is different: the pole at s = would have appeared in the physical region 
as the pole in (p{s) a.t s — A/^ > (mQ + mq)'^. However this does not happen. As well known from e.g. the behaviour of 
the bound-state wave function in the harmonic oscillator potential, the function if{s) is a smooth exponential function 
of s above s = {mq + niq)'^. This means that the would-be pole in ip{s) at s = is completely washed out by the 
interaction. Therefore, for the analysis of the meson transitions <p(s) turns out to be more appropriate than Gy{s). 



D. A simple model for pseudoscalar mesons 



We are now in a position to apply the developed formalism to the analysis of the properties of pseudoscalar mesons 
and to the direct calculation of the decay form factors. To this end we must specify the parameters of the model, i.e. 
constituent quark masses and wave functions of pseudoscalar mesons. 

We consider in this section a simple choice of quark model parameters for pseudoscalar mesons vr, K, D, and B 
which gives a relalistic description of this sector and allows us to illustrate main features of the weak form factors and 
study the transition to the heavy-quark limit. We leave a detailed discussion of criteria for fixing parameters of the 
model, including also vector mesons, and analysis of the corresponding transition form factors for later sections. 

For a pseudoscalar meson built up of quarks with the masses mg and iriq the wave function ip{s) = G^ (s) I {s — M^)) 
can be written in the form 



^ Js^-{ml-mlf ^ ^i/2( 2 2) 

^(^)-^ ^, , \, 4u<kl k^ 4-7r^- (2.93) 

\/2 - {mq - mqY s-^l^ 2y/s 



The normalization condition ( 2.25| ) for Gu is equivalent to the following normalization condition for w 

w^{k)k^dk = 1. (2.94) 



The function w is the ground-state S'-wave radial wave function of a pseudoscalar meson for which we choose in this 
section a simple exponential form 

w{k) = exp (-4afcV/4) (2-95) 

where /ip = rnqruq/ [mq -\- ruq) is the reduced mass. The parameterization of the wave function in the form ( ^.95| ) 
is inspired by the nonrelativistic quantum mechanics and is convenient for the analysis of the dependence of the 
observables on mq, in particular for analysing the case mq oo. 
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In the nonrelativistic quantum mechanics a bound-state wave function is determined by the motion of the particle 
with the mass /xp in the potential independent of masses, and thus a does not depend on the masses as well. 
Relativistic effects destroy this simple feature of the wave function. In QCD the situation is much more complicated 
because additional dimensional quantities such as h.QCD and the condensates appear. So, a should be considered 
as some unknown function of the quark masses. It is possible to obtain the information on the behavior of a as a 
function of mq at fixed vriq = niu.d — 0.25 GeV in the two regions: at small mg and mg — > oo. 



Table 1. The constituent quark masses and the calculated fp for a = 0.02. 



quark 


quark mass,Ge\^ 


meson 


meson mass, GeV 


fp,MeV 


u,d 


0.25 


7r+ (ltd) 


0.14 


130 


s 


0.40 


K+{us) 


0.49 


160 


c 


1.80 


D+{cd) 


1.87 


234 


b 


5.20 


B+{ub) 


5.27 


202 



At mg < 0.5 GeV the value of a can be determined by describing the data in the light-meson sector. The light- 
quark masses given in Table |l] and a-^ = uk = 0.02 provide a good description of the data on fx-, and the elastic 
form factors. Fig. |l^. The meson decay constants and form factors are calculated with the values g\{M'^) = 1 and 
/c(g^) = /c(0), respectively. 




Fig. 10. The tt^ and elastic form factors evaluated for = 0.02. 
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In the region mg — > oo the behavior of a can be found on the basis of the heavy quark symmetry. To this end, let 
us consider the amphtudes of the elastic and inelstic transitions between pseudoscalar mesons consisting of heavy Q 
and light q quarks. 

For the case of the transition between two heavy mesons it is convenient to introduce instead of the dimensionless 
recoil variable uj = vv' , where v and w' are the 4- velocities of the initial and final mesons, respectively, and to analyse 
the transition in terms of the velocity-dependent form factors. 

For the elastic-transition amplitude 

{M,P'\Qj^Q\M,P) = {P' + P)^F,i{q^); < (2.96) 

the recoil takes the form 

- = 1-^>1, (2.97) 

and the velocity-dependent form factor just coincides with the elastic form factor 

Kiiiu) = F.iiq^). (2.98) 

The function hei can be expanded in powers of the variable lo — 1 near zero recoil point a; = 1 as follows 

heiiu;) = 1 - pliiij - 1) + 0{{lu - 1)2). (2.99) 

For the amplitude of the inelastic transition between two pseudoscalar mesons 

{M2,P2\Q27t.Qi\Mi,pi) = [pi +P2)t.F+{q^) + {pl~p2)^.F^{q^)■, 

0<q^ < {Ml - IVhf (2.100) 

the recoil reads 

w = 2 2_ > 1 2.101) 

2M1M2 - ^ ' 

and the velocity-dependent form-factors are related to the form factors F± according to the relation 



In the limit of infinitely heavy quarks Q1.2, both the elastic and inelastic amplitudes are expressed to a l/mg accuracy 
in terms of the single universal Isgur-Wise (IW) function ^(w) |^ 

h+{u) = h,i{uj) = ^{u) + 0(1/toq), 

h^{Lj)^0{l/mQ). (2.103) 

The Isgur-Wise function can be expanded near zero recoil as follows 

i{Lu) = 1 - p'{uj - 1) + 0{{lo - if). (2.104) 

It is important to stress, that the heavy-quark symmetry predicts the absolute normalization of the transition form 
factor in the heavy-quark limit. 

In addition, the heavy-quark symmetry gives the universal relation for heavy-meson decay constants 

VMp/g = const. (2.105) 

The asymptotic relations ( 2.103 ) and ( ^.105| ) are the zero-order terms of the l/mg-expansion which is calculable 
within the Heavy quark effective theory (HQET) A particular form of the IW function depends on the heavy 
meson wave function. In next Sections we discuss in detail the structure of this expansion for the transition of a heavy 
pseudoscalr meson into p seudo scalar and vector mesons. 

The expressions ( ^.1031 ) and (|2.105| ) mean that the HQ symmetry restricts the possible behavior of the meson wave 
function at large mg. 
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Table ^ gives the results on fp and p^j vs mg at ruq 
function of mg for various values of a. 



0.25 GeV, and Fig.]!!] presents the quantity ^mqjp as the 



Table 2. The decay constants /p of pseudoscalar mesons built up of quarks with the masses mg and rriq and the slope of 
h^i at = 1 calculated from {Mq\Q^)j^Q\Mq) as functions of at = 0.25 GeV . 





a = 0.01 


a = 0.02 


a = 0.04 


a = 0.08 1 


mg, GeV 


fp,MeV 


•2 

Pel 


fp,MeV 


Pel 


fp,MeV 


„2 
Pel 


fp,MeV 


•2 

Pet 


0.25 


151 


0.04 


130 


0.06 


104 


0.08 


80 


0.1 


0.4 


190 


0.25 


160 


0.35 


128 


0.5 


97 


0.65 


1.8 


324 


0.6 


234 


0.65 


163 


0.82 


110 


1.0 


5.2 


308 


0.75 


202 


1.0 


132 


1.05 


85 


1.1 


10 


254 


1.0 


162 


1.05 


102 


1.1 


64 


1.25 


20 


195 


1.0 


122 


1.1 


76 


1.23 


48 


1.45 


40 


143 


1.0 


89 


1.11 


55 


1.25 


34 


1.66 


80 


103 


1.0 


63 


1.11 


39 


1.25 


24 


1.66 



IE 



1 

0.9 
0.8 
0.7 
0.6 
05 
0.4 
0.3 
0.2 
0.1 










a=0.01 




^^^^ a=0.02 




a=0.04 


1 III 


0=0.08 





Fig. 11. The quantity mg fp as the function of mg at 



0.25 GeV. 



In the HQ limit, for a finite binding energy of the meson the heavy meson and the heavy quark masses coincide, 
Mg/mg = 1. So, the value of y/mq/p should be independent of the heavy quark mass. Clearly, the asymptotic 
relations (2.103) and ( p.l05| ) are satisfied if the parameter a of the wave function ( p. 95 ) tends to a constant aoo as 
mg CO. 

Thus, the function a{mQ) has the following behavior: it is equal to 0.02 at mg < 0.5 GeV and tends to a constant 
aoo as mg — > oo. For investigating the B and D mesons and their decays we need the information on a in the region 
mg ^2^ 5 GeV. 

The simplest way is to extract a at mg = 2-^5 GeV from the analysis of /d and Jb as we have done for the light 
mesons. In the absence of the experimental data we refer to the results of other models. As one can see, the decay 
constants fp calculated with a from the range 0.02 < aD,ctB < 0.04 cover the regions 160 MeV < /d < 230 MeV 
and 130 MeV < /s < 200 MeV which include the predictions of most of the models. Hence, the values of ao and 
aB related to the true wave functions of D and B mesons are expected to be inside the interval 0.02 0.04. 

However, there is an attractive possibility to specify aD,B more precisely. Namely, it seems reasonable to assume 
a to be approximately constant in the region mg > 1 2 GeV. There are at least two arguments behind this 
assumption. Firstly, a system consisting of a heavy and a light particles behaves like a quasinonrelativistic system. 
And secondly, there are no visible sources within QCD to yield steep changes of a in this region. Then for the B and 
D mesons one expects ao = as = ckoo- The next step is to estimate aoo- We consider the value aoo = 0.02 to be 
both attractive and reasonable: on the one hand, the same parameter describes all ground-state mesons, and on the 
other hand, one finds for aoo = 0.02 
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5.8 Gel^3/2 



in agreement with the value 0.6 -f- 0.7 estimated in P^ . 

Assuming ao = ctB — ctoo , we can estimate the magnitude of the higher order 1 /toq correct ions which d etermine 
the deviations of the calculated fp and p^j at finite mq from the asymptotic relations ( 2.103 ) and ( 2.105 ). Rather 
strong violation of the HQ symmetry for b— and c— quarks (5 -f- 15% at mq = 5 GeV and 20 ~ 30% at mq — 2 GeV) 
can be observed both in fp and p^j at aoo = 0.02 0.04. 

We shall analyze the transition form factors obtained at au.B — 0.02 and 0.04. If our assumption ao — ctB — 
Oioo — 0.02 does not work properly, the form factor calculations for a = 0.02 and a = 0.04 give an interval which is 
expected to include the true value. Table |l| gives the numerical parameters of the model. 

The results on the axial-vector decay constant fp are shown in Fig.^ and Table ^ Assuming a{mq) = aoo at 
mq > 2 GeV , one can see the asymptotic relation y/mqfp = const to work perfectly at mq > 40— 50 GeV, and 
finds essential corrections to the asymptotic relations at lower mq. For = 0.02 one obtains fo — 234 MeV and 
/b = 202 MeV that confirms the expectation fp, ~ fp . These values for the decay constants correspond to the 
constituent quark Icptonic decay constant g'^ — 1. 

Figures 03, 13 and 04 present the form factors calculated with ap, B — 0.02 and 0.04. 




Fig. 12. The D and B elastic form factor F+. Solid line a — 0.02, dashed line a — 0.04. 
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Fig. 13. The form factor F+ for the B — > tt (left) and D ~> n (right) transitions. Solid - Q:_g(D) — 0.02, dotted - the 
monopole fit, dash-dotted - the dipole fit. Dashed - a^j^,) — 0.04. 
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Fig. 14. (a) The form factor F+{q'^) for B ^ D. Solid - ao = ctB = 0.02, dotted - the monopole fit, dash-dotted - the 
dipole fit. Dashed - qd ~ cib ~ 0.04; (b) The form factor F+ for the D K decay. Solid - ao = 0.02, dotted - the monopole 
fit, dash-dotted - the dipole fit. Long-dashed - the Landau singularity contribution, long-dash-dotted - the non-Landau term. 
Dashed - qo = 0.04. 



Figure |l4| (b) shows the relative magnitude of the various contributions to the form factor F4., separately, as the 
function of for the D ^ K transition. Clearly, the full form factor is a monotoneously-rising function of q^, 
whereas the behaviour of its normal and anomalous parts is rather specific: 

The normal part rises up to some value of q^ where it takes the maximal value, then it goes down rather steeply 
and vanishes at the maximal q^ = (m2 — mi)^ corresponding to 'quark zero recoil point'. 

The anomalous part is identically zero at negative q^ but comes into the game as q^ goes into the positive region. 
It is small for small positive 9^, but increases steeply as q^ approaches its maximal value. 

In other words, at 5^ < the contribution of the non-Landau singularity (anomalous term) is absent, and the 
Landau-type singularity (normal term) determines the form factor; in the region < < (to2 — mi)^ both of them 
are essential; at the point q^ — (7712 — mi)^ the contribution of the Landau singularity vanishes, and the non-Landau 
singularity determines the decay form factor at this 'quark zero recoil' point. 

For heavy-to-heavy meson transitions, a specific relationship between the Landau and the non-Landau contributions 
to the dispersion representation is observed: the normal Landau contribution dominates the form factor at all q^ < 
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(m2 — rrii)^ , whereas the region where the anomalous singularity is essential shrinks to a very narrow vicinity of this 
point. So, effectively the transition form factor is determined by the contribution of the Landau singularity only. 
Thus, the HQ symmetry can be formulated in the language of the analytic properties of the transition form factors 
as the dominance of the Landau singularity in the almost whole kinematical region. 

In the case of the meson decay related to a heavy-to-light quark transition, the anomalous non-Landau contribution 
is important in a broad kinematical region. So the relations suggested by the HQ symmetry would not work properly 
in this case. 

Table U shows parameters of the monopole F+{q'^) — F+(0)/(l — 9^A^mon) the dipole F+{q'^) — F+(0)/(l — 
g^/Mjjp) fits to the calculated transition form factors. The dipole formula approximates the calculated form factors 
with better than 1% accuracy. Parameters of the monopole fit agree very well with the vector meson dominance. 

The K TT transition form factor is well approximate by the linear function = F+(0) -I- aq^, i^+(0) = 0.96, 

a — 1.27 GeV^^ in agreement with the results of 



Table 3. Parameters of the monopole and dipole fits to the F+ form factor. Masses of the lowest vector mesons which are 
expected to dominate the form factors are given in brackets. 





Q-B = O-D = 0.02 


OS — OLD = 0.04 1 


Decay 


F+(0) 


Mmon, GeV 


Mdip, GeV 


F+(0) 


Afd.p, GeV 


B ^ D 


0.73 


5.7 


7.7 


0.68 


7.20 


B 


0.23 


5.2 [5.324] 


6.2 


0.22 


6.08 


D-> K 


0.70 


2.22 [2.11] 


3.0 


0.70 


2.95 


D ^ n 


0.55 


2.1 [2.01] 


2.8 


0.59 


2.68 



E. Discussion 

In this section we investigated form factors of hadron transitions within the relativistic dispersion approach based 
on the constituent quark picture and proposed a formalism for a direct calculation of hadron decay form factors. The 
developed approach was applied to the analysis of the electroweak properties and transitions of pseudoscalar mesons. 

The main results of this chapter are as follows: 

1. We analysed elastic and transition form factors for pseudoscalar mesons at spacelike momentum transfers, q^ < 0, 
and obtained for them double spectral representations in the mass variables si and S2, squares of the invariant 
masses of the initial and final qq pairs, respectively. These representations involve double spectral densities of the 
corresponding triangle Feynman diagrams and wave functions of the participating mesons. 

2. We performed the analytic continuation of the form factors to timelike momentum transfers. We observed the 
appearance of the anomalous cut and, respectively, the appearance of the anomalous contribution to the form factor 
for q^ > 0. The anomalous cut appears due to the motion of the singularities of the Feynman triangle graph from 
the second sheet onto the physical sheet through the normal cut. The normal cut is related to the Landau-Cutkosky 
singularities of the Feynamn graphs, whereas the anomalous cut is related to the non-Landau type singularities which 
come into the game in the region of timelike momentum transfers. It is important ot emphasise that both the normal 
and anomalous contributions to the form factors are expressed in terms of the wave functions of the initial and final 
mesons only in the physical region above the threshold, si > (m2 + m^)^ and S2 > {^i + nT-a)^ ■ 

3. We demonstarted the equivalence of our approach based on spectral representations with the light-cone constituent 
quark model for the description of leptonic decays and for transition form factors between pseudoscalar mesons at 
spacelike momentum transfers. It should be noticed that the dispersion approach has important advantages compared 
to the light-cone approach in the region q^ > 0, where the direct application of the latter is hampered by the 
contribution of pair-creation subprocesses. 
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4. For meson decays related to heavy-to-hcavy quark transitions a dominance of the normal contribution over the 
anomalous contribution for almost all from the decay region except for the very vicinity of the zero-recoil point 
has been observed. This allows a formulation of the heavy-quark symmetry in the language of the analytic properties 
of the decay form factors as the dominance of the normal Landau contribution in the almost whole kinematic region 
of momentum transfers. 

5. Electroweak properties and form factors of pseudoscalar mesons have been analysed using a simple parametrization 
of the meson wave function based on the heavy quark symmetry. We have examined the dependence of the axial- vector 
decay constant on the heavy-quark mass, and found fo — 235 MeV and fs — 200 MeV. 

Analysing the dependence of fp and the heavy meson form factor on the heavy quark mass we have found that the 
violation of the HQ symmetry relations can be expected at the 10-20% level for the b and c-quark masses. 

In the next Chapter we consider the 1/toq expansion of the transition form factors given by the double spectral 
representations and match this expansion to the heavy quark expansion in QCD. In particular, we demonstrate that 
the heavy quark expansion of the P ^ P transition form factors is fully compatible with the 1 /mq expansion in QCD 
in the leading and subleading 1 / mq orders. 
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III. HEAVY QUARK EXPANSION AND UNIVERSAL FORM FACTORS IN THE DISPERSION 

APPROACH 



This Chapter presents a detailed discussion of form factors for weak meson decays within the dispersion approach 
following the analysis of Ref. |^ . 

We calculate the double spectral densities for the form factors describing the transition of a pseudoscalar (P) meson 
to pseudoscalar (P) and vector (V) mesons induced by the vector, axial-vector and tensor currents. These spectral 
densities are given in terms of the soft wave functions of the participating mesons and the double spectral densities 
of the corresponding triangle Feynman graphs. 

We then analyse the spectral representations for the form factors in the two specific cases of the heavy-to-heavy 
and heavy-to-light transitions. 



(i) The heavy-to-heavy transition means that the masses of the initial quark m2 and the final quark mi participating 
in the weak transition are much larger than the confinement scale A 



TO2 > mi, m2 



mg > A. 



(3.1) 



In this case the meson transition can be analysed using the formalism of Heavy quark effective theory (HQET) 
an effective theory obtained from QCD for heavy quarks. Expansions of the form factors in powers of 1/toq 
can be obtained in terms of the universal form factors which appear in each order of 1/toq. These process- 
independent form factors contain the information about the long-distance dynamics in the heavy-quark limit 
and can be calculated only within some nonperturbative approach. 



We studied the form factors of the dispersion approach for the quark masses satisfying the relation ( |3.l| ). We 
perform the l/mq expansion of our spectral representations and require its structure to match to the known 
structure of the l/mq expansion in HQET. 



(ii) The heavy-to-light transition means that 



m2 



niQ > mi ~ A. 



(3.2) 



In this case the explicit structure of the 1/toq expansion cannot be obtained directly from QCD by the existing 
methods, but QCD provides relations between the form factors in the region of large near zero recoil ||]. We 
therefore require the form factors of our dispersion approach to obey these relations. 

The conditions (i) and (ii) allow us to determine the necessary subtraction terms in the spectral representations for 
the form factors. As prompted by the structure of the l/mg expansion of the dispersion form factors, no subtractions 
are necessary for P ^ P transition form factors, but subtractions for some of the P ^ V form factors are necessary. 



Section III A gives the definitions of all the necessary form factors for P ^ P and P ^ V weak transitions. The 
structure of the 1 /mg expansion of the transition form factor induced by the vector and axial- vector current has been 
calculated within HQET in |q| . We apply the formalism of HQET to transitions induced by the tensor current and 



report the 1 /mg expansion of the form factors hg^ , hg 



hgg , and /i.. 



(the definitions are given in the next section) in 



the leading (LO) and next-to-leading (NLO) orders in l/mg. In particular, we find the l/mg correction to the value 
of the form factor hg^ at zero recoil to vanish exactly as for the form factors hf^ and hf (the Luke theorem). 

We later use the l/mg expansion of the weak form factors from HQET as a benchmark for testing the l/mg 
expansion of the form factors of our dispersion approach. 



In section [II B we calculate the double spectral densities for form factors and give results for the necessary sub- 



traction terms which are explained later in the sections III C and HI D 



In section III C we consider the case of the heavy-to-heavy meson transition and perform the expansion of the 
spectral representations for the transition form factors in the LO and NLO in 1 /mg . 

The spectral representations without subtractions are found to agree with HQET in the LO for all P ^ P and P 
V form factors. The next-to- leading order analysis shows the necessity of subtractions in the spectral representations 
for some of the P V form factors in order to match to the structure of the HQET. The matching condition provides 
constraints on the subtraction terms. 
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Assuming the strong peaking of the meson soft wave functions in terms of the relative quark momenta with a width 
of order of the confinement scale, we calculate the Isgur-Wise function and the NLO universal form factors in terms 
of the wave function of the infinitely heavy meson. Our dispersion approach leads to the following relations for the 
universal NLO form factors: 

X2(w) = 0; 
X3(w) = 0; 

> 0, C3(l) = (^)/3, (3.3) 
where (z) is an average kinetic energy of the light quark in the heavy meson rest frame. 



In section III D| we discuss heavy-to-light quark transitions and consider the 1 / mg expansion for the wave functions 



and the form factors in this case. We then consider the heavy-to-light meson transitions in which case a small 
parameter Aqcd/'t^Q emerges and analyse the form factors in the leading KqcD/'niQ order. Requiring the fulfillment 
of the Isgur-Wise relatins for the heavy-to-light transitions Q further constrains the subtraction terms providing 
explicit spectral representations with subtractions for the form factors of interest. 



Section III E illustrates the main results with numerical estimates and evaluate the universal form factors for various 



quark model parameters. 



A. Meson transition amplitudes and heavy-quark expansion in QCD 



The amplitudes of meson decays induced by the quark transition q2 qi through the vector = qij^q2, axial- 
vector Aa — gi7u7592, tensor T^i, — qiau,vq2, and pseudo-tensor T^^, = qiau,vlbQ_2 currents have the following structure 

(l/(M2,p2,e)|V;(0)|P(Mi,pi)) - 25(g2)e^,„^e*>^pf, 

(M2, P2, e)\A^ (0) |P(Mi , pi )) = [ f{q^)g^,c. + ai {q^)Plo.Pl^. + ^2 (g'>iaP2^ ] , 
(P(M2,P2)|T^.(0)|P(Mi,pi)) = -2^s(q2) (pi^ps. - Pi.P2m), 
{V{M2,P2,e)\T^M\P{Mi,Pi)) = [.9i((z')ep.a/3p'^ +52(g')eM-/3P"' 

+go{q^)piaef,vpyPiPl], 

{P{M2,P2)\TlM\P{MuPi)) = siq^^t^ucpP^q^. 
{V{M2,P2.e)\T^^M\P{MuPi)) = 9l{q^)KPl^. ~ e>.) + <?2(9')(e>2M - 

+9o{q^){<^*Pl){PluP2t, - Plt,P2u), (3.4) 

with q—pi — P2, P — Pi + P2- Also the following linear combinations of the form factors will be used 

.f± = ^(/l±/2), 

1 , , 

a± = -(ai ± a2), 

5± = ^(5i±52)- (3.5) 
The matrix element of the penguin operator relevant for rare decays has the following structure 
(F(M2,p2,e)k-icTM-9"(l +75)<Z2|P(Afi,Pi)) - -ie,,.o/5e*''p>^2g+ 

-(g^, • qP - q.P^)e*'' {^g+ - ^.g_ j + {e*q) (^g^ - ^P^ (^.9_ - \qP 9^ ■ (3.6) 

We use the following conventions: 
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Accordingly, 

Spi^'ri"-/"!^) = 4*6'^''"^ ^^,7^ = -'-e^^a.pT"''- (3.8) 

The relativistic-invariant form factors contain the dynamical information on the process and should be calculated 
within a nonperturbative approach for any particular initial and final mesons. 

For analysing the transition in the case when both the parent and the daughter quarks inducing the meson transition 
are heavy, i.e. mi ~ TO2 ^ J^qcd it is convenient to introduce a new dimensionlcss variable 



and velocity-dependent form factors h{uj) connected with 4-velocities and not 4-momenta as in ( p.4D in the following 
way 



{P{M2,P2)KmPiMi,pi)) = y/MlM2[hf^{io){vl+V2)^. + hf_{u){vi - ^2^], 
(l^(M2,P2,e)|V;(0)|P(Mi,pi)) - ^MiMa /^3MeM-a/3e*"<«2^, 
(T/(M2,P2,e)|A^(0)|P(Mi,Pi)) = ^e*" .J M^Ah [hf{uj){l+u)g^^ 

-ha^{uj)viaVif_, 

"■02 

{uj)vi 

(P(M2,P2)|r^.(0)|P(Afi,pi)) = -2Z \/MiM2 h,{uj) {VI^V2. - Vi,V2^), 

(T/(M2,P2,e)|T^.(0)|P(Mi,pi)) = ie*''^AhM2 [hg4u;)e^,^p{vi + V2f 

+ hg_ {uj)ef_,^a/3{vi - V2)'^ 

+ hgg {Lu)viaefj,^pjv'^v^ ] . 
(P(M2,P2)|T' (0)|P(Mi,pi)) = 2 VM1M2 /is(c^) eM-/3<«2 . 



(T/(M2,P2,e)|T^'.(0)|P(Mi,pi)) = ^MiMa [/i<,+ + ^2)^ - +^^2).) 



+hg_ (w)(e*(wi - W2)m - - ^2)^) 
-|-/igo(w)(e*ui)(wi,.U2^ - vif_,V2„)] ■ (3.10) 



These form factors are related to the form factors introduced by the relations (3^) as follows 



2VM1M2 

^ ~ 2VMiM2^^' 

1 M2, 



1 



a2 



/= v/MiM2(l+c^)/i/, 
1 

2VM1M2' 
M2 
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The form factors h are convenient quantities as in the leading 1 /mg order all of them are expressed through a single 
universal function of the dimensionless variable lu - the Isgur-Wise function. A consistent heavy-quark expansion 
of the form factors, i.e. expansion in inverse powers of the heavy-quark mass, can be constructed within the Heavy 
Quark Effective Theory based on QCD with heavy quarks. 

The general structure of the l/mg expansion of the heavy quark form factors in QCD for the meson transition 
Ml —^ M2 induced by heavy quark transition TO2 ^ "m-i have the form (omitting corrections 0{as, as/mq, l/mg): 

\mi ni2 J \ 2 

hg= e+ — + — P2 + — 

\ mi 711,2 / 2 mi m2 

{ I 1 \ uj-lX^ 1 If t^-U 

\mi 7712 J LO + i 2, mi 7)7.2 \ 07+1 

K,= ^(-Ae + 2(c^ + l)X2-6), 

777,1 + 1 

foj-l I 1\A If 

\a; + l 777i ?772 / 2 mi \ u + 1 J 7712 

hs^ e+ ( — + — ) + 

\777 l 7772 / \ 2 

= - C - — Pi - — P2, 
7772 n^l 

\77il 7772/ 2 7772 

V= -L(«±f + 2,«). ,3.12, 

7771 \ ti; + 1 / 

In the leading 1 /mg order (LO) all the form factors are represented through the single universal Isgur-Wise function 
£,{ijj) , whereas in the next-to-leading order (NLO) the 4 new form factors pi,p2,X2i a-nd ^3 appear. The universal 
form factors are functions of a single variable lo. 

The form factor ^3 originates from the expansion of the transition quark current, and the form factors pi,p2,X2 
are connected with the nontrivial relationship between the mesonic states in the full and the effective theory. The 
universal form factors satisfy the conditions 

Pi(l)-0, 

P2(l)-0, (3.13) 

whereas no constraints on ^3 and X2 are imposed by the heavy quark symmetry. As found by Luke |^], the l/mq 
corrections to the form factors hf and hf^ vanish due to kinematical or dynamical reasons: 

hf{l) = l + 0{l/ml), 
hu{l) = l + Oil/ml). (3.14) 

The same is found to be true also for the form factor hg^ : namely, 

V(l) = -l + 0(l/m|). (3.15) 



The parameter A in (3.12) comes from the l/mq expansion of the mass of a meson consisting of the heavy quark 
and light degrees of freedom 

A/q = 777q+A + C»(1/777q). (3.16) 

In our notations for heavy quarks and mesons, this gives 
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Ml = 7712 + A + . . . , 
M2 = 7771 + A + . . . , 

for the parent and daughter particles, respectively. 

It is straightforward to derive the following useful relations 

Ml + M2 7772 + ^771 



y/MiM2 \/rnim2 
Ml — M2 m.2 — mi 



VM1M2 y/mim2 



VM1M2 = y/mrm^ 
where the dots denote high er ord er ter ms. 



(1 


1 

+ — 


\ 777 1 


171,2 






V 7771 


^ 7772 


1 


^) 


7771 


7712/ 



/ 77l2 — 7r7l 



A 
2" 



Using the relations ( 3.11 ) and ( 3.18 ), we obtain for the form factors (3.4) the following expansions 

7771 



/i 

/2 



\/mim2 

7772 

1 



91 

92 

9 

ai = 
a2 = 



2^777 x 7772 

mi 



^/mim2 

7772 

^/m^m^ 
1 

2^777 1 7772 
1 



— (Pl+6) + — (^1-6) 
7771 ^772 

C + — + — (pi+6) 

rr7i 7r72 

e + — (pi-6) + — (pi-6) 

mi m2 



7771 m2 

e + — P2 + — (Pl + 6) 

7771 "7-2 

? H P2 H [Pl~ 6) 

7771 ^772 

1 1 



^mim2 7772 + 

1 



- [-Ae + 2(L^ + l)x2-6] , 



^/mim2 



^ + — (Pi - 6) + — (Pi - 6) - —^77^ 
7r7i 7r72 7r7i a; + 1 



-^3 



/ = ^/mim,2{uJ + 1) 



2x2 1 ^ 
mi mi Lu + 1' 



,1,— H \+ — P2 + — pi 

CiJ + 1 V 777 1 777 2 / 777 1 777 2 



UJ-l 
UJ+1 



6 



50 



(77717772)^^/^ 



Ae + 6 

w + 1 



2X2 



(3.17) 



For the following analysis it is worth noting that the behavior of the combination 2piP2 ■ g — mi ■ §2 and / 
NLO coincide, namely 

/ ~ 2pip2 ■ g-mi- g2. 
It is also convenient to introduce the form factor a'2 such that 

a'2 = 0,2 + 2s. 

In what follows we need the expansions of the following linear combinations of the form factors a'2 and ai 

1 



017)7,2 — 027)7,1 



Ol7772 + 027)71 = 



^/mim2 

1 

y/777 l 777 2 



[4X2 - Cs] 

2Ac 

LU+1 



UJ- 1 
'cj + 1 



(3.18) 

(3.19) 
(3.20) 
(3.21) 
(3.22) 
(3.23) 
(3.24) 
(3.25) 

(3.26) 

(3.27) 

(3.28) 
in LO and 

(3.29) 
(3.30) 

(3.31) 
(3.32) 



36 



B. Transition form factors in the dispersion approach 



The results presented in the previous section are strict consequences of QCD in the heavy-quark hmit. The universal 
form factors ^ and pi,P2,X2,C3 can be calculated within a nonperturbative dynamical approach. 

We study the form factors within the dispersion approach based of the constituent quark picture which turns out 
to be an efficient framework for describing meson decays. We start with < and represent the form factors as 
double spectral representations in the invariant masses of the initial and final qq pairs. The form factors at > 
are derived by performing the analytical continuation. 

The transition of the initial meson q{m2)q{'m^) with the mass Mi to the final meson q{Tni)q(jn^) with the mass M2 
induced by the quark transition TO2 — > mi through the current q{mi)J ^q{m2) is described by the triangle diagram of 
Fig.|l^. As explained in the previous section, for constructing the double spectral representation we must consider a 
double-cut graph where all intermediate particles go on mass shell but the initial and final mesons have the off-shell 
momenta pi and p2 such that = si and p| = S2 with {pi — ^2)^ — q^ kept fixed. 




Fig. 15. One-loop graph for the meson transition. 



For the transition B —f D,D* we have 7712 = mi,, mi = mc, and = mu- The constituent quark structure of 
the initial and final mesons is given in terms of the vertices Fi and F2, respectively. The vertex describing the initial 
pseudoscalar meson has the spinorial structure 

Fi = nsGi/^/K. (3.33) 
The vertex for the final pseudoscalar state reads 

F2 = 77562/ (3.34) 
The final vector meson is described by the vertex 

r2^^[AJ^ + Bikl-k3)^,]G2/^/N',, (3.35) 

where for an S'-wave vector meson 



A=-l, B = l/{^ + mi+m3). (3.36) 

As explained above, the double spectral densities / of the form factors are obtained by calculating the relevant traces 
and isolating the Lorentz structures depending on pi and p2. The invariant factors of such Lorentz structures provide 
the double spectral densities / corresponding to contributions of the two-particle singularities in the Feynman graph. 
Let us point out that this procedure allows one to obtain double spectral densities, whereas subtraction terms should 
be determined independently. We determine the subtraction terms from matching the l/mg expanded form factors 
of the quark model to the heavy-quark expansion in QCD. 

Recall that at < the spectral representations of the form factors have the form 

fi^l^^ = TT^ [ ds2(p2{s2) I dsnpi{si) .fi^'^^'^ L ' 
where the wave function ipi{si) = Gi{si)/{si — Mf) and 
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+ , 2\ S2(»7ii + ml — (J ) + g (toi + TO?) — (toi — TO2)(mi — mi) 
■^i (*2,g ) = 7^—2 ^ 



± 



2m| 

A^/^(s2, m§, m^)A^/^(g^, m^, m|) 



2m^ 



Here A(si, S2, g^) — (si + S2 — 9^)^ — 4siS2 is the triangle function. 

Calculating the necessary traces and isolating the coefficients of the relevant Lorentz structures, we come to the 
following expressions for the double spectral densities fi{si, S2, q^) of the form factors (27): 

s = 2 [mia2 + m2ai + m3(l — ai — 02)], (3.38) 

/i — 2mis + 4a2[s2 — (mi — ms)^] — 2m3S, (3.39) 

/2 = 2m2S + 4ai[si - (m2 - ms)^] - 2m3S, (3.40) 

g = -AS - ABf], (3.41) 

gi = Ah-8f3 + 8B{mi+m3)P, (3.42) 

92 = Ah + 8B{m2 - m3)f], (3.43) 

0.20 = -2S + 4BC2ai + ai2Co, (3.44) 

aiD = -4A(2m3 + SCi)ai + aiiCo, (3.45) 

m2 , 2 2n '^l / 2 2\ ^3 / 2 2 2\i 

Jd = -4A[mim2m3 + — (s2 - m^ - m3) + — (si - "^2 - m3) ^(q - m^ - ma)] 

+CoP, (3.46) 

5oD = -8Aai2 - 8B [-?7i3ai + (m3 - m2)aii + (m3 + mi)ai2], (3.47) 



where 



Co = ~8A(m2 - mg) + 4BC3, (3.48) 

Ci = S2 - (mi + m3)^ (3.49) 

C2 = Sl - (to2 - TO3)^ (3.50) 

C3=q^- (mi + m2)2 - Ci - C2. (3.51) 



The coefficients ai, 012, 0:11, CH2, and /3 are defined according to the relations 

dkidk2dkz5 [pi — k2 — k3)6{p2 — k^ — ki)5{k\ — m\)6{kl — m\)5{k2 — m1)k3^ 
= (aiPiM + a2p2ti) A(si, S2, 9^), 

^ J dkidk2dk3S{pi - k2 - k3)6{p2 - /s3 - ki)5{kf - m,1)S{k2 - m|)J(fc| - ml)k3^k3i, 

= {ailPlfiPli^ + ai2(PlMP2i. + Pli^P2p) + a22P2MP2<. + /JffA"') A(si, S2, 9^), 



given 


by Eq. (^.79|). Exphcitly, they 


read 


ai = 


[{si + S2 - q^){s2 - m\ + ml) 


- 2s2(si - m2 + mg)] /A(si,S2,g 


a2 — 


[{si + S2 - q^){si ~ml + ml) 


- 2si(s2 - m^ + mg)] /A(si,S2,g 




— [2ml ~ ~ "^2 + "^3) ^ 


a2(s2 - m^ + TO3)] , 


ail = 


a\ +4/3s2/A(si,S2,9^), 




"12 = 


Q!lQ!2 - 2/3(si + S2 - (7^)/A(si, 


S2,g^), 



(3.52) 
(3.53) 

(3.54) 

(3.55) 
(3.56) 



In Eqs ( 3.44 - 3.47 ) we label with a subscript 'D' the true double spectral densities obtained from the Feynman graph 
for those form factors which require subtractions. We fix this subtraction procedure by requiring the l/mg expansion 
of the form factors to have a proper form in accordance with QCD in leading and next -to-leading orders for the case 
of a meson transition caused by heavy-to-heavy quark transition. It is convenient to introduce for such form factors 
the modified double spectral densities which include properly defined subtraction terms as follows 
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f = fD + - 5i) + (M| - S2)]~g, 

J^2 J 2' 



02 = 0'2D 



9a = 5oD + 



{uj + l)m2 V 
(<Ii + l)mi V ^sT ^/sj / 2' 

{u + l)mim2 V \/s2 / 2 



(3.57) 
(3.58) 

(3.59) 

(3.60) 



As the analytical co ntinu ation to the timelike region is performed, in addition to the normal contribution which 
is just the expression ( |3.37 ) taken at > 0, anomalous contribution emerges. As we have shown in the previous 
section, the origin of the anomalous contribution is connected with the motion of the zero of the triangle function 



X{si, S2,q^) located at sf = (y^ — )^ from the unphysical sheet at < onto the physical sheet at > 
through the normal cut between the points sj~(s2,(?^) and sf{s2,q'^)- Pinching of the points and sj~ occurs at S2 
given by Eq. (2.88). The spectral representation of the form factor at > takes the form 



IGtt^ J J X^'^(Si,S2,q^) 



(3.61) 



One should also take into account that if the double spectral density / is singular at the point sf , which is the case for 
the form factors describing transitions to vector meson, then a properly defined anomalous part contains subtraction 
terms which appear as the contribution of a small circle around s{*. Obviously, these subtractions have quite different 
nature from subtractions in the spectral representations at < performed to match the structure of the heavy 
quark expansion in QCD. The corresponding expression was derived and can be found in |p6|| . 

The normal contribution dominates the form factor at small timelike and vanishes as q^ = (rn2 — m-i)^ while the 
anomalous contribution is negligible at small q^ and steeply rises as q^ —^ {1112 — wi)^- 

Both for pseudoscalar and vector mesons with the mass M built up of the constituent quarks mq and lUq, the 
function ip is normalized according to the relation 



87r2 



Ai/2(. 



(mg - mqf] = 1. 



(3.62) 



This equation corresponds to the normalization of the elastic charge form factor at = 0. 

Recall the spectral representation for the leptonic decay constant of the pseudoscalar meson fp 



ds ip{s) 



Stt^s s 



(3.63) 



In the subseque nt section s we analyse th e form factors given by the dispersion representation (3.61) with the 
spectral densities ( ^.38 -3.43) and (3.57-3.6C) and demonstrate them to have the following properties in the case of 
a heavy parent meson: for the transition induced by the heavy-to-heavy quark transition they satisfy the LO and 
NLO relations |6| of the 1 /mg expansion in accordance with QCD provided the functions (pi are localized near the qq 
threshold with the width of order Aqcd- For the meson decay induced by the heavy-to-light quark transition they 
satisfy the LO relations between the form factors of the vector, axial- vector and tensor currents |q|. 



C. Heavy-quark expansion in the dispersion approach for heavy-to- heavy transitions 

In this section we consider the form factors of the dispersion quark model in the case when both rn2 and toi are 
large. We calculate the universal form factors and demonstrate that requiring the structure of the 1 / mg expansion in 
the quark model to be consistent with the structure of such expansion in QCD allows us to determine the subtraction 
terms. 
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Soft wave function and normalization condition 



First, we need to specify the properties of the soft wave function of a heavy meson. A basic property of such soft 
wave function (p{s,mQ,ing, A) is a strong peaking near the qq threshold. For elaborating the I/uiq expansion, it is 
convenient to formulate such peaking in terms of the variable z such that s = (mg + + z)'^ (hereafter we denote the 
mass of the light quark as TO3). The region above the qq threshold which contributes to the spectral representation 
corresponds to z > 0. A localization of the soft wave function in terms of z means that the wave function is nonzero 
as 2; < Aqcd- In the heavy meson case mg ^ ~ z ~ A. Let us notice that for a heavy meson the localization in 
terms of z is equivalent to the localization in terms of the relative momentum in the meson rest frame 



z(z + 2m3) + 0(l/mg). 



(3.64) 



The normalization condition ( p.62| ) which is a consequence of the vector current conservation in the full theory provides 
an (infinite) chain of relations in the effective theory. Namely, expanding the soft wave function in 1 /mg as follows 



:0o(2;,TO3, A) 



V3(s,mg,m3, A) = - 

we come to the normalization condition in the form 

dz(?!)g(z)V^(z + 2to3)3/2 



Amr. 



Xi(z,m3,A) + 0(l/m|) 



1 + 7; Xi{z) - 7, ^ 



2mr 



2mr 



(3.65) 



(3.66) 



This exact relation is equivalent to an infinite chain of equations in different l/mg orders. Lowest order relations 
take the form 



/ dz 02(^)yj(z + 2m3)3/2 = 1, 

/ dz 4'l{z)Mz + 2to3)'/'xi(2) = 1. 



(3.67) 
(3.68) 



The variables u and u) 



In the description of the transition processes the dispersion formulation of the quark model has the following 
feature: since the underlying process is the quark transition, the relevant kinematical variable for the description of 
the dynamics of the process is the quark recoil Cu which is defined as follows 

q^ = (m2 — mi)'^ — 2niim2{(jj — 1). (3.69) 

The relationship between lo and uj is given by the condition that the spectral representation for the form factor is 
written at fixed value of q^ and hence 

g2 ^ (^Mi - - 2MiM2{uJ - 1) = (to2 - mif - 2mim2{u} - 1). (3.70) 

In the case of heavy particle transitions these quantities are related to each other as 

w = w + A( — + — ) (tj-l) + 0(l/m^). (3.71) 
\mi m2 J 

We shall obtain the representations of the form factors as functions of the variable lo. The variables lu and (D are 
different by the terms of order l/mg at nonzero recoil. On the other hand, the quark and meson zero recoil points 
coincide with 1 /wg accuracy. This means that in the analyses of the 1 / mg expansion at nonzero recoil the difference 
between the ld and Co might be ignored in the Isgur-Wise function, but gives nontrivial contribution to the NLO form 
factors. At the same time, at zero recoil the difference between uo and Co might be neglected both in the leading and 
next-to-leading orders. Namely, the quark model provides the expansion of the form factor in the following form 

h = ho{Cj) + —h\^\co) + —h\'^\co) + ... (3.72) 
nil m2 

= ho{io) + h'^{uo)K (— + + —h^i\Q) + —h''i\cj) + .... (3.73) 
\mi m2 / mi m2 
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As we shall see later, among the NLO form factors only pi^2 are affected by the the difference between lo and w 
whereas ^,■^31X2 are not. 

Relative magnitudes of the normal and the anomalous contributions 



We are going now to demonstrate that the anomalous contribution comes into the game only in close vicinity of 
the zero recoil point whereas beyond this region is negligible. 

Let us study the behavior of the anomalous contribution in the region 



1 



-(2+e) 



(3.74) 



Introducing the variables Zi and 22 such that si = (7712 + m^, + zi)^ and S2 = {nii + 7713 + wc find that the 
magnitude of the anomalous contribution is controlled by the value of ^2 1*^) such that S2 = {mi + 7713 + which 
is the lower boundary of the Z2 integration. If Z2(fi) becomes large, i.e. of the order 777q, the anomalouis contribution 
is suppressed by the fall-down of the soft wave function. This suppression is at least stronger than l/r7^g. This means 
that the anomalous contribution is nonvanishing only if 



In the region uj — \ 



7)7,17712 



7T7l7772Ct; — 77lJ 



\/ m\ + 77l2 — 27771777,20; 



7771 +0(7773) ~ A. 



(3.75) 



777g' ■ , one finds Z2(^) 



Hence the anomalous contribution comes actually into the 



game only in the 0(1/?77q) vicinity of the zero recoil point but otherwise might be neglected. On the other hand, at 
the quark zero recoil point Q = 1, the normal contribution vanishes and the form factor is given by the anomalous 
contribution. 

We shall calculate the form factors in the region u) — \ = 0{1) where only the normal contribution should be taken 
into account in leading and subleading orders. 



The LO analysis 



To perform the LO analysis of the form factors let us start with the integration measure. With l/mq accuracy it 
can be represented in the form 



1 dsids20 (s2 > (7771 + 7773)^ 



(Sl <Sl<S{ 



167r2 Al/2(si,S2,92) 

and the expression for zi reads 



-^dz2^/z2{z2 + 2m3)^9iz2)e{rf < 1), 



Zi = Z2U! + 7773(0; - 1) + T]y^ Z2{Z2 + 27773) \/o;2 - 1 + 0(l/r77Q). 



(3.76) 



(3.77) 



Let us point out that the LO integration measure is sy m metr ic in zi and Z2. 

Next, we need expanding the spectral densities ( 3.3g )-( 3.47 ). To this end we must take into account that under the 
integral sign zi and Z2 are localized i n the r egion z < A due to the soft wave functions 4>{z). 

In LO the kinematical coefficients ( ^.48 )-( p.52| ) in the region a; — 1 = 0(1) simplify to 



A(si, S2, 9^) = ^m\ml{uj'^ - 1), 
1 



ai 



a2 



13^ 



7772(0; + 1) 
1 

7771 (01 + 1) 
1 



7773 + Z2 + 

m3 + zi + 



Z2 - Zl 

Co- 1 

Zl - Z2 

Q - 1 



an = ai 



o; + 1 

/3 



(7773 + Zl)(7r73 + Z2) + 



{zi - Z2) 



21 



m\{ijj'^ — 1) ' 



0-12 = Q;iQ!2 



(3u) 



7(7i 7772(o;2 — 1) 



(3.78) 
(3.79) 

(3.80) 

(3.81) 

(3.82) 

(3.83) 
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B = 



1 



2mi 

Co = -Am2{u)- 1), 

Ci = 2miZ2, 

C2 = 27712(2:1 + 27713). 



(3.84) 

(3.85) 
(3.86) 
(3.87) 



One finds the LO behavior of the form factor densities (3.38-3.44) is determined by the term proportional to s. The 
latter reads in the LO 



S ~ 2 771, 



zi + Z2 + 2m3\ 



w + 1 



The LO expression for / takes the form 



/l5 = (CJ + 1) 7773 + 



Zl + Z2 + 27)75 
(I) + 1 



(3.88) 



(3.89) 



The spectral densities ai and g^o vanish in the leading order. 

Hence the LO relations (3.1E-3.28) are fulfilled with the Isgur-Wise (IW) function 



= j d.Z2(f>o{z2)\/ Z2{z2 + 2m3) J ^(l)o{zi) ^7773 + — ^ 



27773 + Zl + Z2 



(3.90) 



I n (|3.90 ) we used the equality of oj and uj with l/mq accuracy. The normalization condition of the LO wave functions 
(|3.67) yields ^(1) = 1. For the slope of the IW function at zero recoil, — — f'(l), one finds 



P^^l f dz^{z + 2m:,f'^{dp'o{z)fz{z + 27713) 



(3.91) 



Let us point out that the subtraction terms in the spectral densities do not contribute in the LO relations. As we 
shall see later, they are important in the NLO analysis. 

The NLO analysis of the form factors s, f\, f2, gi, g2, dnd g 



First, let us concentrate on the NLO relations (3.19-3.24). It is convenient to analyze the linear combinations of 
the form factors which do not contain the LO contribution. These combinations are 



1 



1 



9~3 

91 + fi 

92 + /2 

/l - 27771S 



[P2 - Pl+ 6] 



2^ 777 l7n.2 7771 

, [-P2 +P1+ £.3] , 

VTO17772 

7772 1 



7771 ^777 1 7772 

26 



[P2 - Pl+^3 



The spectral densities of the form factor combinations in the l.h.s. of cqs. ( 3.92-3.94) read 









9 


~ s 


7771 


91 - 


f /i 


= -4/3, 


92 ~ 


f /2 


= 4^/3, 

7771 



Comparison with the eqs. ( 3.92 - 3.94 ) leads to the relation 

Pl{uj) = P2{uj), 



(3.92) 
(3.93) 
(3.94) 
(3.95) 

(3.96) 

(3.97) 
(3.98) 

(3.99) 
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For the form factor ^3 we come to the representation 



dz2 00 (^2 ) \/ 22 (2:2 + 21713) 
2 



Y0o(zi)2 



— (to3 + Zi)(to3 + Z2) + 



1 



(3.100) 



with Zl given by the expression (3.77). In (3.100) we have neglected the 0(1/toq) difference between uj and oj. 
On the other hand, the equation (3.95) yields the representation for the form factor ^3 in a different form 



6(w) = / dz2cf)o{z2)\/ Z2{Z2 + 2TO3) / ^^0(21) 



Z2 + 2r7i3 



Z2 + TO3 



(zi - Z2) 
LO - 1 



m3 
2 



m3 



2m3 + Zl + Z2 
1 + 



(3.101) 



One can check that for the soft wave functions providing convergency of the integrals and nonsingular at z = the 
representations (3.100) and ( |3.10l| ) are equivalent. At zero recoil one finds 



^3(1) = / dz2 ^2{Z2 + 2m3)3/2 



(3.102) 



where (z) is an average kinetic energy of the light constituent quark inside an (infinitely) heavy meson in its rest 
frame. It is worth noting that the function ^3 is positive for all uj. 

The universal form factor pi = p2 can be found from the expansion of s{ijj) with a NLO accuracy. In this case 
the NLO terms in the XjmQ expansions of the integration measure, the wave function, and the spectral density s 
contribute. Namely, we can write for the form factor the expression 



hsiuj) 



TO2 



m3 



4m 
1 1 
mi m2 / 



~f-\ , I 

s{uj) H 1 

mi m2 



mg 



^M^i) 1 + 7— xi(2i) 00(^2) 1 + 7— xi(^2) 



4m2 



CW(c^) + 



+(— + — ) K^'' (^) + ?'(^)A(c. - 1)] 
mi m2 



(3.103) 



and hence pi(w) = ^^-^^(w) + ^'{u!)A{lu — 1) — ^3{uj) — ^^{io)- We do not present exphcit expression for pi. However, 
let us consider pi at zero recoil. The analysis of the anomalous contribution at uj = 1 gives 

hJoj = l) = 1- 



1 1 

mi m2 

2 



3 



(3.104) 



Using the relations hs{Lu = 1) = /is(a; = 1) + ©(l/mg) and ^^3(1) — (z) /'S we obtain pi{l) = just as required by 
HQET. 
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The NLO analysis of the form factor f 



First, let us demonstrate that the dispersion representation of the form factor / requires subtraction. To this end 
consider the anomalous contribution at a) = 1. For the form factor fjj constructed from the spectral density fjj 
through the double dispersion representation without subtractions we find |^ 



= 2v/MiM2 



l + l ( — + — ) ((z>+m3-A) 
2 \ mi m2 ' 



(3.105) 



This value contradicts the Luke theorem which requires the l/mg corrections to vanish at zero recoil. Let us 
demonstrate that the form factor / constructed from the spectral density with the included subtraction term 



f ^ fD + iMi-si+Mi-s2)g 
= fD + (2piP2 - 2piP2)g, 



(3.106) 



satisfies the NLO relation (3.27). Here 2pip2 — si + S2 — q ■ 



We have noticed above that to LO and NLO the expansion of the form factor / and the combination 2piP2- g — mi- g2 
coincide ( 3.29] ). Hence checking the NLO relations for the form factor / is equivalent to checking with the NLO 
accuracy the relation 



fo — /2pi 



P2-g 



mi • 32- 



The spectral density of the r.h.s. of eq. (3.107) can be written as 



2piP2 ■ g — mi ■ g2 — 2mim2{uj + l)s + [2mi(z2 + m^) + 2m2(2:i + m^)]s 
— 8(mi + m2)/3 — 4m2(<D — l)/3. 



(3.107) 



(3.108) 



For checking the expression (3.107) in NLO we need the expansion of the spectral density s in LO and NLO which 
has the structure 



1 

r 



miQ!2 + m2Q:i + m3(l - ai - 012) 

zi + Z2 + 2to3 ^ Si ^ S2 
Co + 1 mi m2 



ma 



with 



1 Zi — Z2 (zi + Z2 + 2m^)^ 
Si ^ -22(2:2 + 2m3) - zims - {z2 + 2mj,)— — — . 

2 bJ — 1 UJ +1 



(3.109) 



(3.110) 



and 5*2 is obtained from Si by replacing zi and Z2. The spectral density reads 

/£, = 4mim2((I' + 1) (ma + fl_L^H_L?!!j!. j -|- 2m222(22 + 2m3) + 2miZi(zi + 2m3) - 4m2(tl' - l)/3. 



(3.111) 



Notice that for checking the NLO relation ( [3.107] ) between the form factors we do not need explicit expression for the 
integration measure in the NLO: in the LO the spectr al densities are equal and he nce th e NLO contributions from 
the integration measure into both sides of the eq (3.107) are equal too. Finally, Eq. (3.107) is satisfied if the following 
relation is valid 



/ dz2 \/ Z2{Z2 + 2m3)0o(22) / ^0o(2l) 



-Z2(z2 + 2m3) + 2(zi + mg) ( + ^ ) - 4/3 + 2(^5 + \)Si 

+ 1 



(3.112) 



Hereafter we denote as ao a form factor reconstructed from its double discontinuity through dispersion representation without 
subtractions. 
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with zi given by (3.77). One can check this relation to be true for a ny fun ction (j)o{z) regular at 2; = 0. Hence, the 
form factor / calculated with the subtracted double spectral density ( 3.106 ) satisfies the HQET relations in LO and 
NLO at all uj. 

Strictly speaking, the NLO analysis does not allow us to uniquely specify the subtraction term: namely, any spectral 
density of the form 

f^fD + (Ml - 51 + Af2 - S2)pf (3.113) 
has a proper NLO behavior in accordance with ( p. 27 ) provided the spectral density p j behaves in the LO as 

zi + Z2 + ^rtiz 

UJ+l 



As we demonstrate in the next section, the analysis of the heavy-to-light LO relations requires the identification 
P/ = 5- 

Let us notice that the form factor obtained within the light-cone approach [n4j can be represented as dispersion 
representation with the spectral density 



he = ^/d + M2 

\/S2 



Si - S2 - Ml - M| - g2 ai + 02 



2M2 



At zero recoil one finds 



/ic(l) = 2v'MiM2 
that contradicts the Luke theorem. 



^-\(- -]{{z)+m,-K) + 

2 \ mi 7712 



(3.114) 



(3.115) 



The NLO analysis of the form factors ai.2 and X2 



The double spectral densities of the form factors ai and 03 = 02 — 2s in the LO read 

= 4 [(z2 + 27773)01 - m2aii(cj - 1)] , 



0-20 = 4 



"^2 / 

(Z2 + 27773)01 - m2ai2[uj - 1) 

7771 



(3.116) 
(3.117) 



The quantity 02 is more convenient than 02 for calculations because the LO term of the heavy quark expansion of 02 
is zero. 



The unsubtracted spectral representation for ai(a>)77i2 ~ a2(w)77ii in combination with eq (3.31) gives 

1 

X2{^) = ^si^^) ~ ^ dz2(t)o{z2)\/z2{z2 + 27713) / -^Mzi) [aiTO2 - a2mi] . 

-1 



Substituting the representation ( 3.100| ) for ^3 we find 

X2(w) = 0. 



(3.118) 



(3.119) 



Let us now consider the linear combination ai{uj)m2 + 02(^)7771. As a first step, show that the unsubtracted 
dispersion representation is not compatible with HQET. To this end calculate the unsubtracted form factors oid and 



a'2D at zero recoil: 



aiD(l) = 
4d(1) = 



1 



1 



^77717772 777 2 
1 1 



^777 1 7772 777 1 



2 ms' 

— (z) H 

3^ ' 2 

I , . m^i' 

-(z) H 

3^ ' 2 



(3.120) 
(3.121) 



and hence 
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1 



aiD{l)m2 + a2£,(l)?7ii = [{z) + ma] 



in contradiction with the HQET result eq. ( 3.32| ) 



A 



(3.122) 



(3.123) 



This fact suggests a necessity of subtraction in the quantity aim2 + a^mi. Let us write the spectral density with 
subtraction in the form 



0,17712 + AjTOi = aiDm2 + a2o7ni + 



K (Ml-Si , Af|-S2\ Pa 



with 



Pa ^ 2 ma 



+ 1 V 2^51 

zi + Z2 + 2m3 ' 
(I) + 1 



Then corresponding representation for the form factor reads 



ai(cj)m2 + a^(cj)mi = 4 /„Ln, / dz24'Q{z2) Z2{z2 + 2TO3) / ^(/)o(^i) 
aiD"^2 + a2Z5"^i ^ + ^2 + 2m3) (^77 



(3.124) 



(3.125) 



(3.126) 



zi+Z2+2m-i 
u> + l 



kA 



According to the HQET relation ( |3.32 ) this quantity should be equal to 

1 r 2A . ^ a'^-i 



^/mi77l2 



UJ+1 



(3.127) 



The term proportional ^(lo) yields k — A. One can check that this value also makes other parts of both expressions 
equal. Hence we arrive at the subtracted spectral density 



aiTO2 + a2TOi — aiDm2 + a2]j7ni 



1 (Ml~si M|-S2 



Pa- 



U +1 \ ^Si y/S2 

The resulting spectral densities of the form factors ai and 02 with the built-in subtraction terms take the form 

1 fMl-si , Ml-S2 \ Pa 



(3.128) 



oi = aiD 



0-2 = a2D 



{lu + l)m2 



1 /Mf-Si , Mi-S2\ Pa 



(w + l)mi \ ^sT 



(3.129) 
(3.130) 



The NLO analysis of gg. 

The unsubtracted spectral density goD has the form 

4 



90D 



At zero recoil one finds 



m2ai(m3 + m2ai + miQ;2) 



jz) + ms (f) 

(mim2)3/2 2 6 



/3 



u + 1 



On the other hand, taking into account our earlier finding X2 — 0, the HQET result (3.28) reads 



5od(1) 



1 



(mim2)'^/2 [2 6 



(3.131) 



(3.132) 



(3.133) 
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Our experience obtained in considering Oim2 + a'2mi hints that the subtraction procedure adds a term proportional 
at zero recoil to A — (z) — ms, and hence we expect subtraction to work properly also in the case of go- 
As a matter of fact, the spectral density 



90 = .9013 + 



1 



{uj + l)?7lim2 \ 



satisfies the HQ expansion ( p.2S| ) provided the function pg„ behaves in the LO as 

Zl + Z2 + 2to3 ' 
LU+1 



(3.134) 



(3.135) 



The scalar form factor 

Let us consider the form factor of the transition between two pseudoscalar mesons induced by the scalar current 



{P{M2,P2)\qiq2\P{Ml,Pl)) ^ fs{q^) 
For the analysis of the heavy-quark transition we introduce the velocity-dependent form factor as follows 

{P{M2,P2)\qiq2\PiMi,pi)) = ^/M^{1 + u;)hf^{u;) 

such that 



fsiq^) = ^AhM2il+io)hf^{L^). 
From HQET we find the following expansion 



uj + l \2 



An important consequence of this expansion is the relation 

hfSl) = l + 0{l/ml), 



Pi 



(3.136) 
(3.137) 
(3.138) 

(3.139) 

(3.140) 



which shows that the deviation of hf^ (1) from unity emerges only at the l/rriQ order, as well as for the form factors 
hf, hf^, and hg^ ( 3.14 ) . 

Notice that in the leading and subleading orders in 1/toq the following quantities have the same expansion 



fs ~ 2{pip2 - mim2)s + r7i2/i + "^l/2 
~ 2(pip2 - MiM2)s + Mi/i + Af2/2. 



The form factor can be obtained according to the relation 

1 



Mq') 



q>^{P{M2,P2)\qilMP{Mi,Pi)), 



m2 — mi 

which follows from the equations of motion for quark fields. We then find 

1 



fs{q') 



7712 — rai 



[(Mf-M|)/+(g^) + 'zV-(g^)] 



and 



7772 — "nil 



[hf^ (Ml - M2)[u + 1) - hf_ (Afi + M2){lo - 1)] . 



(3.141) 



(3.142) 



(3.143) 



(3.144) 



In this case the structure of the expansion (3.13E) is satisfied automatically. 

In the dispersion approach the double discontinuity of the form factor fs is found from the Feynman graph 
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/f = 2to3 [{mi - - q^) + 2rn2 (s2 - ("t-i - ^t^s)^) + 2rni (si - (7712 - mzf) - 4mi77i2TO3- (3.145) 

The double discontinuities of the form factors /+, /_ and fa obey the classical equation of motion 

1 



fs {si,S2,q 



7772 — 7771 



(si - S2)/+(si,S2,g )+q f-{si,S2,q ) 



(3.146) 



As becomes clear from comparing Eqs.( 3T4^ ) and ( 3.146| ), in order to satisfy the stucture of the expansion ( 3.139 ) 
the form factor requires subtraction. The spectral density which includes the proper subtraction term reads 



/. ^ /f + [{Mi - si) - (M| - S2)] 



7772 ~ "T-l 



(3.147) 



Equivalently, the behavior of the form factor fs can be brought in accordance with HQET in the NLO by defining 
the subtraction as follows 



provided pf^ 



fs = /f + [{Mi - si) + {Ml ~ S2)] Pf^ 
Pg, ~ Pa, Eqs. (|3.125|) and (^.1351). 



(3.148) 



Concluding this section let us summarize our main results: we have calculated the universal form factors and 
demonstrated the dispersion representations with relevant subtractions in the case of heavy-to-heavy transitions to 
reproduce the structure of the heavy-quark expansion in QCD in the leading and next-to-leading orders. However, 
we have not been able to fix uniquely these subtractions. As we shall see in the next section, the heavy-to-light 
transitions provide further restrictions on the form of the subtraction terms. 



D. Heavy— to— light meson transitions 



In this section we discuss meson decays induced by the heavy-to-light quark transitions in which case Mi = 
m2 + 0(1) is large, while M2 — mi is kept finite. As found by Isgur and Wise 1^, in the region uj — \ = 0{1) the 
form factors of the tensor current can be expressed through the form factors of the vector and axial-vector currents 
in the leading l/m2 order as follows 

^(9')-^/i(g'), (3.149) 
g2{q^)^~2Mig{q^), (3.150) 
5o(9') = 7V[25(«') + «2(<z')], (3.151) 



M- 



51(9') - + 2P1P2 • g{q^)]. (3.152) 

We address the two issues: (i) perform the leading order l/m2 expansion of the form factors and show the fulfillment 



of the relations (3.149-3.152) and (ii) discuss the scaling behavior of the form factors of the transition of different 



heavy mesons into a fixed final light state. 

The LO l/mg expansion of the form factors in the quark model 

In the case of heavy-light transitions one observes the appearance of the two scales: the light scale M2 ~ 7771 ~ 
7773 ~ A, and the heavy scale M2 — mi, and we may expand the form factors in inverse powers of the small parameter 
A/7772. The kincmatical coefficients in the leading A/7772 order in the region Cu — 1 = 0{1) simplify to 

A(si, S2, q^) ~ 47772 [(^1 + "^3 + ojmiY — S2] , (3.153) 
/3 = 0(1), (3.154) 

1 {Zl + 7773 + (I)77ll)(s2 — 777^ + 77l3)/2 — S2(^l + '^s) 



ai 



7772 (Zi + 7773 + a;777i)2 - S2 
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1 

7712 



a2 



ai2 



an 



{zi +1713+ ujmi){zi + ms) — (s2 



mi)/2 



Oil): 

a\a2 ■ 



(zi +1713 + OJmiY — S2 

[3{zi + 1713 + WTOl) 
2to2 [(^^l + TO3 + LUTTliy — S2] 



Ol-L 

TO2 



/3S2 



TO2 [{Zl + 7173 + LOTUiy 



S2 



o 



For the double spectral densities fi these expressions yield the following LO relations 

1 



s 

goD 
91 



2m2 
2^2.9, 
1 



/2, 



m2 
1 

TO2 ' 



-(2.9 + a2D), 
[-Id + 2pip2 ■ g] 



(3.155) 
(3.156) 
(3.157) 

(3.158) 

(3.159) 
(3.160) 
(3.161) 

(3.162) 



with 2pip2 = si + S2 — ■ Notice that these relations hold also for Co — 0(to2)- 

First two equations directly give the LO equality of the corresponding form factors. The relation (3.161) between 
the unsubtracted spectral densities is more interesting and requires the LO identity of the subtraction terms 



Pa - Pgo ■ 



(3.163) 



The choice pa — Pgg = 5 is acceptable al thoug h there are no firm backgrounds to justify this very choice. 

The most informative is the relation (3.162): this relation not only suggests a necessity of subtraction in the form 
factor / but also determines the subtraction ter m. Wh en considering the heavy-to-heavy transitions we have found 
that the subtracted spectral density of the form (3.106) 



f = fD + {M( - s, + - S2)pf 
matches the HQET expansion if ~ g in LO. The eq. ( p.l62| ) prescribes pf — g. 

ScaliTig of the form factors with tuq 



(3.164) 



Let us point out that in the case of the heavy-to-light transition there is no parametrical suppression of the 
anomalous contribution to the form factors compared with the normal one as it was in heavy-to-heavy transitions 
and both parts contribute on equal footings. We shall present all the results for the normal part but the same holds 
also for the anomalous part. 

In the region uj — 0{1) one finds for the normal contribution 

J ^--^ _ 1 f ds2(p{s2)X^^'^{s2, mj, 711,1) f dr]y/u;'^ -~ 1 (I)q{zi)I{zi, Z2,mi,m2,m3,u;) ^^^^ 

Vm^ J 167rmi J ^(mi((I) + 1) + zi + 22 + 2rn3)(mi((I' ~ 1) + zi - Z2) 



where 



S2 = {mi + 7713 + Z2) , 

/ 27773 +2;2\ , N A^/^(S2, 777?, 777o) , 

Zl = Z2W 1 + + m3{u> - 1) + T]V^^ - 1 ^ ' ' ^' + 0(l/r772) 
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Using the relations ( 3.153| ) one finds that the spectral densities scale at large m2 as 

= m^'- Pi{uj,mi,m3, zi, Z2). (3.166) 

Namely, 



/l 


= o(i), 


h 


-0(to2), 


9 




ai 


= 0{l/m2), 


0.2 


= o(i), 


f 


= 0(m2), 


91 




92 


= 0(^2), 


90 


-0(1/to2), 


s 


= 0(1). 



Hence, the form factors have the scaling behaviour of the form 

fi = m2''^^^ri{uj;mi,m3]ip2,4>o). (3.167) 
The variables u and a) are connected with each other as follows 

M9 A 

ui = —u; + 0(l/m2), (3.168) 

nil nil 



and hence the ratio 



/,(a-)/m'^'+^/^ (3.169) 



is universal for the transition of any heavy meson into a fixed light state. This behaviour reproduces the results of |g] 
up to the logarithmic corrections which arise from the anomalous scaling of the quark currents in QCD. 

To summarize, matching the LO l/m2 relations between the form factors in the quark model to the corresponding 
relations of ||^ allowed us to determine the subtraction term in the form factor / and equated to each other subtraction 
terms in the form factors 01^2 and go. 



E. Numerical estimates for the universal form factors 



In this section we apply the derived results to the model-dependent estimates of the universal form factors. 

1. First, let us study the dependence of the IW function on the parameters of the quark model. We adopt the 
exponential parametrization of the radial wave functions in the form w(/c^) ~ exp(— fc^/2/3^). For obtaining the LO 
wave function (/)o we must take into account the relationship between the parameters z and A;^ which is found from 
the following equation 



This yields the relation 



Vs — \/ k"^ + niQ + y k"^ + m\ — mg + 



= z{z + 2m3) + 0(1/toq). 



Using the eq. ( 2.93 ) we come to the following form of the LO wave function 

, f . j z + m3 ( z(z + 2m3)\ 

where /3o is the LO harmonic oscillator parameter 



(3.170) 



(3.171) 



(3.172) 
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/3(TOQ)=/3o + 0(l/mQ). 



(3.173) 



For calculations we need to specify (3o. The exact value of Po is not known but extrapolating the known values of 
/3d, /3b, and Pd' we find Po — 0.4 both in the WSB [|o| and the 1SGW2 Q models. In the model for pseudoscalar 
mesons discussed in the previous Chapter we have found the following approximate relation 



PiniQ) ~ 2.5 



(3.174) 



which gives Po — 2.5 . Tab le ^ presents the relevant numerical parameters. The results of calculating the IW 
f unct ion through the eq. (3.90) are shown in Fig. |l^. Table |^ demonstrates the values ^'(1) calculated with the eq. 
( 3.91 ) and the parameters of the quadratic fit of the form 



^io) = I - pUu^ - I) + S^iu; - 



(3.175) 



obtained by interpolating the results of calculations in the range 1 < w < 1.5. One can observe the value of — C'(l) to 
be considerably larger than the parameter obtained by the interpolation procedure. 



Table 4. Parameters of the quadratic fit to the IW function and the NLO form factor ^3 in various quark-model versions. 



Ref. 


ma 


Po 




p' 


S 


C3(l) 


Pi 




Set 1 
Set 2 1 
Set 3 i 


IC 

~ 


1 
1 
1 


0.35 
0.33 
0.25 


0.4 
0.4 
0.63 


1.47 
1.43 
1.37 


1.34 
1.31 
1.12 


0.78 
0.76 
0.62 


0.077 
0.08 
0.17 


1.49 
1.4 
1.2 


0.88 
0.84 
0.65 




Q 5 I \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 

1 1.1 1.2 1.3 1.4 1.5 

CO 



Fig. 16. The Isgur-Wise function calculated with various quark-model parameters: dashed line - set 1, solid line - set 2, 
dotted line - set 3. 

The IW function obtained from the dispersion quark model agrees with the results of the light-cone quark model 
applied to the elastic P —i- P transition for (infinitely) heavy mesons C(l-5) = 0.62 and turns out to be a bit 
smaller than the SR resuh ^(1.5) = 0.66 ||. 

The NLO form factor ^3(0;) is foi md to be very sensitive to the light-quark mass. The parameters of the quadratic 
fit to ^3(0;) calculated through eq. ( |3.100 ) in the form 

h,{Lj) = h,il) [1 - _ 1) + s.iuj - 1)2] (3.176) 

for various sets of the quark-model parameters are shown in Table ^. One can observe an approximate relation 
£^3{uj)/^{uj) ~ 0.08 GeV for the light quark mass ~ 0.35 GeV and ^3(a;)/^(w) ~ 0.17 GeV for the light quark mass 
7713 ~ 0.25 GeV in the whole region I < uj < 1.5. This ratio is to be compared with the SR result with the 0{as) 
corrections omitted j^] 

^3(tj)/^(w) =1/3-0.16 061^, A~0.5GeK (3.177) 
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2. We are in a position to estimate the higher order corrections as we can calculate the form factors at finite masses 
and the leading-order contribution separately. For the ISGW2 parameter set {rric = 1.82, mi, = 5.2, = 0.45, 
Pb = 0.43, and = 0.38), the results on hf^ and hf a re sho wn in Fig. and Table || presents the results of 
interpolating in the range 1 < a; < 1.5 with a quadratic fit (i.l7t). 



Table 5. The form factors calculated with the parameters of 1SGW2 model. 





hu 


hf 




p' 
s 


0.93 
0.87 

0.37 


0.96 
1.06 
0.53 


1.0 
1.25 
0.7 




Fig. 17. The velocity-dependent form factors calculated for the ISGW2 parameters: solid line - the IW function, dashed 
line - hf{ijj), dotted line - hf^{uj). 

The form factor h f can be seen to agree with the results of a combined fit to the ALEPH, ARGUS, CLEO, DELPHI, 
and OPAL data ^ 

= 1.07 ±0.20. (3.178) 

One observes sizeable difference between the form factors hf and hf_^ and the Isgur-Wise function both in the value 
at zero recoil and the slope parameter. This indicates sizeable higher order corrections. At zero recoil the l/mg 
corrections to hf vanish and writing Q 

/i/(l)==l + ,5i/,„2 (3.179) 

we find (5i/m2 — —0.04 that agrees well with the value —0.055 ± 0.025 0. At the nonzero recoil the higher-order 
corrections yield a more flat w-dependence of the form factor h f compared with the IW function such that the slope 
of the IW function in our model turns out to be considerably larger compared with the slope oi hf-. An approximate 
relation p]y^ ~ + 0.2 is found. This is different from the sum-rules result p]^ ~ — 0.2 Q. Let us also point 
out that h'{l) obtained by the direct calculation (the theoretical value) turns out to be considerably larger than the 
parameter of the quadratic fit obtained by the interpolation over the range 1 < uj < 1.5 which could be measured 
experimentally. The same happens for the Isgur-Wise function. 
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F. Discussion 



A detailed analysis of our dispersion approach to meson transition form factors was performed in this chapter. 
We calculated the double spectral densities of the form factors at spacelike momentum transfers from the Feynman 
graphs, and determined the necessary subtractions by considering the two following cases: 

(i) by performing the 1 /toq expansion of the form factors in the case of heavy-to- heavy transitions and matching 
them to HQET, and 

(ii) by considering the relations between the form factors of the vector, axial-vector and tensor currents in the case 
of the heavy-to-light transition and matching them to the general relations in QCD obtained in ref. 

This procedure led to the spectral representations with appropriate subtractions at spacelike momentum transfers. 
The form factors at timelike momentum transfers were then obtained by performing the analytical continuation in 
. The analytical continuation yielded the appearance of the anomalous contribution to the form factors at > 0. 

As a result, we have come to an explicit model of the form factors describing the weak meson transitions which 

develops the correct structure of the heavy-quark expansion in agreement with HQET. 

Let us emphasize that the obtained representations allow to calculate form factors both at spacelike and timelike 

momentum transfers in terms of the wave functions of the participating mesons. 

The following results were obtained in this chapter: 

1. Applying the method of ref. to meson amplitudes induced by the heavy-quark tensor current we have found 
the expansions of the relevant form factors in LO and NLO in HQET. The 1/toq corrections in the form factor hg^ 
vanish at zero recoil just as in /i^^ and hf. 

2. The l/mg expansion of the soft wave function was constructed. The normalization condition of the wave function 
yields an infinite chain of normalization conditions for wave function components emerging in various 1/toq orders. 
The LO normalization condition provides the proper normalization of the Isgur-Wise function at zero recoil. 

3. Form factors of heavy-to-heavy meson transitions were considered. 

First, the leading l/mg-order analysis of the spectral representations for the form factors was performed. It 
was shown that to this accuracy all the transition form factors are represented through the Isgur-Wise function in 
accordance with HQET. The Isgur-Wise function was calculated in terms of the leading-order component of the heavy 
meson wave function. Subtraction terms in the spectral representations for the form factors do not contribute to this 
accuracy. 

The anomalous contribution which emerges when the analytical continuation to the timelike region is performed 
comes into the game only in the vicinity of the zero recoil point. It is otherwise suppressed by at least the 

second inverse power of the heavy-quark mass. 

We then analysed spectral representations for the form factors with next-to-leading order accuracy. 

We found that the P — > P form factors do not require any subtractions to this accuracy and calculated the universal 
form factor ^3 in terms of the LO wave function. 

In the case of the P ^ V transition the representations without subtractions perfectly match HQET for the form 
factors (7, (71, ^2, and mia2 — 771201 and leads to the relation pi = p2 (or xs = 0) and X2 = 0. At the same time, 
spectral representations without subtractions for the form factors /, go, and 771102 -I- TO201 do not agree with HQET 
to this order. 

Matching of the dispersion representations for these form factors to HQET shows the necessity of subtractions and 
restricts the form of the subtraction terms. 

For performing the heavy-quark expansion only a strong localization of the meson momentum distribution with a 
width of order of the confinement scale was necessary. No other constraints on the soft wave functions or numerical 
parameters of the model have emerged. 

4. We analysed heavy-to-light meson transitions and found that the Isgur-Wise relations between the form factors of 
the tensor and vector and axial-vector currents further constrain subtraction terms in the spectral representations 
for the form factors. 

5. We observed a discrepancy between the predictions of the various versions of the quark model on the NLO universal 
form factors : Namely, the analysis of the Wirbel-Stech-Bauer model [|o| performed in |^] resulted in pi 7^ p2 7^ and 
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P3 7^ 0, but p4 = 0; the quasipotential quark model predicts all the NLO form factors to be nonzero; a consistent 
relativistic treatment of the qq intermediate states in Feynman graphs performed in this paper gives pi = P2, P3=0, 
and /34 7^ 0. 

The diagreement between our relativistic approach and the WSB model [|l0| can be traced back to the not fully 
consistent treatment of the quark spins in the latter. On the other hand, the origin of the discrepancy between 
the relativistic quasipotential approach of and our dispersion approach is not fully understood and should be 
considered in more detail. 

6. Numerical results from our dispersion approach using several parameter sets (i.e. constituent quark masses and 
wave functions) demonstrate only a moderate dependence of the Isgur-Wise functions on the value of the light-quark 
mass and the shape of the soft wave function. 

Namely, the light-quark mass ma = 0.25-;-0.35 and the wave-function width /3o — 0.4-^0.65 give ^^(1.5) = 0.55-^0.6 
that is a bit smaller than the SR result ^(1.5) = 0.66 

The form factor h f is found to behave in agreement with predictions of other models and experimental results. 

We observed a sizeable difference between the form factor h f and the Isgur-Wise function both in the value at zero 
recoil and the slope parameter. This difference might be important for the analysis of the experimental data for the 
B {D,D*)lv decays. 

The size of the higher-order l/mg-corrections to hf{l), (5i/,„2 = —0.04, obtained using the ISGW2 parameters 
agrees favorably with the sum rule estimates. 

The slope of the Isgur-Wise function in our model turns out to be considerably larger compared with the slope 
of hf. An approximate relation pj^ ~ pf^^ + 0.2 (see Fig. |l^ is found. This is opposite to the sum-rules result 
Pj^y ~ pf^^ — 0.2 [0. We would like to notice that the slope parameter is very sensitive to the interpolation procedure: 
namely, the value h'{l) turns out to be considerably larger than the result of the quadratic interpolation over the 
range 1 < w < 1.5. 

Our dispersion approach to the transition form factors can be further improved by performing the heavy-quark 
expansion in higher orders and matching to HQET order by order. It should be taken into account that the quark 
model effectively describes the whole I/ttlq series, but the short-distance corrections are not contained in the model 
and should be considered separately. 
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IV. CALCULATION OF THE WEAK FORM FACTORS 



In this Chapter the dispersion approach developed in previous sections is apphed to form factors of weak transitions 

|. We present form factors for weak decays of B/^g) and 
mesons to hght pseudoscalar and vector mesons obtained by the recent detailed analysis | |6l[ |. 

The use of this dispersion approach based on the constituent quark picture allows us to reveal the intimate connection 
between different decay modes and to perform the calculations in the full physical g^-region. In fact quark models 
provide for the only approach which leads to relations between the decays of various mesons through the meson wave 
functions and gives the form factors in the full range. However, quark models arc not closely related to the QCD 
Lagrangian (or at least this relationship is not well understood yet) and therefore have input parameters which are 
not directly measurable and may not be of fundamental significance. 

We reduce the usual disadvantages of quark models related to ill-defined effective quark masses and not precisely 
known meson wave functions by fitting the quark model parameters to lattice QCD results for the B p transition 
form factors at large momentum transfers and to the measured total D {K,K*)lv decay rates. This allows us to 
predict numerous form factors for various decay channels and for all kinematically accessible values. 

As demonstrated in the previous Section, the calculated form factors have the correct structure of the long-distance 
corrections in accordance with QCD in the leading and next-to-leading l/mq orders, if the radial wave functions wikF') 
are localized in a region of the order of the confinement scale, fc^ < A^. We assume a simple gaussian parametrization 
of the radial wave function 

w;(fc2) cx exp(-fcV2/3^), (4.1) 
which satisfies the localization requirement for [3 ~ ^qcd- 



of heavy mesons. This issue was studied in 
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A. Parameters of the model 



We consider the slope parameter (3 of the meson wave function (4.1) and the constituent quark masses as fit 
parameters. The relevant values for the -B, p, and tt mesons are determined from the requirement to reproduce 
the observed value of the pio n de cay constant f-n = 130 MeV and from the fit to the lattice results on the form factors 



T2{q^) and Ai{q ^) (see Eq. ( [l.3| ) below) a t q^ — 19.6 and 17.6 GeV^ fs^. The spectral representations for the decay 
constants (3.63) and for the form factors ( ^.37 ) are used for calculations. 
The values obtained for m;,, m„, and Pb, /3p, P-k are displayed in Table ^. 



Table 6. Constituent quark masses and slope parameters of the exponential wave function (in GeV). 





nis 


nic 


mt 


P^ 


Pk 


Pd 


Pb 






Pb^ 


Pp 


Pk* 


Pd- 


P^ 


0.23 


0.35 


1.45 


4.85 


0.36 


0.42 


0.46 


0.54 


0.45 


0.48 


0.56 


0.31 


0.42 


0.46 


0.45 



A few comments are in order: 



As it was already noticed, the spectral representation ( ^.37] ) takes into account the long-distance contributions 
connected with the structure of the initial and final mesons. To describe additional long-distance effects in the 



channel, Eq. (3.37) should be multiplied by the constituent quark form factor fq^^q^^iq^) which contributes 
to the resonance structure in the q^ channel. However, in the region of calculation < (m2 — toi)^, the wave 
functions provide already for a rise of the form facotrs with q^ which is well compatible with a properly located 
meson pole. Thus an additional quark form factor is not needed there, but we shall use a proper extrapolation 
formula when considering the vicinity of the poles. 

The quark model double spectral representations take into account long-range QCD effects but not the short- 
range perturbative corrections. However, by fitting the wave functions and masses to reproduce the lattice 
points, these corrections are effectively taken care of: Corrections to the quark propagators correspond to the 
appearance of the effective quark masses. Corrections to the vertices at the relevant values of the recoil variable 
CO = {Mg + — q^) /2MbMt^ should be small as found in form factors of other meson transitions p9[ . 

The value obtained for the 6-quark mass mb = 4.85 GeV is close to the one-loop pole mass which in fact is the 
relevant mass for quark model calculations. 
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We now need to fix the parameters describing the strange and charmed mesons. The charm quark mass can be 
determined from the well-known 1 /mg expansion of the heavy meson mass in terms of the heavy quark mass and the 
hadronic parameters A, Ai and A2. Using the recent estimates of these parameters one finds 



3.4 GeV. 



(4.2) 



This provides rric — 1.45 GeV. For one expects ~ 350 — 370 MeV taking into account that niu — 230 MeV. 

A stringent way to constrain the parameters rris, (3k, Pk', and Pd is provi ded b y the measured integrated rates 
of the semileptonic decays D {K,K*)li/. In addition we apply the relation (3.63) which connects Pk with nis by 
using the known value of the if-meson leptonic decay constant /k — 160 M eV. The parameter values found this way 
are displayed in Table The corresponding form factors and decay rates are given in Tables | and 0. 

The polarization of the K* in the D K*lv decays turns out to be in good agreement with the experimental 
result (Table |ll|), and the calculated D meson decay constant fo = 200 MeV corresponds to the expectation for the 
magnitude of this quantity. 

The parameter (3d* cannot be found this way, but it should be close to (3d because of the heavy quark symmetry 
requirements. We therefore set (3d' = (3d- 

Also listed in Table |^ are the parameters which describe strange heavy mesons. They are discussed in subsection 

D. 



Table 7. 
Table p. 



Leptonic decay constants of pseudoscalar mesons in MeV calculated via ( 3.6S ) for the quark-model parameters of 



u 


/a- 


Id 


!b 








132 


160 


200 


180 


183 


220 


200 



Table 8. Meson masses in GeV from PDG @ 





Mk 


Mr, 




Md 




Mb 


Mb, 


Mp 


Mk* 


M^ 


Md* 


Mot 


0.14 


0.49 


0.547 


0.958 


1.87 


1.97 


5.27 


5.37 


0.77 


0.89 


1.02 


2.01 


2.11 



The knowledge of the wave functions and the quark masses allows the calculation of the form factors in Eq. (3.4). 
It is however more convenient to present our results in terms of the dimensionless form factors F+, Fg , Ft, V, Aq, 
Ai, A2, Ti, T2, which are the following linear combinations of the form factors given in Eq. (3.4): 



F+ 

Fo 

Fx 
V 

Ai 

A2 

Ao 

Ti 

T2 



q 



Pq 



(Ml 
(Ml 



1 



M2)S, 

M2)9, 



Ml + M2 • 
-{Mi+M2)a+, 



Pq 
Pq 

9- - ^50- 



9-, 



Pq ■ a_, 



(4.3) 



^In |2^ a different set of the parameters was used which also provided a good description of the available experimental data 
on semileptonic B and D decays. However, the corresponding form factors have a rather flat g^-dependence and do not match 
the lattice results at large . 
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Each of these form factors contains contributions of g^-channel resonances only of the sam e sp in. This is the 



advantage of the form factors of the set ( [4.3| ) compared with the form factors defined by the Eq. (3^), some of which 
contain contributions of resonances with different spins. 

1. The form factors F+, Ft, V , Ti, and Aq 

The form factors F-^-, Ft, V, Ti contain a pole at — My = and Aq contains a pole at = Mb = M^_. 
The residues of the form factors at these poles are given in terms of the product of the weak and strong coupling 
constants. These coupling constants therefore determine the behavior of the form factors at q^ near the resonance 
poles (beyond the decay region). We consider as an example the case of the B ^ iT,p transition. 

Weak decay constants 



Weak decay constants of mesons are defined by the following relations 

(0|g7p75&|S(g)) = ^f^p''\^. 
mi,b\B*{q))=e<^')MB^f\^'^ 

ma^,h\B*{q)) - ^(e;(^*)g, - €<^''> q^)f^/'\ (4.4) 
where e/i is the B* polarization vector. In the heavy quark limit one has 

AB) _ AB') _ AB') 
J P ~ Jv ~ JT 

Strong coupling constants 

Strong coupling constants are connected with the three-meson amplitudes as follows 

{n{p^)B*{q)\B{p,)) ^ -]-gB.B.P,ef'> 



{p{p2)B*{q)\B{p,)) = Lp^,e(P^e\pP^q, 



{p{p2)B{q)\B{p,)) - ^gBBpP^elrK (4.5) 

where q = pi — p2, P = Pi + P2 and is the polarization vector of the vector meson. 
In the heavy quark limit 

gs'Bp — gsBp- 

The form factors F+, Ft, V, Ti contain pole at q^ — M^, due to the contribution of the intermediate B* (1^ 
state) in the channel. The residue of this pole is given in terms of the product of the weak and strong coupling 
constants such that in the region q^ ~ Mg, the form factors read 



Ft 



V = 



AB') 
9B-B7rJv 


1 




2Af* 


1 ~ q^/Ml, 


+ ..., 


AB') 
9B'B-nJT 


Mb + M^ 


1 


2M* 


M* 1 


-<flMl, 


gB'BpJv 


Mb + Mp 


1 


2M* 


Mb- 1 




AB') 
QB'BpJt 


1 




2Af* 


1 - q^/Ml, 


+ ... 



(4.6) 



Here . . . stand for the terms non-singular at q = M^, 
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Similarly, contains the contribution of the B (0 state). In the region of ~ M'^ it can be represented as 
follows 

. ^ gsBpff^ Mb 1 

° 2Mb 2Mp 1 - 92 /Ml ^ - ' ^ 

Let us notice that the residues of the form factors are not all independent and are connected with each other as 
follows: 

Res{FT)Res{V) _ Mb + Mp Mb + . 
Res{F+)Res{Ti) Mb- Mb* ' 

This relation can be used as a cross-check of the consistency of the extrapolation for the form factors. 
The coupling constants are related to the residues of the form factors according to the relations 



AB') 
gB'B-nJv 



2Mb* 



2Mp 



ResiF+), 
2Res{Ao), 



f^""'^ _ ResiFr) Mb* 



f(B*) Res{F+) Mb + M^' 

93* Bp _ Res{V) Mb* 
9B*B^ ~ Res{F+) Mb + Mp ' 



(4.9) 



2. The form factors Fq, Ai, A2, T2 and T3 



The remaining form factors, Fq, Ai, A2, T2 and T3, do not contain contributions of the lowest lying negative parity 
states (for instance, Fq contains a contribution of the O"*" state, and Ai contains that of the l"*" which have considerably 
higher masses). As a result they have a rather flat behaviour in the decay region, whereas the form factors F+, 
Ft, V, Ti, Aq are rising more steeply. 



From the spectral representations ( 3.37 ) together with the parameter values of Table I the form factors are obtained 
numerically. For the applications it is convenient, however, to represent our results by simple fit formulas which 
interpolate these numerical values within a 1% accuracy for all values in the region < q^ < (7712 — mi)^. Also, 
they should be appropriate for a simple extrapolation to the resonance region. 

Let us start with the form factors F^, Ft, V, Ti, Aq. If we interpolate the results of the calculation with the simple 
three-parameter fit formula 

f(n^) = (4 10) 

■'^^ ' {l-q^/M^){l-qy{aMy)' ^' ' 

the least-x^ interpolation procedure leads in all cases to a value of the parameter M which is within 3% equal to the 
lowest resonance mass. We consider this fact to be an important indication for the proper choice of the quark-model 
parameters and for the reliability of our calculations. We therefore prefer to fix the pole mass M to its physical value. 



The fit functions (4.10) represent the results now with an accuracy of less than 2%. To achieve the accuracy of less 



than 1% in all cases we take the form 

f(n^) = (4 11) 

■'^^ ' il-qyM^)[l-aiqyM^ + a2qVM^y ^' ' 

where M = Mp for the form factor Aq and M = My for the form factors F^, Ft, V, Ti. In the Tables below we 
quote numerical values of (T2 only if an accuracy of better than 1% cannot be achieved with (T2 = 0, and take (72=0 
if this accuracy can already be achieved with the two parameters /(O) and cti. 

For the heavy-to-light meson transitions the masses of the lowest resonances are not very much different from the 



highest q^ values in the decay. Eq. (4.11) then allows an estimate of the residues of these poles. These residues can 



be expressed in terms of products of weak and strong coupling constants. The errors for these constants induced 
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by changing ai and (T2 in our fitting procedure (keeping to the 1% requirement) do not exceed 10% . Moreover, the 



residues of the form factors at the meson pole are not independent and satisfy certain constraint (i^), which provides 
a consistency check of the extrapolations. The mismatch in ( |4.§| ) is always below 10%, and in most of the cases much 
lower. 

For the form factors Fq, Ai, A2, T2 and T3 the contributing resonances (0+, 1+, etc) lie farther away from the 
physical decay region and the effect of any particular resonance is smeared out. For these form factors the interpolation 

f{q^) = /(0)/[l - TiqVM^ + <T2qVM^]. (4.12) 

If setting (72 = allows us to describe the calculation results with better than 1% accuracy for all q^, a simple 
monopole two-parameter formula is used. 

The values of /(O), cti, and CT2 are given for each decay mode in the relevant sections. 



* One should note that the parameters ai and (T2 i n the fit formula (i.l2) for the form factors Fq, Ai, A2, T2, and T3 are 
introduced in a different way than in the fit formula (4.11) for the form factors F+, Ft, V, Ti, and Aq. 
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B. Charmed meson decays 



1. D-^ K,K* 



The D ^ K, K* decays are induced by the charged current c ^ s quark transition. As described in the previous 
section, the measured total rates of these decays are used for a precise fit of the parameters of our model. With the 
parameters of Table ^ we obtain the form factors listed in Table ^. Table |l^ compares the form factors at = with 
the results of other approaches and Table |l^ presents the decay rates. 

Table 9. The D ^ K,K* transition form factors. My = M^,- = 2.11 GeV, Mp = Md, = 1.97 GeV. For the form factors 
F+,Ft,V,Ao,Ti the fit formula Eq. (1.11) is used, for the other form factors - Eq. (4.12) 





K 


D ^ K* 1 




F+ 


Fo 


Ft 


V 






M 


Ti 


Ta 


Ts 


f(0) 


0.78 


0.78 




0.75 


1.03 


0.76 


0.66 


0.49 




0.78 


0.78 


0.45 




0.24 


0.38 




0.27 


0.27 


0.17 


0.30 


0.67 




0.25 


0.02 


1.23 


0-2 




0.46 










0.20 


0.16 






1.80 


0.34 



Table 10. Comparison of the results of different approaches for the semileptonic D ~* K, K* form factors at (f 



Ref. 


F+{Q) 


Ft(0) 


\/(0) 


Ai(0) 


A2(0) 


This work 


0.78 


0.75 




1.03 




0.66 


0.49 


WSB [he 




0.76 






1.3 




0.88 


1.2 


Jaus'96_| 




0.78 






1.04 




0.66 


0.43 


SR If 




0.60(15) 






1.10(25) 


0.50(15) 


0.60(15) 


Lat (average 


il 


0.73(7) 






1.2(2) 


0.70(7) 


0.6(1) 


Lat M 




0.71(3) 


0.66(5) 












Exp g 




0.76(3) 




1.07(9) 


0.58(3) 


0.41(5) 


Table 11 


The D {K,K*)lu decay rates in 10^° 


s ^ obtained within different approaches, \Vcs \ = 0.975. 


Ref. 


r(D ^ K) 


r(D 


^ K*) 


r(A"*)/r(A") 


Tl/Tt 


This work 




9.7 


6.0 




0.63 




1.28 


Jaus'96_p4| 




9.6 


5.5 




0.57 




1.33 


SR &~ 




6.5(1.5) 


3.8(1.5) 




0.50(15) 




0.86(6) 


Exp Q 


9.3(4) 


5.7(7) 


0.61(7) 


1.23(13) |8| 


Extrapolating 


the form factors to 


g2 = Ml, (or = 


= Ml 


br Ao) gives 


the following 


estimates of the coupling 



constants 



ADs 
gD.DK'Jp 



2Mk' 
ADD 

JT 

ADt) 
Jv 

gO'DK 



= 1.05 ±0.05, 
= 1.7±0.1, 
= 0.95 ±0.05, 
= 1.1 ±0.1. 
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2. D -^TT,p 



These decays are induced by the c d charged current. Since all the necessary parameters have already been fixed, 
this mode allows for parameter-free predictions. Table |l^ presents the results of our calculations. In Tables 13 and 
[l^ we compare our results with different approaches and with experimantal data. 



Table 12. The calculated D ^ n,p transiti on fo rm factors. My = Md* = 2.01 GeV, Mp = M p = 1.87 GeV. For the form 
factors F+, Ft ,V, Ao,Ti the fit formula Eq. (1.11) is used, for the other form factors - Eq. (1.12). 





D 


D^p 1 




F+ 


Fo 


Ft 


V 


Ao 


Ai 


A2 


Tl 


Ta 


n 


,/(0) 


0.69 


0.69 


0.60 


0.90 


0.66 


0.59 


0.49 




0.66 


0.66 


0.31 


(Tl 


0.30 


0.54 


0.34 


0.46 


0.36 


0.50 


0.89 




0.44 


0.38 


1.10 


(72 




0.32 
















0.50 


0.17 



Table 13. Comparison of the results of different approaches for the semileptonic Z) — > tt, p form factors at q'^ — 0. 



Ref. 


F+iO) 


-Ft(O) 


V{0) 


Ai(0) 


A2{0) 


This work 


0.69 


0.60 


0.90 


0.59 


0.49 


wsB 


0.69 




1.23 


0.78 


0.92 


Jaus'96 [p-4| 


0.67 




0.93 


0.58 


0.42 




0.50(15) 




1.0(2) 


0.5(2) 


0.4(2) 


Lat(ave_)_P5| 


0.65(10) 




1.1(2) 


0.65(7) 


0.55(10) 




0.64(5) 


0.60(7) 









Table 14. The D {TT,p)lf decay rates in 10^° , \Vcd\ = 0.22. 



Rcf. 


r{D tt) 


r{D ^ p) 


r(p)/r(7r) 


Fl/Tt 


This work 


0.95 


0.42 


0.45 


1.16 


WSB m 


ol 


0.68 


0.67 


1.0 


0.91 


Jaus'96j 


ij 


0.8 


0.33 


0.41 


1.22 


SR pg 




0.39(8) 


0.12(4) 




1.31(11) 


Mehkhov'97j27^] 


0.62 


0.26 


0.41 


1.27 


Exp fS^- 


S9| 


0.92(45) 


0.45(22) 


0.50(35) 





The form factors at = are close to the predictions of the relativistic quark model of Ref. |1J], but the 
dependence is different such that our model and |Q predict different decay rates. Although the experimental errors 
are very large and nearly all theoretical results agree with experiment, we notice perfect agreement of our decay rates 
with the central values. 

For the coupling constants we get the following relations 

AD') AD) AD') 

9D'D.Iv ^105^0 05 9DDpJp ^2.1 ±0.2, ^g g^O^^ 5£:£p ^ 1.3 ± 0.2. 

2Md' 2Mp jiD') g^,^^ 

Taking /^^*^ ~ 220 MeV, we find 

Qd'Dtt = 18 ± 3 

in good agreement with a calculation of go'D-K based on combining PCAC with the dispersion approach []60| . This 
value also comapres very well with the first recent measurement of go'DTv reported by CLEO | ]90| 

gjj^ij^ = 17.9 ±0.3 ±1.9. 

Notice a serious disagreement between our result and the sum rule estimate [ pl| 

= 12.5 ±1.0. 

The experimental result clearly speaks in favour of our prediction. 
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C. Beauty meson decays 



1. B ^ D,D* 



These decays arise from the heavy-quark b ^ c transition. Here one has rigorous predictions for the expansion of the 
form factors in terms of the heavy-quark mass Q . Namely, the main part of the form factors can be expressed through 
the universal form factor - the Isgur-Wise function. However, different models provide different g^-dependences of the 
Isgur-Wise function as well as different subleading 1/toq corrections. 

We recall that our spectral representations of the form factors explicitly respect the structure of the long-distance 
QCD corrections in the leading and the subleading orders of the heavy-quark expansion. Thus, we expect reliable 
predictions for the form factors. Our numerical results are summarized in Tables |l^, |l^, and p^ . 



Table 15. The B ^ D,D* tra nsitio n form factors. My = Mb* ^ Mp 
F+, Ft ,V, Ao,Ti the fit formula Eq. (i.ll) is used, for the other form factors 



Me 



i.4 GeV. For the form factors 



Eq. (4.12) 





B ^ D 


B ^ D* 1 




F+ 


Fo 


Ft 


V 




^1 




Ti 


Ta 




/(O) 
CTl 

0-2 


0.67 
0.57 


0.67 
0.78 


0.69 
0.56 


0.76 
0.57 


0.69 
0.58 


0.66 
0.78 


0.62 
1.40 
0.41 


0.68 
0.57 


0.68 
0.64 


0.33 
1.46 
0.50 



Table 16. Comparison of the results of different approaches for the semileptonic B ^ D, D* form factors at — 0. 



Ref. 


F+(0) 




i?i(o) = v(o)Mi(o) 


i?2(o) = ^2(o)Mi(o) 


This work 


0.67 


0.66 


1.15 


0.94 


Jaus'96 (li 


0.69 


0.69 


1.17 


0.93 


Neubcrt j|7[ 






1.3 


0.8 


Exp 1^ 






1.18±0.15±0.16 


0.71±0.22±0.07 



Table 17. The B {D, D*)lv decay rates in | ps~^ . 



Ref. 


r(B ^ D) 


T(B D") 


T(D'')/T{D) 


Vl/Vt 


This work 


8.57 


22.82 


2.66 


1.11 


Jaus'96 n 


9.6 


25.33 


2.64 




Mer97 If 


8.7 


21.0 


2.65 


1.28 


Exp 


(1.34 ±0.15)10"^ ps'^ Q 


(2.98 ±0.17)10"^ ps-' 


2.35(1.3) 


1.24(0.16) ^ 



For the coupling constants we find 



QB'BDJv 

ABa) 

QB^BD-Ip 

ab:) 

JT 

ABi) 

Jv 

gB'BD' 
9B'BD 



= 1.56±0.15, 
= 3.3 ±0.3, 
= 0.9 ±0.1, 
= 1.05 ±0.05. 
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2. B ^ K, K* 



These decays are induced by the h s Flavor Changing Neutral Current (FCNC). We recall that the B ^ Tr,p 
form factors at large have been used to fix the parameters of the model. Thus we expect that the predictions 
for the B ^ K, K* form factors which in fact differ from the former mode only by SU(3) violating effects should be 
particularly reliable. Table ^ presents the calculated form factors and Fig |l8| exhibits our predictions together with 
the available lattice results at large . The good agreement shows that the size and the sign of the SU(3) violating 
effects are correctly accounted for. 



Table 18. The calculated B 
form factors 



K, K* transition form factors. Mi 



V = Mb* = 5.42 GeV, Mp = Me 
Ft, V, ^0, 21 the fit formula Eq. ( 4.11 ) is used, for the other form factors - Eq. ( 1.12 ) 



5.37 GeV. For the 





B ^ K 


B ^ K* 1 




F+ 


Fa 


Ft 


V 




^1 




Ti 


Ta 


Ts 


/(O) 


0.36 


0.36 




0.35 


0.44 


0.45 


0.36 


0.32 




0.39 


0.39 


0.27 




0.43 


0.70 




0.43 


0.45 


0.46 


0.64 


1.23 




0.45 


0.72 


1.31 


0-2 




0.27 










0.36 


0.38 






0.62 


0.41 



Table 19. Comparison of the results of different approaches for the B ^ K, K* form factors at = 0. 



Rcf. 


F+(0) 


Ft(0) 


vm 


Am 


^2(0) 


Ao(0) 


Ti(0) 


T3(0) 


This work 


0.36 


0.35 


0.44 


0.36 


0.32 


0.45 


0.39 


0.27 


SR@ 


0.25 




0.47 


0.37 


0.40 


0.30 


0.38 




Lat+StechJM 






0.38 


0.28 




0.32 


0.32 




LCSR'98 ^ 


0.34 


0.374 


0.46 


0.34 


0.28 


0.47 


0.38 


0.26 


Lat H 


0.30(4) 


0.29(6) 
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0.8 — 
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Fig. 18. Form factors of the B K* transition vs the lattice results. 



For the coupling constants we obtain 



a^'^^^fy'^ ^ 0.65 ± 0.05, 9B.BK'.flf'^ ^ ^ gj. ^ 
2Mb. 2Mi^. 

It qS'BK' 

= 0.95 ± 0.05 ^ ^1.15 ±0.05. 
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3. B ^TT,p 



The B p transition has been used for determining the parameters of our quark model in the u, d and b sectors 
by fitting the quark- model form factors to available lattice results on T2 and Ai at large see Figure 19, The 

parametrizations for the B n, p form factors are given in Table 

The form factor Fq at large lies below the central lattice values but nevertheless agrees with lattice results within 
the given error bars. Notice however that in our model the form factor Fq is calculated as a difference of /+ and /_ 
and at large q^ turns out to be much more sensitive to the subtle details of the pion wave function, than /+ and /_ 
separately. A simple Gaussian wave function which works quite well for /+ and /„ , might not be sufficiently accurate 
for Fo. 



Table 20. The calculated B ^ Tv,p transiti on fo rm factors. My = Mb* = 5.32 GeV, Mp = M b = 5.27 GeV . For the form 



factors Fj^,Ft,V,Aq,Ti the fit formula Eq. (i.ll) is used, for the other form factors - Eq. (4.12) 





B 


B^p 




F+ 


Fo 


Fx 


V 


^0 


^1 


A2 


Ti 


T2 


Ts 


/(O) 


0.29 


0.29 


0.28 


0.31 


0.30 


0.26 


0.24 


0.27 


0.27 


0.19 


0"! 


0.48 


0.76 


0.48 


0.59 


0.54 


0.73 


1.40 


0.60 


0.74 


1.42 


C2 




0.28 








0.10 


0.50 




0.19 


0.51 
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Table 21. Comparison of the results of different approaches on weak B ^ tv, p form factors at — 0. 



Ref. 


F+(0) 


FriO) 


1/(0) 


Ai(0) 


A2{0) 


Ao{0) 


ri(o) 


7^3(0) 


This work 


0.29 


0.28 


0.31 


0.26 


0.24 


0.29 


0.27 


0.19 


Jaus'96 H 


0.27 




0.35 


0.26 


0.24 








LCSR'98^ 


0.305 


0.296 


0.34 


0.26 


0.22 


0.37 


0.29 


0.20 




0.28(4) 


0.28(7) 















Table Ell compares the results obtained from the quark model of Ref. ||5^ with results from the quark model 
of Jaus [0 and latest light-cone sum rule (LCSR) results |3^]. One observes very good agreement between the 
quark model of Jaus, LCSRs, and our approach. The only visible difference with the LCSR method occurs in 
the form factor Ao{0), which is caused by small differences of the two methods in ^i(O) and ^2(0) (recall that 
Ao(0) = ((Ml +M2)Ai(0) - (Ml - M2)A2(0))/2M2) This discrepancy exceeds the 15% error bar quoted for the 
LCSR results only marginally. If the LCSR results at small and lattice results at large are correct, our approach 
surely provides a realistic description of the form factors at all kinematically accessible values. 

Extrapolating the form factors to the poles, we obtain 



AB'} 

gB'B^fv 



2Mb' 



2Mp 

AB'} 

JT 

AB') 

Jv 



0.6 ±0.05, 
1.4 ±0.2, 

0.97 ±0.03, 
1.2±0.1. 



It is convenient to define the quantity g according to the relation 



-gBB'TT- 



Using /^* 



200 MeV gives the estimate 



2^MbMb 



gBB'TT = 32± 5, 
g = 0.4 ±0.06 



(4.13) 



The latter value is inj 
than g = 0.5 ± 0.02 |6C 



;ood agreement with the lattice result g = 0.42 ± 0.04 ± 0.08 
] based on combining PCAC with our dispersion approach. 



(4.14) 

and is only slightly smaller 



Table 22 compares the calculated semileptonic decay rates from different theoretical approaches. 



Table 22. The B —» {■Kp)lv Decay rates in units |K,6|^ps ^ ■ 



Ref. 


T(B Txlv) 


r{B p£iy) 


Tl/Tt 


This w( 
ISGW2 01 
Lat fc 
LCSR [ 
LCSR [ 


3rl 
VI 

37 
31 




1 
1 


9.6 

8.85 


15.8 ± 2.3 
14.2 

13.5± 4.0 
17.7 


0.88 ±0.08 
0.3 

0.52± 0.08 
1.12 
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Combining the experimental data reported by CLEO B{B° p~l^v) = {2.^7 ±Q.2^°_:1\q ± 0.41) • 10"^ (96) with 
our prediction for T{B plv) leads to the following estimate for the central value of Vub'- 

iKbl =3.25-10-3. (4.15) 

Using the B° tt'I+v measurement from Belle B{B° -> tt-I+v) = (1.28 ± 0.20 ± 0.26) • 10"'* ^ and our prediction 
for T{B — > nil/) leads to the central value 

iKbl =3.22-10-3. (4.16) 
The perfect agreement between these values speaks in favour of the reliability of our predictions for the form factors. 



Summing up our results on the decays of the nonstrange heavy mesons, we found no disagreement neither with the 
exisiting experimental data nor with the available results of the lattice QCD or sum rules in their specific regions of 
validity. The only exception is the form factor Fq at large in. B n and D —> n decays, where our results are 
lying slightly below the lattice points. However, this disagreement can be related to a strong sensitivity of Fq at large 

to the details of the pion wave function. Small changes in the pion wave function, which only marginally affect /+ 
and /_, can change the form factor Fq. But such subtle effects are beyond the scope of our present analysis. 

In the next section we apply our model to the decays of strange heavy mesons for which a few new parameters have 
to be introduced which are specific to the description of weak decays of strange heavy mesons to light mesons. 
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D. Decays of the strange mesons Ds and Bs 



Before dealing with these decays, we must first specify the slope parameters of the Bs and the Dg wave function. 

) and using fs^/fB = 1-1 and fo^l fo — 1-1 in agreement with the 
lattice estimates for these quantities ^5|] . The resulting values of the slope parameters are listed in Table ^. Since all 
other parameters have already been fixed the calculation of the form factors is straight forward. The only exceptions 
are the decays into the 77, 77', (/) final states. For these decays we need to know the (j) wave function, the mixing angle 
and the slope of the radial wave function of the ss component in 77 and rj' . Our procedure of fixing these parameters 
are discussed in the relevant subsection. 



We obtain these parameters by applying (3.63 



1. Ds 



K,K* 



These meson transition are driven by the charged-current c —f d quark transition. The results of the calculation 
are given in Table psl. The predictions for the semileptonic decay rates are displayed in Table |2^. 



Table 23. The calculated Ds K,K' transition form factors. My = Md' = 2.01 GeV , Mp ^ Md ^ 1.87 GeV. For the 



form factors F+, Ft,V, Ao,Ti the fit formula Eq. (4.11) is used, for the other form factors - Eq. ([1.12 





Ds^K 


Ds^K* 1 




F+ 


Fo 


Ft 


V 


Ao 


Ai 


A2 


Ti 


T2 


T3 


/(O) 


0.72 


0.72 




0.77 


1.04 


0.67 


0.57 


0.42 




0.71 


0.71 


0.45 




0.20 


0.41 




0.24 


0.24 


0.20 


0.29 


0.58 




0.22 


-0.06 


1.08 


0-2 




0.70 










0.42 








0.44 


0.68 



Table 24. The Ds {K, K*)lu decay rates in 10^° s"\ \Vcd\ = 0.22. 



Ref 


V{Ds ^ K) 


T{Ds -» K*) 


V{K'')/V{K) 


Tl/Vt 


This work 


0.63 


0.38 


0.6 


1.21 



For the coupling constants we obtain 



AD') 



2Md' 



0.95 ±0.05, 



AD) 

SD^^^'fp ^1.85 ±0.15, 

2Mk- 



AD') 

JT 

AD') 

Jv 

QD'DsK 



= 0.9 ±0.1, 



1.15 ±0.15. 
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2. ^r,,r,',<l> 



These decay modes are induced by the charged current c — > s quark transition. The pseudoscalar mesons rj and rj' 
are mixtures of the nonstrange and the strange components with the flavour wave functions rjn = "^^"^ and rjs — ss, 
respectively, 



77 = cos{ip) 77„ - sin(v5) ij^ 
rj = sin(.^) ry„ + cos(v5) 77^ 



(4.17) 



with the angle (/? ~ 40° Hg]. 
The decay rates of interest are 



* rjlv) ^ fi\-n^{ip)T{Ds 



(4.18) 



Let us give a brief explanation of these formulas: The semileptonic decay rates are determined by the form factor 
The spectral representation of this form factor does not involve the final meson mass explicitly. This means that for 
the ss component of both rj and rj' we have to deal with the same form factor, which can be expressed through the 
radial wave function of this component. On the other hand, the phase-space volume of the decay process is determined 
by the physical meson masses, as indicated in (4.18). It should be clear, however, that the rjs is not an eigenstate of 
the Hamiltonian and does not have a definite mass. 

Assuming universality of the wave functions of the ground state pseudoscalar 0^ nonet, the radial wave function 
of the nonstrange component ^^n^ coincides with the pion radial wave function | ]69[ |. The radial wave function 
should be determined independently. From the analysis of a broad set of processes the leptonic decay constant /s of 
the strange component ry^, has been found to lie in the interval fs = (1.36 ± 0.04)/^ [Q. This allows us to determine 
the slope parameter in such a way that the calculated value of fs lies in this interval, and the calculated ratio 
T{Dg rj)lT{Dg rj') agrees with the experimental data for — 40°. This procedure yields for the slope parameter 
(3^1^ — 0.45 Q For the slope parameter of the wave function of the (/)-meson, which is the vector ss state, we expect 
a value close to /?,,^ . 

In fact, — 0.45 GeV leads to the Bg ^ (/> transition form factors which agree well with the LCSR results at 
— (see subsection 4) . With all other quark model parameters fixed from the description of the nonstrange heavy 
meson decays and by taking a simple Gaussian form of the radial wave function, the decay rate T{Ds (plu) is a 
function of This function has a minimum at the value (3^ = 0.45 GeV; nevertheless, the corresponding value of 
the decay rate is la above the central experimental value. 
The results of our calculations are given in Tables Eq and Eq. 



rjs,(l) transition form factors. My ~ Af 



D* 



Table 25. The calculated D, 
form factors F+, Ft, V, Aq, Ti the fit formula Eq. (4.11 ) is used, for the other form factors 



2.11 GeV, Mp 
Eq. 



Mr 



1.97 GeV. For the 



(4.12) 







Ds - 




Ds - 




-x^ 1 




F+ 


Fo 


Ft 


Fo 


Ft 


V 


Ao 




A2 


Ti 


T2 


Ts 


/(O) 


0.78 


0.78 


0.80 


0.78 




0.94 


1.10 


0.73 


0.64 


0.47 


0.77 


0.77 


0.46 


CTl 


0.23 


0.33 


0.24 


0.21 




0.24 


0.26 


0.10 


0.29 


0.63 


0.25 


0.02 


1.34 


0-2 




0.38 




0.76 
















2.01 


0.45 



^Another procedure of taking into account the SU(3) breaking effects to obtain 5',,^ from 'I',,^ has been proposed in [69|. 
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Table 26. The Ds {rj,r]' ,(l))lv decay rates in 10^" s ^, |Vcs| — 0.975. The experimental rates are obtained from the 
corresponding branching ratios using the Ds lifetime td, ~ 0.495 ± 0.013 ps from the 1999 update [p4| 



Ref 


T{Ds ^ ,7) 


V{Ds ^ r^') 


r(D, ^ <j>) 


This work 
Exp |9| 


5.0 
5.2±1.3 


1.85 
2.0±0.8 


5.1 

4.04±1.01 



For the coupling constants we obtain 



9D',Dsris Jy 

2Md- 

ADs) 
gD,D,(liJp 



2Ah, 



ad:) 

JT 

ADi) 

Jv 
9d^Ds<p 



= l.OiO.l, 

1.6 ±0.3, 

0.93 ±0.03, 
1.08 ±0.04. 
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3. Bs K, K* 



This mode is driven by the b ^ u charged current tr ansiti on. The only additional new parameter needed here is the 
slope of the Bs wave function. We obtain it by using (3.63) and taking fsj fs = l-l- The results of our calculation 
are given m Table 



Table 27. The calculated B. 
form factors _F- 



K* transition form factors. 



Mv = Mb' = 5.32 GeV, Mp = Mb = 5.27 GeV . For the 
Ft,V,Aq,Ti the fit formula Eq. (4.11) is used, for the other form factors - Eq. (4.12). 





Bs-^K 


Bs K* 




F+ 


Fq 


Ft 


V 




^1 




Ti 


Ta 


Ti 


/(O) 


0.31 


0.31 




0.31 


0.38 


0.37 


0.29 


0.26 




0.32 


0.32 


0.23 




0.63 


0.93 




0.61 


0.66 


0.60 


0.86 


1.32 




0.66 


0.98 


1.42 


0"2 


0.33 


0.70 




0.30 


0.30 


0.16 


0.60 


0.54 




0.31 


0.90 


0.62 



These form factors lead to the following relations 

AB') 

gB'B.KJv 



2Mb' 

AB) 
gSB^K'Jp 



2Mk' 

AB") 

JT 

AB') 

Jv 

gB'B^K- 
gB'B^K 



= 0.44 ±0.04, 
= 1.3±0.1, 
= 0.95 ±0.05, 
= 1.2 ±0.1. 



The form factors at g = are compared with the LCSR predictions in Table 28, We observe some disagreement 
between our predictions and the LCSR calculation which gives smaller values for all the form factors. A closer look 
at the origin of this discrepancy shows that its source is the strength and sign of the SU(3)-breaking effects. They 
lead to opposite corrections in the two approaches. 



Table 28. Comparison of the QM and LCSR results on the Bs K, K* form factors at = 0. 



Ref. 


vm 


Ai(0) 


^2(0) 


A)(0) 


Ti(0) 


T3(0) 


This work 


0.38 


0.29 


0.26 


0.37 


0.32 


0.23 


LCSR'98 |§ 


0.262 


0.19 


0.164 


0.254 


0.22 


0.16 



To discuss these SU(3) breaking effects, let us start with i? — > p, which in fact differs from the Bg K* only by 
the flavour of the spectator quark, and move to Bg K* by accounting for the SU(3) violating effects: 

Within the LCSR method there are two changes which affect the form factors: first, the change fs^ Jb leads 
to an increase of the Bg — > K* form factors; second, the change of the symmetric twist-two distribution amplitude of 
the p-meson to the asymmetric one of the K* meson leads to a decrease of the form factors. The second effect turns 
out to be much stronger than the first one with the result of an overall decrease of the form factors. 

In the quark model, the same SU(3) breaking effects take place: The change of the spectator mass (it determines 
the increase of /bJIb) and the change of the K* meson wave function (due to the change of both the quark mass 
and the slope parameter of the light meson wave function). Here the influence of the slope of the heavy meson wave 
function upon the form factor is only marginal. Therefore, the resulting effect of these changes leads to an increase 
of the form factors. 

We want to point out that the difference between the results of the two approaches does not arise from specific 
effects (higher twists, higher radiative corrections etc) which are present in the LCSRs but absent in the quark model. 
The observed difference is only due to the different strength of the SU(3) violating effects at the level of the twist-2 
distribution amplitude. As was discussed in , this distribution amplitude can be expressed through the radial soft 
wave function of the meson. The change of the quark-model wave function caused by SU(3) violating effects does not 
induce a strong asymmetry in the leading twist-2 distribution amplitude. 

In view of the discrepancy between our results and the LCSR it would be interesting to have independent calculations 
of the Bg K* form factors at small from the 3-point sum rules, as well as a lattice calculation for large q^. 
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4- Bs v,v',<t> 



These weak meson transitions are induced by the FCNC b ^ s quark transition. The resuhs of the form factor 
calculation are given in Table ^ and compared with the LCSR predictions at = in Table ^ The agreement 
between the two values is satisfactory at least within the declared 15% accuracy of the LCSR predictions. This allows 
us to expect that also the Dg —>■ 4>, ?7, v' form factors and the corresponding decay rates (given earlier in subsection 2) 
are calculated reliably. 



Table 29. The calculated Bs ^ ris,(j) transition form factors. My ~ Mb* = 5.42 GeV, 
form factors F+, Ft ,V, Ao,Ti the fit formula Eq. (4.11) is used, for the other form factors - 



Mp = Mb, = 5.37 GeV. For the 
Eq. (p2|). 







Bs risiM,,) 


Bs - 


-^rUM^,) 


Bs^<t> 1 




F+ 


Fo 


Ft 


Fo 


Ft 


V 


Ao 


^1 


A2 


Ti 


Ta 




/(O) 


0.36 


0.36 


0.36 


0.36 




0.39 


0.44 


0.42 


0.34 


0.31 


0.38 


0.38 


0.26 


CTl 


0.60 


0.80 


0.58 


0.80 




0.58 


0.62 


0.55 


0.73 


1.30 


0.62 


0.83 


1.41 


0-2 


0.20 


0.40 


0.18 


0.45 




0.18 


0.20 


0.12 


0.42 


0.52 


0.20 


0.71 


0.57 



Table 30. Comparison of the QM and LCSR results on the Bs (j> form factors at = 0. 



Ref. 


V{0) 


Ai(0) 


^2(0) 


^o(O) 


ri(o) 


T3(0) 


This work 


0.44 


0.34 


0.31 


0.42 


0.38 


0.26 


LCSR'98 |§ 


0.433 


0.296 


0.255 


0.382 


0.35 


0.25 



0.6 ±0.05, 

1.5 ±0.1, 

0.95 ±0.05, 
1.13 ±0.06. 



For the couphng constants we obtam 



9B'B,rts Jy 

2Mb* 

gB,B,4,Jp 

ABl) 
It 
ABt) 
Jv 

9B'B,tj> 
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E. Discussion 



In this Chapter we have calculated numerous form factors of heavy meson transitions to light mesons which are 
relevant for the scmilcptonic (charged current) and penguin (flavor-changing neutral current) decay processes. Our 
approach is based on evaluating the triangular decay graph within a relativistic quark model which has the correct 
analytical properties and satisfies all known general requirements of long-distance QCD. 

The model connects different decay channels in a unique way and gives the form factors for all relevant values. 
The disadvantage of the constituent quark model connected with its dependence on ill-defined parameters such as the 
effective constituent quark masses, have been reduced by using several constraints: the quark masses and the slope 
parameters of the wave functions are chosen such that the calculated form factors reproduce the lattice results for the 
B — > p form factors at large and the observed integrated rates of the semileptonic D ^ K, K* decays. 

Our main results from this Chapter are as follows: 

• In spite of the rather different masses and properties of mesons involved in weak transitions, all existing data 
on the form factors, both from theory and experiment, can be understood in our quark picture. Namely, all 
the form factors are essentially described by the few degrees of freedom of constituent quarks, i.e. their wave 
functions and their effective masses. Details of the soft wave functions are not crucial; only the spatial extention 
of these wave functions of order of the confinement scale is important. In other words, only the meson radii are 
essential. 

• The calculated transition form factors are in all cases in good agreement with the results available from lattice 

QCD and from sum rules in their specific regions of validity. The only exception is a disagreement with the 
LCSR results for the K* transition where we predict larger form factors. This disagreement is caused 

by a different way of taking into account the SU(3) violating effects when going from B — > p to Bs — > K* 

and is not related to specific details of the dynamics of the decay process. We suspect that the LCSR method 
overestimates the SU(3) breaking in the long-distance region but this problem deserves further clarification. 

• We have estimated the products of the meson weak and strong coupling constants by using the fit formulas for 
the form factors for the extrapolation to the meson pole. The error of such estimates connected with the errors 
in the extrapolation procedure is found to be around 5-10%. 

We cannot provide for definite error estimates of our predictions for the form factors because of the approximate 
character of the constituent quark model. However from the fine agreement obtained in cases where checks are 
possible, we believe that the actual accuracy of our predictions for the form factors is around 10%. Since some 
parameters have been fixed by using lattice results and have also been tested using the sirni rule predictions, further 
improvements of the accuracy of our predictions will follow if these approaches attain smaller errors. Of course, each 
precisely measured decay will also allow a more accurate determination of the parameters of the model and thus can 
be used to diminish the errors at least for closely related decays. 
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WEAK ANNIHILATION IN THE RARE RADIATIVE B 

FACTORS 



fry DECAY AND THE B "flv FORM 



Results obtained in the previous Chapter allow us to fully describe semileptonic i?(s) and decays. We now 
discuss rare radiative decays of heavy mesons induced by the flavour-changing neutral currents (FCNC) such as 
B 
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K*^ and B ^ decays. This Chapter is based on Refs. [|64 
The investigation of rare semileptonic B decays induced by the flavour-changing neutral current transitions b s 
and b ^ d represents an important test of the standard model or its extentions. Rare decays are forbidden at the tree 
level in the standard model and occur at the lowest order only through loop diagrams. This fact opens the possibility 
to probe the structure of the electroweak sector at large mass scales due to contributions of virtual particles in the 
loop diagrams. Interesting information about the structure of the theory is contained in the Wilson coefRcents in the 
effective Hamiltonian which describes the 6 — > s, d transition at low ener gies . These Wilson coefficients take different 
values in different theories with testable consequences in rare B decays [ 100| -107|. 

Among rare B decays the radiative decays b sj and b dj have the largest probabilites. The 6 — s- 57 transitions 
are CKM favoured and have larger branching ratios than the b —^_drf transitions. The b sj transitio n ha s been 
observed by CLEO in the exclusive channel B K*j in 1993 |108| and measured inclusively in 1995 [109|. The 
B ^ pj decay will be extensively studied by BaBar and BELLE. 

Any reliable extraction of the perturbative (short-distance) effects encoded in the Wilson coefficients of the effective 
Hamiltonian requires an accurate separation of the nonperturbative (long-distance) contributions, which therefore 
should be known with high accuracy. 

In rare exclusive semileptonic decays one faces several types of long-distance contributions which arise in the meson 
transition amplitude of the effective Hamiltonian, which contains penguin and 4-quark operators (see |110| for more 
details). 

The dominant long-distance effects for the radiative decays arise from the electromagnetic penguin operator and is 
given in terms of the form factor Ti{q^ — 0). 

The 4-quark operators generate many different contributions: 

One of them is the so-called long-range charming penguins, i.e. intermediate cc states in the g^-channel which 
include the cc resonances (■0, V'', ■••) and the charmed hadronic continuum. The charming penguins are crucial for 
rare semileptonic decays b {s,d)l'^l~ , but provide only a small contribution for the radiative decays b — > ( s, d)y 
[ 111 - 113 |: In the factorization approximation the cc contribution precisely vanishes due to the gauge invariance [ 113 |, 
whereas the nonfactorizable effects turn out to be only at the level of few percent [112|. 

Another effect of the 4-quark operators is the so-called weak annihilation which corresponds to the annihilation of 
quarks in the channel of the initial or final mesons. 

In the B K*"f decay the weak annihilation is negligible compared to the penguin effect: it is suppressed by two 
powers of the small parameter A ~ 0.2 of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. In _B ^ p7 both effects 
have the same order in A and must be taken into account. 

The penguin contribution is known rather well. On the other h and^ the WA in B py has been studied in less 
detail: the relevant form factors were analysed within sum rules |114,115| and perturbative QCD |11(;]. However, 
some contributions to these form factors were neglected. These contributions may be relevant if precise measurements 
become available. 

We now apply the dispersion approach to the weak annihilation in the B^ — > p^j decay. 

In the factorization approximation |117| the weak-annihilation amplitude can be represented as the product of meson 
leptonic decay constants and matrix elements of the weak current between meson and photon. The latter contain 
the meson-photon transition form factors and contact terms which are determined by the equations of motion. The 
photon can be emitted from the loop containing the b quark which is described by Bj transition form factors. It can 
also be emitted from the loop containing only light quarks described by the pj transition form factors (Fig Gu) . 





Fig. 20. Diagrams describing the weak annihilation process for B ^ p7 in the factorization approximation: (a) The photon 
is emitted from the loop containing the b quark, (b) The photon is emitted from the loop containing only light quarks. 
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We calculate the -B7 form factors within our relativistic dispersion approach which expresses these form factors in 
terms of the B meson wave function. We demonstrate that the form factors calculated by the dispersion approach 
behave in the limit nii, ^ 00 in agreement with perturbative QCD. The B meson wave function was previously tested 
in the B — > light meson weak decays and is known quite well, allowing us to provide reliable numerical estimates for 
the Bj form factors. 

The /97 transition form factor is related to the divergence of the vector and axial- vector currents. In the case of 
the axial-vector current it is proportional to the light-quark masses if the classical equation of motion is applied. 
Because the quark momenta in the loop are high, these masses have to be identified with current quark masses. For 



this reason the corresponding p7 form factor was neglected in previous analyses |114,115|. We find however that this 
argument is not correct: this form factor remains finite in the limit of vanishing light quark mass m and behaves 
like ~ Mpfp/Mg which means the violation of the classical equation of motion. Here we present the result for the 
P7 transition form factor. A detailed discussion of the anomaly appearing for the matrix element {p^\d^A^\0) can be 
found in Ref. |6| 

Finally, we provide numerical estimates of the weak-annihilation amplitude taking into account the Bj, transition 
form factors, and the contact term contributions. 



In section VA the effec tive Hamiltonian for the b ^ d transition and the general structure of the amplitude are 



presented. In Section VB we discuss the photon emission from the B meson loop and obtain the Bj transition form 



facto rs within the dispersion approach. Section V C contains resu lts for the p7 transition form factors. In Section 



V D the numerical estimates are given. The concluding Section V E summarises our results for the weak annihilation. 



A. The effective Hamiltonian, the amplitude and the decay rate 

The amplitude of the weak radiative i? — s- p transition is given by the matrix element of the effective Hamiltonian 
for the b ^ d transition 

AiB ^ p7) = (7((zi)p(g2)|i/cff (fo - d)\Bip)), (5.1) 

where p is the B momentum, 52 is the p momentum, and qi is the photon momentum, p ~ qi + q2, qf = 0, = M^, 
p^ = Mg. The effective weak Hamiltonian has the structure [ 11§| -T21|: 



H,s{b ^ rf) - ^ ^tC7^{fi)Or^ - {Ci{p)Oi + ^2(^)02) , (5.2) 

where only operators relevant for the penguin and weak ahnnihilation effects are listed. Gp is the Fermi constant, 
— V*j^Vqb, CiS are the Wilson coefficients and O^'s are the basis operators 

g _ 

= T-^ daapumb{p) (1 + 75) ba Fpu, (5.3) 

01 = da^i,{l - 75)Wa M/37i'(l " 75)&/3, 

02 = da^i,{l - lb)u[j ?2/37i/(l - 75)&a- (5.4) 

with e — ^JAnacra- Other notations are given in Eq. ( |3.7D . 
The amplitude can be parametrized as follows 

A{B^ p"7) = [egie^^aeJ-Fpc + i^l" {gv^pqi - P^J,qlu) -Fpv] , (5.5) 

where Fpc and _Fpv are the parity-conserving and parity- violating invariant amplitudes, respectively. £2(61) is the 
yO-meson (photon) polarization vector. We use the short-hand notation Cabcd — ^aPp.iyCL^b^c'^d'^ for any 4-vectors 
a, 6, c, d. 

For the decay rate one finds 

TiB- ^ p-7) = %2!!iM| (l - M,Vm|)^ (|FpcP + iFpvP) . (5.6) 
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1. The penguin amplitude 

The main contribution to the amphtude is given by the electromagnetic penguin operator Or-yi 

^pcng(S -> p7) -^^^tCr^TiiO) (egic-g^e; + iel''e*/{g^f,pqi - Pf,qi^)) , (5.7) 

where Ti is the form factor of the B p transition through the tensor curr ent 0j2|,|^,|l|. The corresponding 
contribution to the invariant amphtudes is therefore 

F^c' = F^v' = -^tCr^rm- (5.8) 



2. The weak annihilation amplitude 

The radiative B ^ p"/ transition also receives contribution from the 4-fermion operators Oi and 02- For the 
charged B^ — > /3^((j'2)7('?i) transition the corresponding amplitude reads 

AwAiB- ^ p--f) = -^^„(p(g2)7(qi)|J7,(l - 75)u . U7,(l - l5)b\Bip)), (5.9) 

In what follows we suppress the lable WA. Neglectin g th e nonfactorizable soft-gluon exchanges, i.e. assuming vacuum 
saturation, the complicated matrix element in Eq. (^.9|) is reduced to simpler quantities - the meson-photon matrix 
elements of the bilinear quark currents and the meson decay constants. The latter are defined as usual 

{piq2)\dj,u\0) ^ e;^Mpfp, fp>0, 
(0|w7,756|B(p)) = ip.fB, fs > 0. (5.10) 

It is convenient to isolate the parity-conserving contribution which emerges from the product of the two equal-parity 
currents, and the parity-violating contribution which emerges from the product of the two opposite-parity currents. 

The parity-violating contribution 

The parity-violating contribution to the weak annihilation amplitude has the form 

ApyiB ^ Pi) = {(p7M7."|0)(0|u7.75fo|S) + {p\d-,,u\0){-/\ui,i5b\B)} . (5.11) 

Here ai is an effective Wilson coefficient, which we take as ai = C\ + C2/NC at the scale ~ 5 GeV. 
It is convenient to denote 

AV = {p-{q2)\diuu\G){i{qi)\ui,i5b\B-{p)) (5.12) 

and 

AT = {p-{q2)l{qi)\diM^){Q\u-i,i^h\B-{p)). (5.13) 
1. Let us start with . One can write 

{l{qi)\ui.l^h\B- {p)) = ee^T^f, (5.14) 

where 



T^Ap.Ii) = ^ / dxe'^^^{Q\T{J^(x),ui,i^b)\B-{p)), (5.15) 

and 



2 _L - i- 

Jt^{^) = -^u-f^u - -d'^f.d - 2^7m^ (5.16) 



75 



is the electromagnetic quark current. 

The amphtude in general contains 5 independent Lorentz structures and can be parametrised in various ways. 
A possible way is to write as follows 



rr-iB rTi_L 



pqi 



pqi 



C2, 



(5.17) 



where 



+ i (qiPt^ ~ pqi qip) qivF2A{ql) 



is transverse with respect to gi^, 



~Wi qit?jPyF^A{q 



(5.18) 



(5.19) 



The Lorentz structures in this expansion are chosen such that they contain no singularity for q\ ^ Q. 

The invariant amplitudes Ci and C2 in the longitudinal structure can be determined using the conservation of the 
electromagnetic current 9^ = 0, which leads to the relation 



qit^Tj^Mqi) = -(0|[Q,d>75"]|5~(p)) 

= iQeIbPu- 



This gives 



Ci — QbJb, 
C2 = 0. 



(5.20) 



(5.21) 



In particular, C = —fs for the B~ meson. Notice however, that the longitudinal part of the amplitude in this form 
is not a contact term^ This means that some of the invariant amplitudes in the transverse part are mixtures of the 
true form factors (i.e. quantities with definite and known analytic properties) and remnants of the contact terms. 

After setting q\ = and multiplying by e*'^ (this very combination is needed for A^^) only some of the Lorentz 
structures contribute and we find 



^1 J-t,Jqf=o 



Ve 



Mj 



(5.22) 



where we have used the relatio n qip — ^(A/^ — M^) for q\ — 0. 

The parametrization of Eq. (5.17) is of course not unique and there are other possibilities. There is however only 
one choice when the longitudinal part is a pure contact term. Such a choice, which is also the most natural one, is 
prompted by the structure of the Feynman diagram: 

Let us rewrite the usual electromagnetic coulping of the quark as follows {qi = k' — k): 



(m + k')-f^{m + k)^{ni + k') I 7^ 



. ql^i 
qi — 

qi 



(m + k) 



ql 



{k'^ - vi^){m + k)- {k^ - w?){m + k') 



The first term is explicitly transverse with respect to qi^ and leads to T^. 

The second term, containing the factors (fc^ — m?) and (fc'^ — m^), leads to the contact term 

Wu — JB- 

ql 



(5.23) 



^We recall that according to the usual definition, a contact term is a quantity which is represented by the delta-function and 
its dirivatives in the coordinate space. One can then easily check, that the Fourier transform of the expression -^^^ is not a 
contact term. 
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Finally, we come to the following parametrization 



^fiv — ^iiu ^ ~2 ^1 "I l2 '-'21 



9i 



9? 



(5.24) 



with 



9? 



9iM j 9i^^2a(9i) 
giM )p^^-P'3A(9l)• 



Notice, that this is the only parametrization of the amplitude, which provides a distinct separation of the amplitude: 
form factors in the gauge-invariant transverse part of the amplitude, and contact terms in the longitudinal part of 
the amplitude. 

Clearly, = Fia, whereas other form factors of the sets F and F' are different. Now the Lorentz structures 
contain explicit singularities for qf = which must cancel each other in the amplitude T^. This leads to the 
constraints on the form factors at = 



Equation (5.20) gives 



and hence 



qipFUO)-qipFiAO) + C!, = 0, 
-qipF^AO) + C[=0. 



^2=0, 



(5.26) 



(5.27) 



i^i'^(0) = -F^^(0), 



2,/s 



(5.28) 



By virtue of the relations (5.26), for the amplitude at = we find 

A'C = lefpMpel^e*/ {g^.TOiFi^(O) + ^^^1.^2^(0) - ^^^.^3^(0)} 



ie ^pMpt*^t*^ \ (g^^pqi - Pf,qiu)FiA{0) - p^^qiu 



Ve 



Mi 



Af2 



(5.29) 



Notice that the contact term does not contribute to the amplitude directly, but indirectly determines the value of the 
form factor ^3^(0). 



2. Let us now turn to Ajp^ . 

Using the equation of motion for the quark fields {Qu = 2/3 e, Qd ~ 

i^vd^qix) = mq{x) - QqA^^''q{x), 
id''q{x)^y = ~mq{x) + QqAyq{x)^'" 



-1/3 e) 



(5.30) 



one obtains for ^ 



^2 ^iPi^fsip j\d-f^u\0) 
^ feip-imdluum 
= -{Qd-Qu)e{p-j\A''d-f,u\0) 
= ee]^e*2''gpyfpMpfB + 0(m„, m^). 



(5.31) 



Now the full contribution is due to the contact terms only, since the transverse form-factor part of the amplitude 
(p~7|(i7i/w|0) vanishes after the multiplication by Pi,. 
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3. For the sum we find 



PV I ^PV 



Ve 



Ml 



M2 



(5.32) 



Summing the contributions of the photon emission from the i?-meson loop and the p-meson loop gives the amplitude 
Apv which can be represented in the form Apy = e*^Tp with q^T^ = as required by gauge invariance. The 
weak-annihilation contribution to the form factor -Fpv for the — > decay takes the form (after redefining 
F^A = 2Fa/Mb) 



2Fa 
Mb 



Ml Af 2 



(5.33) 



Notice that in addition to the form factor contribution Fa, there is a contribution proportional to fs- The latter 
is determined by the contact terms: namely, the contact terms which are present in the amplitudes of the photon 
emission from the B meson loop and from the p-meson loop do not cancel each other, but lead to a nonvanishing 
contribution]^ In the amplitude of the B^ —> decay the contact term proportional to fs changes sign, and also 
= —F^ as can be obtained from charge conjugation of the amplitude ^pv- For the B^ — > p^j decay the 
contact term proportional to fs is absent. 

The parity- conserving amplitude 

This amplitude reads 



ApciB ^ Pi) = -^Cuai{{p\d^,u\0){j\u-/,b\B) + (7pM>75u|0)(0|jl7,75&|B)} 



(5.34) 



The _B ^ 7 amplitude from the first term in the brackets contains the form factor Fyiq^ 



Ml 



(7(gi)|u7.6|i3(p)) 



e e 



Mb ' 



(5.35) 



The second term in ( ^.34 ) can be reduced to the divergence of the axial- vector current and contains another form 
factor, Gv- namely. 



'|u7^756|B)(7p|(i7^ 75^10) = fB{lp\d,yd'^^j5u\Q) 

= e fBGvf^q^alqze^- 



(5.36) 



Thus, the weak annihilation contribution to Fpc reads 



-^pc ~ iu^iMpfp 



2Fv 
Mh 



Mpfp 



(5.37) 



Summing up, within the factorization approximation the weak annihilation amplitude can be expressed in terms of 
the three form factors Fa, Fy, and Gy- 



^ Here we disagree with the r ecen t statement of [122] that contact terms cancel each other and do not contribute to the final 
amplitude. The conclusion of |122 is based on choosing the longitudinal part of the amplitude in the form |(?^^ 92^92,/ 
which is not a contact term according to the usual definition of the latter. As a result, the invariant amplitude of 



at gi = corresponds to the combination Fa + 

B 

results from different approaches later in this section. 



We take this relationship into account when comparing numerical 
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B. Form factors for the B — > ^Ev transition 



In this section we derive formulas for the form factors Fa,v describing the weak transition B —>■ ^Iv within the 
dispersion approach to the transition form factors. Recall that the form factors FAy describe the transition of 
the B-meson to the photon with the momentum gi, q\ = 0, induced by the axial- vector (vector) current with the 
momentum transfer (72, q| = M^. For technical reasons, it is convenient to treat the form factor F^iy) ^ describing 
the amplitude of the photon-induced transition of the _B-meson into a hu axial-vector (vector) virtual particle with 
the corresponding factor l/(s — g^) the dispersion integral. Then we can directly apply the equations obtained in 
Chapter II for meson-to- meson transition form factors. 



1, The form factor Fa 



The form factor Fa is given by the diagrams of Fig |2l|. FigJ 
the process when the h quark interacts with the photon; Fig \ 
quark u interacts while h remains a spectator. One has 



shows f'^\ the contribution to the form factor of 
describes the contribution of the process when the 



Fa^F^ 



F 



(u) 



(5.38) 





= (-) 




^TvTs 



a. b. 

Fig. 21. Diagrams for the form factor Fa- a) F^\ b) F^ 



s — Mg ) . The B-meson wave 



The meson is described by the vertex 6(fcb) ij^uj ku) G (s) / \/Wc, with G{s) — (fB{s){s — Mg 
function tps is normalized according to the relation ( 3.62| ). 

It is convenient to change the direction of the quark line in the loop diagram of Fig pTjj. This is done by performing 
the charge conjugation of the matrix element and leads to a sign change for the 71/75 vertex. 

(1)/ 



Now both diagrams in Fig Ell a,b are reduced to the diagram of Fig which defines the form factor Fj^ (mi , 7712) 



Setting mi = to;,, m2 = to^ gives F 



F 



(b) _ 



QbF'^\mb,mu). 



Similarly, setting mi — to„, m2 = m^ gives F)^ , 

F^A^ ^-QuF['\mu,mb). 




Fig. 22. The triangle diagram for (mi, 7712). The cut corresponds to calculating the imaginary part in the variable p^. 
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For the diagram of Fig pa (quark 1 emits the photon, quark 2 is the spectator) the trace reads 



-Sp ij5{m2 ~ k2)-f„j5{mi + k[)^^{mi + ki) 

= 4i(fci + k[)f^{mik2 + ^2^1)^ + Ai{g^^qa - g^aq,^)imik2 + m2ki)c 



(5.39) 



The spectral density of the form factor i<"^^-'(mi,TO2) in the variable p^, p = fci +fc2, is the coefficient of the structure 
g^i, obtained after the integration of the trace over the quark phase space. Performing necessary calculations, we 
arrive at the following single dispersion integral 



Mb 



-F 



(1) _ 



47r2 



mi log 



ds ipb{s) 
{s ~ Ml) 



m| + A^/2(s, 771^,7712) 



+ (m2 - mi) 



Ai/2(s,mg,m2) 



1 



(Ml 



2s 
M2)/2. 



Ai/2(s, nrf, m^) 

s + m\ — m\-\- A-^/2(s, m\,m\) 



nil log 



s + mi — m2 — A-'^/2(s, mi, rn^) 



(5.40) 



For qf = one has pqi — yiv-i^ ^-^p) 

Now, let us analyse the behaviour of the form factor in the limit m^ — > cx). To this end it is convenient to rewrite 
the spectral representation ( ^.40 ) in terms of the light-cone variables as follows (see |^ for details) 



2 

Mb 



-^F^2\mi,m2) 



47r2 



dxidx2dk\ 

x\x2 



5{l — xi — X2) 



S-M2 



X \^miX2 + m2Xi — (mi — m2)fc^/p(7i} . 
Here Xi is the fraction of the _B- meson light-cone momentum carried by the quark i, and 

s = m\/xi + m\/x2 + k\/xiX2- 



For the form factors f'^^ and F^^' one obtains 



Mb 
2 
M~B 



-F 



Qu 
Qb- 



Nc f dxdk\ fBi.s) 



47r2 
47r2 



X^Xh s - M2 
dxdk^ (Pb{s) 
xlxu s~ Mj 



ruu Xb + mb Xu + 



■nib Xu + muXb + 



2{mb - rnu)k']_ 



Ml 



M2 



2(m„ - m{,)fcj 



Ml 



M2 



with 



(5.41) 



Xb 



•^11. 



(5.42) 



Recall that the i?- meson decay constant is given in terms of the wave function by Eq. ( 3.63 ). Due to the wave 
function ipb{s), the integral in the heavy quark limit is dominated by the region Xu = ^/rnb^ Xb ^ ^ — A/rub, where 
A is a constant of order Mb — mb- This leads to the following expansion of the form factors in the l/m^ series 



(u) 



Kruh 



Mf 



-F 



(b) _ 



n 

Wb — 2 



(5.43) 



Clearly, the dominant contribution in the heavy quark limit comes from the process when the light quark emits the 
photon, whereas the emission of the photon from the heavy quark gives only a 1/m;, correction. The expressions 
(5.43) for the form factor F^' agrees with the result of Ref. |116|, while we have found a different sign forF^'). 
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2. The form factor Fv 



The consideration of the form factor Fy is very similar to the form factor Fa ■ Fy is determined by the two diagrams 



shown in Fig 23: Fig |23| a gives Fy \ the contribution of the process when the b quark interacts with the photon; Fig 
describes the contribution of the process when the quark u interacts. One has 



(&) 



(5.44) 




a. b. 

Fig. 23. Diagrams for the form factor Fv- a) -F^?\ b) F^\ 

It is again convenient to change the direction of the quark hne in the foop diagram of Fig Hb describing i^^"^ by 
performing the charge conjugation of the matrix element. For the vector current 7,^ in the vertex the sign does not 
change (in contrast to the 71^75 case considered above). 

Then both diagrams in Fig ^ a, b are reduced to the diagram of Fig |^ which gives the form factor Fy\mi,m2)- 
Setting mi — rub, m2 — ruu gives Fy'': 

F^^^ = QbFl,'\mb,m^). 

Setting mi = m„, m2 = rrib gives Fy'^ 

fI;^ = QuFl,^\mu,mi,). 




Fig. 24. The triangle diagram for Fy\ini,m2). The cut corresponds to calculating the imaginary part in the variable p^. 



The trace corresponding to the diagram of Fig 4(1- active quark, 2 - spectator) reads 

- Sp i-f5im2 - hhi^imi + k[)jf_,{mi + ki) = -4e^^agi ("^1^2 + "^2^1)0- 
The spectral representation for the form factor takes the form 



Mb 



-F, 



(1) 



47r2 



ds ipBjs) 



S \S + TOj — to| — Al/2(s, TOi, TO2) 



(5.45) 
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To analyse the heavy quark hmit rrih — > oo we again represent the form factor in terms of the light-cone variables 



Mb' 



x\x2 



■S{l-xi 



Ml 



(5.46) 



For Fy and Fy the corresponding expressions read 



Mb ^ 



(b) 



Nc f dxdk\ if sis) 



47r2 
47r2 



dxdk^ ^b{s) 



{rUu Xb + mb Xu} ; 



Ml 



{mbXu +mu Xb} , 



By the same procedure as used for Fa, we obtain in the limit rub 00 



Mb 



-F, 



(«) 



Anib 



n r y ~ 2 



(5.47) 



The dominant contribution in the heavy quark limit again comes from the process when the light quark emits the 
photon. Now both fo rm fa ctors F,^/^ and Fy'^ in ( 5.47 ) perfectly agree with the expansions obtained in 

As seen from Eqs. (5.43) and (5.47), one finds Fa — Fy in the heavy quark limit, in agreement with the large-energy 
effective theory 
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C. The form factor Gv 



The form factor Gv is determined by the divergence of the matrix element of the charged current between the 
vacuum and the p^^ state, 



The corresponding diagrams are shown in Fig Eq. 



(5.48) 





Fig. 25. Diagrams describing the amplitude (7(gi)p~(g2)|d7,,75g|0), p = qi + g2, p'^ ~ Mg. The diagram (a) is multiplied 
by the d-quark charge, and the diagram (b) is multiplied by the u-quark charge. The cut corresponds to the calculation of the 
imaginary part in the variable 



If the classical equations of motion are applied, the form factor Gy is proportional to the light-quark masses, 
is the reason why this form factor was neglected in previous analyses 



This 



110|. However, a proper calculation shows 



that this argument is not correct: in fact the classical equations of motions do not hold and the divergence contains 
the anomaly. 

The anomalous behavior of the divergence of the axial- vector current in the chiral limit is a well-known phenomenon 
discovered in the two-photon amplitude (77|i9i/97i/75<?|0) |123|. A very clear way to demonstrate the anomaly is to 
start with the matrix element of the axial vector current and to calculate the spectral representations for the relev ant 
form factors. The anomaly is then obtained by performing the divergence at the final stage of the calculation |124|. A 
similar treatment applied to the matrix element (7p~ |c?7j/75u|0) leads to the form factor Gv which does not vanish in 
the limit m — s- 0. We will present a detailed discussion of the anomaly in radiative B decays in a separate publication 
Jsst . Here we only quote the final result: considering the spectral representation in the variable g| and introducing 
the appropriate p-meson radial wave function 'ijjp{s), one finds the following expression for the form factor Gv 



M| 
47r2 



ds i,p{s) {s - Ml) 
(s - Ml - zO)2 



(5.49) 



i^pis) is normalized according to the relation (for massless quarks) 



oo 



= 1. 



(5.50) 



Clea rly, G y is finite for m = which means a violation of the classical equations of motion for the axial- vector current. 
Eq. ( 5.49| ) describes the anomaly which takes place for the 7p final state as well as for the 77 one. There is however 
an important difference between the two cases: For the 77 final state the divergence remains finite in the limit m ^ 
and = Ml — > 00. For the final state the divergence is finite for m ^ but decreases as l/p^ for 00. It is 

convenient to introduce the parameter n such that 



Gv = [Qu + Qd) K 



Mpfp 



M% 



(5.51) 



with K staying finite for m = and Mb 00. 
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D. Numerical estimates 



Table presents formulas for the penguin and the weak-annihilation invariant amplitudes for the p^^ and 

i?" (o'^Tdecays. 

Before giving numerical estimates, let us look at the scaling behavior of the form factors in the limit Mb oo: 



ri(0) - m/^^ [35], Fa^b - Mb'^': 



(5.52) 



such that 



^1/2 



pWA 71 .r- 



3/2 



(5.53) 



The terms proportional to Jb (Jb ~ l/V-^s) l/Afs-suppressed compared with the terms containing the form 
factors Fa,v- As we see below numerically this leads to a suppression by a factor of 4 — 5. 



Table 31. The penguin and the weak-annihilation invariant amplitudes for the B — > p 7 and B°{bd) — » ^"7 decays. 



Quantity 



B- 



P 7 



B° 



r^WA 
-f^PV 



iua\Mpfp 
^uaiMpfp 



2F': 



Mb 
Mb 



(0) 

■ + 



2fB 
Ml-M'i 



{Qu + Qd) K 



fB 



^ua2MpfpO J -jg. 

^ua2MpfpO ( _ 2 Q„ 



We now proceed to numerical estimates for the i3- meson decay. The s cale- dependent Wilson coefficients Ci (fi) and 
ai [ji) take the following values at the renormalization scale ~ 5 GeV ]l2l| : 



Ci==l.l, C2 = -0.241, C7^ = -0.312, ai = Ci + Ca/A^c =i 1-02. a2 ^ C2 + Ci/N^ ~ ±{).l. (5.54) 



The penguin form factor was previously calculated within the dispersion approach with the result (0) = 

0.27 ± 0.3, see Table ^ in Chapter IV. Using the same parameters and the B meson wave function as in Chapter IV, 

we obtain the form factors Fa,b shown in Table |3^. Our result for the form factor Fy is in good agreement with the 

estimates from other approaches. The form factor Fa agrees well with the constraints from perturbative QCD and 

turns out to be considerably larger than the corresponding sum rule estimate. 

The value of the Gy is sensitive to the details of the p meson wave function i/jp. The reason for that is the 

t2\ 



tne p 

presence of the term (s — M^) in the integrand in Eq. ( 5.49| ) which changes sign in the integration region. Assuming 
'tpp{s) ~ /3^/(/3^ -I- s)^ and setting (3 = 0.8 GeV, which gives a good description of the p meson radius, leads to 
K, — —1.8. Conservatively, we take —k < 2.0 and use this result for further estimates. The form factor Gy does not 
contribute more than a few % to the full amplitude, but can sizeably correct the weak-annihilation part, especially 
for the B'^ —> p'^j decay. 



Table 32. The weak-annih ilatio n form factors Fa, Fy and Gv- The accuracy of our estimates is about 10%. The sum rule 

^ ""Mg/fp- + /b/A/s and -F#« = Fl^MB/fp-. The 
for A 



results are recalculated from [115 according to the relation 



results from [116[ are recalculated according to —Fa.v = 
Gv is found in the range — k < 2.0. 



1 116 



0.5 GeV. According to our estimates, k related to 





Disp Approach 


SR 


1151 


pQCD 1 




G"' 


0.120 
0.092 
0.002 n 


0.110 
0.091 


>0.09 
>0.09 


rA 

G"" 


0.055 
0.043 
0.0055 K 


0.037 
0.046 
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Using the obtained form factors and the decay constants fs = 180 MeV, fp = 210 MeV, /" = 210 MeV, we arrive 
at the results for the penguin and weak annihilation invariant amplitudes shown in Table |33|. 

Our results for FP*^"*; ^nd Fpy' amplitudes agree with the sum rule results [ 
that the anatomy of F'^^ in our analysis is different. Namely, we have found Fa considerably larger than the sum 
rule result. But after including the contact term ~ fs which we have come to Fpy^ close to the sum rule estimate 
111. 

If the form factor Gv is not taken into account, our estimate for F^^^ agrees with the sum rules. Notice, however, 
that the contribution of Gy is enhanced for the decay compared to the B~ decay. Therefore, if the value 
K ~ —1 -. — 2 is confirmed by the future analysis, the F^q^ for the p^'-f decay might appear sizeably increased. 

For a definite conclusion a more accurate study of the form factor Gy is necessary. 



Table 33. Results for the invariant amplitudes and their ratios for the B ^ P7 and B" p^'y decays. In the Standard 
Model |^u/Ct| — 0.4. According to our estimates — k < 2. 



Quantity 


B~ p"7 


-> p<'7 


i^P'="s^ in MeV 

F^^, in MeV [this work] 

F^^, in MeV [this work] 


20 6 
— 5.6ai5„ 
-6.0(1 + 0.06 k) ai C„ 


14 6 
1.9a2^u 
2.5(1 - 0.3 k) aaCi. 


^WAy^pcng ^^j^j^ ^^^.y 


-0.28aiCn/6 
-0.33(1 + 0.06 K)aiCn/6 


0.14a2 iu/it 
0.18(1 -0.33 k) aa C«/6 


^pWA/j,pcng 1 125,115 
^WA/j,pcng j 125,115 


— 0.3ai iu/it 
-0.3 ai iu/it 


0.15 aa iu/it 
0.15 aa iu/it 



115]. We would like to notice, however. 



The calculated form factors and the ratios of the weak-annihilation to the penguin amplitudes is one of the necessary 
ingredients for the calculation of the branching ratios of the B ^ p'y decays, Isospin and CP Asymmetries. In addition 
to the above ratios, these quantities contain the phase induced by the strong interactions and th e CP -violating phase 
of the CKM matrix (see [ |125[ for details). The corresponding analysis was done recently in |125| using the value 
jT^WA/^pong _ _Q 3 ^^j^^ ^WA/^pong _ q Qg f^^. ^- ^ ^-^ _^ ^0^ dccays, respectively. 



E. Discussion 



We have analysed the weak annihilation for the radiative decay i? — s- p7 in the factorization approximation. 

1. Making use of the dispersion approach, we have calculated the form factors Fa and Fy describing the radiative 
semileptonic weak transition B — > "flv. They are given by the diagrams with photon emitted from the B meson loop. 
We have performed the 1 / rrib expansion of the spectral representations for these form factors and demonstrated them 
to exhibit a behaviour in agreement with the large-energy limit of QCD. 

2. We have analysed the contribution to the weak annihilation amplitude from the diagram when the photon is 
emitted from the loop containing only light quarks. For the parity-conserving process this quantity is related to the 
divergence of the axial-vector current 

{Qu + Qd)^^^pMp|Ml. (5.55) 

with K staying finite in the chiral limit m ^ and estimated to be — k < 2. The anomalous behavior of this 
amplitude has the similar origin as the anomalous behaviour of the matrix element (7(gi)7((72)|9;^^/j|0). Let us stress 
this important result: the axial- vector current is not conserved in the chiral limit opposite to the common belief. 

3. We have also included contact terms which were missed in some of the previous analyses. Numerical estimates for 
the weak annihilation contribution to the B ^ p'y amplitude are given in Table B2. 
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VI. NONFACTORIZABLE EFFECTS IN THE B° - B° MIXING 



In this Chapter we analyse corrections to factorization for the B — B mixing amplitude assuming that non- 
factorizable soft gluon exchanges can be described in terms of the local gluon condensate. The Chapter is based 
on Ref. [g. 

Within this approximation the (asGG) -correction to the mixing amplitude can be expressed through specific B- 
meson transition form factors at zero momentum transfer. First, the correction is demonstarted to be strictly negative 
independent of the particular values of the form factors. 

As a next step, we study these form factors making use of the relativistic dispersion approach. 

The study of the oscillations in the system of neutral B mesons provides important information on th e pattern of 
the CP violation in the Standard model and its cxtentions (a detailed discussion can be found in [|l26| ). The main 
uncertainty in the theoretical description of the process arises from the nonperturbative long-distance contributions 
to the mixing amplitude. In the B — B mixing there are two kinds of such contributions: first, effects related to 
the presence of the i?-mesons in the initial and final states, and, second, corrections to the weak 4-quark amplitude 
due to the soft gluon exchanges between quarks of the initial B and the final B meson. By neglecting the second 
contribution the amplitude factorizes. In this case all nonperturbative effects connected with the meson formation 
are reduced to only one quantity - the leptonic decay constant Jb- 

Soft-gluon corrections to the weak 4-quark amplitude give rise to non-factorizable contributions, such that the 
influence of the _B-meson structure is no longer described by the decay constant only. The understanding of the actual 
size of the non-factorizable effects and their adequate description in the heavy quark systems is an important and 
challenging task. 

The theoretical analysis oi B ~ B mixing was first performed within the QCD sum rules 1 127- 129 1 and later using 
lattice QCD, see [130,131] and refs therein. The results from the sum rules are in all cases well compatible with 
factorization (vacuum saturation), within rather large errors. Recent lattice analyses definitely reported negative 
non-factorizable effects, around 5-10% for the central values at the scale ^ ~ rub. Still, the errors of the calculations 
have a similar size. Results from these approaches are shown in Table B4. 



Table 34. Comparison of results from various approaches 



Ref. 




Bb^ (mb) 




BsJBBa 


SR [ 
SR [ 


127,|12^ 
12C| 


1.0 ±0.15 
0.95 ±0.1 








Lat 
Lat 


m] 


0.92(4)+^ 
0.93(8) 


0.91(2):^g 
0.92(6) 


1.41(6):» 
1.38(11) 


0.98(3) 
0.98(5) 


This work 


0.94 ± 0.035 


0.95 ± 0.03 


1.4 ±0.05 


1.01 ±0.01 



We analyse corrections to factorization in the B — B mixing due to soft gluon exchanges in the leading Qg-o rder , 
assuming the main effect of such exchanges to be described by the local gluon condensate introduced in Ref. [132|. 
Then the parameter ABb which describes correction to factorization, is determined by the meson transition form 
factors at zero momentum transfer. 

We demonstrate that in the local gluon condensate approximation the value of ABb turns out to be negative^ 
independently of the particular values of these form factors. 

As a next step we calculate the relevant form factors making use of the dispersion approach. We obtain spectral 
representations for the form factors in terms of the i?-meson wave function and the propagator of the color-octet 
6g-system. 

Parameters of the model, such as the quark masses and the wave function have ben fixed, but for numerical estimates 
we also need the propagator of the color-octet g6-system in the nonperturbative region, so we discuss the possible 
form of this propagator. The lack of reliable information about its details in the nonperturbative region leads to the 
main uncertainties in our numerical estimates. 

Th e Chapter is organised as follows: The next Section presents basic formulas for the B ~ B mixing. In Section 
VI B we calculate the non-factorizable contribution to the mixing amplitude induced by the local gluon condensate 
and show that the correction is neg ative. In Section VI C we calculate the i?-meson transition form factors within the 
dispersion approach. Section VI D gives numerical estimates. 
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A. Effective Hamiltonian and the structure of the amplitude 



The effective Ham iltonian which summarizes the perturbative QCD corrections for the B'^ — B*^ mixing has the 
following form 133 



jjAB=2 ^ ^F^'^V (y*^Vuf Cifl) dO^b ■ dO^b + h.c. 



V2 



(6.1) 



where Gf is the Fermi constant and Co- = 7cr(l ^ 75)- 

The parameter /i stands for the renormalization scale which separates the hard region taken into account pcrtur- 
batively and the soft region. The Wilson coeffic ient C(/i) includes perturbative corrections above the scale fi. The 
explicit expression for C(/i) can be taken from |133|| . The soft contributions are contained in the i?-meson matrix 
elements of the operators in the effective Hamiltonian 

A = {B°\dO„b ■ dOab\B°). (6.2) 

The diagrammatic representation for the amplitude A of Eq. (^^) is shown in Fig. p6|. 



A=2 




Fig. 26. Different representations for the B° — mixing amplitude. The enclosed quark lines show the order of contraction 
of the spinorial indices. The circles stand for the initial (final) B (B) meson. The second line is obtained from the first line by 
performing the combined color and spinorial Fierz rearrangements. 



Using the language of the hadronic intermediate states one finds that the contribution of the hadronic vacuum to 
this amplitude gives 



(i?°|i/f,f=^|B°) 



3 {^tbvtd) G(^)/b, 



with fs the leptonic decay constant of the B meson defined according to the relation 

(0|d7^756|B(p)) ^^fBP^., {B{p)\d^^j5b\0) = -ifsPt.. 
It is convenient to parametrize the full amplitude as follows 



{B°\H^,^r^\B") 



8 Gl Ml 



V2 



{v:,VtdYC{^l)fi 



Bp 



(6.3) 



(6.4) 



(6.5) 



such that the quantity Bb — 1 = ABb measures contributions of the non- vacuum intermediate hadronic states. 

Using the language of quarks and gluons, the corrections to the amplitude of Fig. 26 due to soft gluon exchanges 
are obtained by inserting the (soft) gluons between the quark lines in Fig. Soft gluon exchanges between the 
quarks of the same loop lead to the as-corrections either to the meson vertices or to the quark propagators. They 
only contribute to the leptonic decay constant but do not lead to non-factorizable effects. 

The non-factorizable effects originate from the soft gluon exchanges between quarks of different loops. To lowest 
cts-order these effects are described by the 4 diagrams shown in Fig. 




Fig. 27. Soft-gluon exchanges leading to the non-factorizable effects 
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We assume that the main effect of the soft gluon exchange can be described by the local gluon condensate. In the 
next section we demonstrate that in this approximation the (asGG')-correction to the factorization is negative. 



B. AB in terms of the local gluon condensate 



Let us consider the exchange of a soft gluon between different quark loops assuming the dominance of the local 
gluon condensate |132]. In this case a typical graph of Fig. describing the soft-gluon contribution is reduced to the 
product of the three-point diagrams as shown in Fig. Gq. 



A1 



V 




4% 



Fig. 28. Diagrams for TfJ-'^\p,q) (a) and T!^T\p,q) (b). 




The corresponding amplitude has the form 



dx dx' 
(2^(2^ 



-^i^+^i'^\Ai(x)Ai,{x')), 



where Taa^^\p,q) with the superscript / stands for the amplitude of Fig. |28| a with the B in the initial state, and 

T^^y (p, q') with the superscript r for the amplitude of Fig. |28|b with the B in the final state. 

It is convenient to use the fixed-point gauge Xf^A'' (x) = [ |134|| . In this case the gluon potential reads 



(6.6) 



Aiix)^-xpGi0 + ..., 



(6.7) 



where dots stand for terms with higher order derivatives. The average over the hadronic vacuum is performed 
according to the relation 



5^ 
"96" 



(6.8) 



with 6"" = - 1. Here (GG) is the positive-valued gluon condensate (^GG) = 0.012 GeV* [p2 |. 
The amplitude of Fig. then takes the form 



d 



where we have denoted 



9=" = o /TT^cra.pyP)- 



(6.9) 



(6.10) 



Notice that only the parts of the amplitudes T'^'^''^ antisymmetric in indices a and /3 give a nonvanishing contribution. 
These antisymmetric parts have the following general Lorentz structure: 



"0^,13 

rpA{r,q) ^ rpA{l,q) 
<7a,0 cra,f3 ' 

rpA{r,q) ^ rpA{l,q) 
o-a,/3 (Ja.JS ' 



fl _ -CQ f 9 _ _ 

Jl ~ Jl J J2 ^ J2 



(6.11) 
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The superscript q{q) in these expressions denotes the flavour of the quark (antiquark) to which the soft e xtern al gluon 
is attached in diagrams of Fig. The quantities /f '| are real constants. The last three lines in Eq. ( 6.11 ) are the 
consequence of the C-invariance of the strong interaction 5- matrix. As we discuss in the next section, the constants 
fl'2 can be represented as specific B-meson transition form factors at zero momentum transfer. 

We have to collect now contributions of the four diagrams in Fig. For instance, the contribution to the amplitude 
of the subprocess in which the soft gluons are attached to the b quark in the left loop and the b quark in the right 
loop reads 



'^al-f ^P)tS' (f""' 5/3/3' 



gall' 9a' 13) 



(6.12) 



Similar expressions are easily obtained for other diagrams using the relations ( 6.11 ). Taking into account the rel- 
evant symmetry factors as shown in Fig. we finally arrive at the following expressions for the factorizable and 
nonfactorizable contributions to the mixing amplitude: 



(6.13) 



with the color factor Cp ~ ■ Finally, for A_Bb we find the expression 



Pb 



(6.14) 



The quantities /i and /2 depend on the renormalization scale /z. 

Obviously, the correction to the factorizable amplitude due to the local gluon condensate is negative. Let us point 
out that the relation (6.14) is the direct consequence of the only assumption of the locality of the gluon condensate, 
dominating the nonfactorizable correction. It is completely independent of any details of the i?-meson structure. 



C. Form factors in the dispersion approach 

We now proceed to the calcuation of the form factors /i and /2. We start with the amplitude Tal°'^{p,q)- It is 
diaginal in color indices, so we write Tal°'^{p,q) — -^^Ta2ip,q)- The triangle diagram for Talip,q) is shown in Fig. 

. Throughout this section we shall omit the superscript The quark structure of the _B-meson is described by the 
vertex igi 7552^(5) / \/ Nc- One should set mi = mt, m2 = rud for the calculation of and mi = m^, TO2 = mi, for 
fid). 




Fig. 29. Diagrams for Taa{p,q) (a) and Taa,fj{p) = -^Taa{p,q)\q=o (b). 



We need to calculate T„a,i3{p) — -^T^aip, q)\q=o ( ^.ICj ). For the T„a.{p, q) shown in Fig. |29|a, its derivative T„a,i3{p) 
is given by the diagram in Fig. E9b. The trace corresponding to the Feynman diagram of Fig. E3b reads: 



Saa;f3 — Sp 



ilb{'fn2 - k2)Oc,{mi + ki)"fa{mi + ki)j/3{mi + ki) 



(6.15) 



where = 7(t(1 — 75)- Since this expression is further multiplied by Gaf3, only its antisymmetric part in a and f3 is 
necessary. The latter reads 
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Saa;l3 = {^1 - kl)Sp (^75(^2 " ^2)7^(1 " l5)lp{mi - ki)j 

= {ml - kl) [4e/3QCTi/(TO2fci + mifc2) + •iimikl^igaagfiu - gpaQav)] 
The corresponding expression for the Feynnian amphtude takes the form 

1 /■....„,.. , -1 



(6.16) 



(27r)4i 



dkidk2S{p — fci — 



{mj-kf)^ml-ki) 



We define and fs by the relations 

^Jdk,dk,5{p-k,-k,)^-^^' 
dkidk2S{p — kl — k2) 



(6.17) 



u2-)2,'^2 
1 



{mi - k^ 



{m\ — klY{m\ — kl 



2Y =P''fv{P ), 

■=fs{p')- 



After the fc-integration in Eq. (6.17) one recovers the structure of the amplitude from Eq. ( |6.11 ) 



with 



Taa-aiP) = -4^" [/i e^cCTi/ + i/2 {gaagiSu " ff/3a5aiv)] 



/2 = '71i(/5 - /.„). 



(6.18) 



(6.19) 



(6.20) 



Let us obtain now spectral representations for the form factors and . We first notice that the expressions ( ^.1^ ) 
correspond to the triangle Feynman diagram at zero momentum transfer q = 0. This condition means that both 
g2 = and = p'^, where p' = p — q, so we can write ff {p^) = ff {q^ = 0,p^,p'^ = p^), with i = v,s. 

It is convenient to start with the case ^ and p^ ^ p'^. Then the form factors fi can be written in the form of 
the double spectral representation: 



rF r 2 2 I2\ 
h [Q ,P ,P ) = 



ds 



ds' 



s — p^ s' 



7jA,(s,s',9^), 



(6.21) 



where the double spectral density A.; can be calculated for any of the form factors from the Feynman integral. 

As the next step, we set q ^ — , but still treat p^ and p'^ as independent variables. In this case the double spectral 
representations of the form ( 6.21 ) simplify to the following single spectral representations 



f^iP^p")- 



1 



167r2 
1 

167r2 



ds- 



1 



ds 



s — p^ 
1 



log 



s + m\ — m\ — \^/'^{s, ml, m^) 



1 



p^ s ~ p'- 



(6.22) 



The representations for the form factors in the form (6.22) are valid however only in the region of p^ and p'^ (far) 
below the threshold. The reason for that is the following: the representations (3.21) and ( |6.22 ) are based on the 
Feynman form of the quark propagators. This form is however valid only for highly virtual particles, while in the soft 



region it is strongly distorted by nonperturbative effects. In particular, the pole at k'^ 



in the propagator of a 



color object, like quark and gluon, is absent. Recently, this was confirmed in a lattice study of the gluon propagator 

Hi. 

The modification of the quark propagator in the soft region leads to the change of the spectral representation for 
the form factors in the region of p^ and p'^ near and above the b q thr eshold. 

Notice that the quantities l/(s — p^) and l/(s' — p'^) in Eq. ( |6.21 ) are the propagators of the initial and final bq 
states with virtualities s and s', and the squared masses p^ and p''^, respectively. The nonperturbative effects which 
modify the quark and gluon propagators in the soft region, modify the propagators of the bq states as well. 

We know the character of these changes in the p^.channel corresponding to the B-meson: soft interaction of quarks 
effectively replace the factor l/(s — M^) with a regular B-meson soft wa ve fun ction ^pb{s) in the integrand of ( |6.22[) . 

In the color-octet p'^-channel the Feynman propagator l/(s— p'^) in Eq (6.22) is replaced by the propagator D{s,p ) 
which involves proper modifications in the soft region. We therefore find the following representation for the form 
factors /,(p'2) = /,(g2 ^ o,p2 = Ml,p'^) 
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1 



167r2 
1 

167r2 



We know that D{s,p'^) ^ l/{s - 
we do not know details of D{s,p' 



ds (t>B{s)D{s,p' 
''2) at large values of s 



Ai/2(s, ml, ml) 



(6.23) 



and also that D(s,p'^ ) is f inite at s = p'^. However, 
in the soft region. Motivated by the discussion in |136|, we assume here that the 



p- 



nonperturbative effects can be described by the following modification of the propagator 

1 



D{s,p'- 



pi. 



Mi 



(6.24) 



where Mq is a mass parameter. In order to guarantee the absence of the pole in D{s,p''^) in the heavy quark limit, 
Mq should scale with the heavy quark mass as follows: AIq = 0{Aqcd 'mq). So it is convenient to write Mq in the 
form 



M2 



wrrid mf,. 



(6.25) 



where ~ Age d is the constituent mass of the light quark, and w is a parameter of order unity.^ 
In the next section we make use of Eqs. (|6.23|) and (3.24) to analyse /u,s(p'^) for Bd and Bg mesons. 



D. Numerical results. 



The parametrization for the _B-meson wave function and values of the quark masses have been already fixed. Recall 
that we have determined the wave function by fitting the lattice data on the weak transition form factors for large 
momentum transfers at the normalization point /i ~ 5 GeV. Therefore, the i3-meson soft wave function and the form 
factors fs,v[p''^) at this scale are determined. 

Whereas the B-meson wave functi on is known quite well, good information about the details of D{s,p'^) is lacking. 
We therefore use the simple Ansatz (6.24) for D{s,p'^) in the full range of s and p'^ and treat w as a free parameter 
of order unity. We expect that the variation of w in the interval w = 0.5 2 provides reasonable error estimates for 
the form factors and AB. Notice that the value of AB corresponding to w = 1 agrees favourably with the recent 
lattice estimates, see Table ^ 

The form factor /i''' which gives the main contribution to ABb is shown in Fig. ^ vs p'^. 

Table gives numerical results for the form factors at = 0. The relative magnitudes of these form factors can be 



easily understood taking into account their scaling propertiers in the heavy quark limit: /, 



(d) 



,-1/2 j.{b) 



-3/2 



Table 35. Form factors and ABs, with q — d for Bd, and q = s for Bs. The factors /i,2(p'^) are calculated at p'^ 
for Bd and at p'^ — Mg for Bs. The errors correspond to the interval ^^;=0.5-^2. 





Bd 


Bs 


Ab) _ Aq) 
Jv Jv 


9±3 X 10"^ 


8 ± 3 X 10"^ 




1.1 ±0.2 X 10'^ 


9 ± 3 X 10"^ 




1.4 ±0.5 X 10"^ 


1.0 ±0.4 X 10"^ 


M 
Ji 


7 ± 1.5 X 10"^ 


8±3 X 10"^ 


Ab) 


5.5 ± 1.5 X 10"^ 


5.5 ± 2.0 X 10"^ 


Ad) 

Jl 


7.5 ±2.5 X 10"^ 


7.5 ± 2.5 X 10"^ 


Ad) 


3.5 ± 1.0 X 10"^ 


3.0 ± 1.5 X 10^^ 


ABb, 


-0.06 ±0.035 


-0.05 ±0.03 



* Assuming that for the color-octet light qq system Eq ( [6.2^ ) remains valid with A/g = uim 
for vj = 1. This is close to the value of the gluon propagator D{k^ = 0) ~ 18 GeV~^ as found in [13. 



we find Dik^ = 0) ~ 15 GeV~^ 
This agreement seems to 



be reasonable since one can expect the propagator of the light qq color-octet system to have a structure in the nonperturbative 
region similar to the structure of the gluon propagator. 
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Fi s- 30. The form factor fsij)'"^) for the Bd meson vs p'^ . Sohd curves are results of the calculation via the Eqs. (6.2J) 
and (5.24) for various values of the parameter w. upper curve - w = 0.5, middle curve - w = 1.0, lower curve - w = 2. The 
dash-dotted curve in the region p'^ < (rub + rndf corresponds to w = 0. The vertical dotted line corresponds to p'^ = Mg. 



Fig. 30 shows the form factor /i'''' whi ch giv es the main contribution to ABb- 



The value of AB is obtained from Eq. ( |6.14D taking (^GG) = 0.012 GeV* from [132). Using the Wilson coefRcent 

at the scale /i — nib we obtain the renorm- invariant factors Bb^ ^^'^ (the definition are given in |131[ |) listed in 
Table |4[ 

Our results are in good agreement with the lattice results, except for a slightly different estimate of the SU(3) 
violating effects in the Bd — Bd and Bg — Bg cases. This nice agreement gives support to our Ansatz for the description 
of the nonperturbative effects in the propagator of the c olor-o ctet q q system. 

We would like to notice that taking the Ansatz (6.24) and (6.25) for the propagator of the sq system and setting 
Mg — w iris rrig, with w ~ 1 allows a parameter-free estimate of the K'^ — mixing. In this case the result is 
much more stable with respect to the particular value of w and therefore to the details of the propagator in the 
nonperturbative region. We obtained this way AS/f (1 GeV) = —0.21 ± 0.04. 



E. Discussion. 



In this Chapter we have considered corrections to factorization in the BB mixing amplitude due to soft-gluon 
exchanges, assuming that the main effect of such exchanges can be described in terms of the local gluon condensate. 

1. We have proved that within this approximation the correction of order ^^(GG) to factorization is negative. It can 
be expressed through the specific B-meson transition form factors at zero momentum transfer. 

2. We calculated these form factors within our dispersion approach. Spectral representations for the form factors 
in terms of the B-meson soft wave function and the propagator of the color-octet g6 system were obtained. The 
behaviour of this propagator in the nonperturbative region was discussed. The obtained numerical estimates for ABb 
favourably compare with the recent lattice calculations. 

Let us point out that the proposed approach can be extended to the description of more complicated problems. In 
particular, it can be applied to the analysis of the non-factorizable effects in non-leptonic B decays. 
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VII. QUARK-BINDING EFFECTS IN INCLUSIVE DECAYS OF HEAVY MESONS 



As we have discussed in Chapter II, the dispersion approach allows us to represent amplitudes of the hadron 
interactions as relativistic spectral representations in terms of the hadron wave function. In previous Chapters we 
analysed in detail form factors of meson decays and meson wave functions. 

In this Chapter we apply the dispersion approach to the analysis of quark-binding effects in inclusive semileptonic 
decays of heavy mesons. We derive spectral representations for various differential distributions, such as electron 
energy spectrum, q^- and Mx-distributions, in terms of the same meson wave function [ p7| . 

Using the quark-model parameters determined in the previous Chapters, we provide numerical estimates for various 
distributions describing inclusive B — > Xcive decays. 

Inclusive decays provide a promising possibility to determine those CKM matrix elements which describe the mixing 
of b quark. This is due to the fact that a rigorous theoretical treatment of these decays, including nonperturbative 
effects, is possible. The consideration based on the Operator Product Expansion (OPE) and the Heavy Quark (HQ) 
expansion Q allows to connect the rate of the inclusive B meson decay with the rate of the b quark decay. An 
important consequence of the OPE is the appearance of quark-binding effects in the integrated rates (both total 
and semileptonic) of heavy meson decays only in the second order of the 1 /rng expansion |5| . These second order 
corrections are expressed in terms of the two hadronic parameters, Ai and A2. The latter are the mesonic matrix 
elements of operators of dimension 5 which appear in the OPE of the product of two weak currents. The differential 



distributions can be obtained in the form of expansions in inverse powers of the heavy-quark mass mg |137-139|. 

Whereas providing quite reliable results for the integrated semileptonic decay rate, the OPE method encounters 
difficulties in calculating various differential distributions. For instance, before comparing the OPE-based results for 
the differential distributions with the true distributions a proper smearing over duality interval is necessary. 

There are several reasons for complications arising in the calculation of differential distributions in the resonance 
region near zero recoil, namely: 

1. duality- violating I/tuq effects (i.e. the difference between the true distributions and the smeared OPE results) in 
the differentia l dis tributions near zero recoil which originate from the delay in opening different hadronic channels, 
as noticed in |l40| . Although these effects are cancelled in the integrated semileptonic rate, they can considerably 
influence the kinematical distributions near the zero recoil point; 

2. the convergence of the OPE series for the differential distributions persists only in the region where the quark 
produced in the semileptonic decay is sufficiently fast. This means that the OPE cannot directly predict distributions 
in some kinematical regions, such as: 

• the photon energy spectrum dT/dE-y in the radiative B Xg^: the window in the photon energy between 



raqjl and Mq/2 turns out to be completely inaccessible within the OPE formalism [139|; 



• the lepton energy spectrum dT/dEi at large values of Eg in semileptonic or rare semileptonic decays; 

• the lepton g^-distributions in semileptonic B — > Xc (Xu) and rare B ^ Xg decays at large near zero recoil; 
in this region one encounters both the quark-binding and duality-violating effects. 

Problems related to the quark-binding effects can be solved in principle by performing proper resummation of 
the nonperturbative corrections which in practice however leads to the appearance of a priori unknown distribution 



functions [|14l|,|42 



The inclusion of the quark-binding effects in heavy meson decays was first done in |143], where an unknown 
distribution function of a heavy quark inside the heavy meson was introduced. Evidently, this distribution function 
is connected with the wave function of the heavy meson which also determines the exclusive transition form factors. 
To put this connection on a more solid basis, it is reasonable to consider the inclusive and exclusive processes within 
the same approach. 

We show in this Chapter that the dispersion approach based on the constituent quark picture can be used as an 
efficient tool for calculating differential distributions in inclusive decays of heavy mesons, covering also kinematical 
regions where OPE cannot provide a rigorous treatment. Our dispersion approach allows one to take into account 
quark-binding effects in inclusive heavy meson decays in terms of the meson soft wave function. The latter describes 
the heavy meson properties both in exclusive and inclusive processes and thus allows one to consider on the same 
ground long-distance effects in various kinds of hadronic processes. 

Quark-model calculations of inclusive distributions are essentially based on the evaluation of the box diagram (see 
Fig. |3^ later on) by introducing the heavy meson wave function in one way or another. To illustrate the basic features 
of such an approach as well as its advantages and limitations it suffices to consider the case of a nonrelativistic 
potential model with scalar currents. Inclusion of relativistic effects can be then performed. 
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Let us consider a weak transition induced by the scalar current J = c6, where both b and c are heavy. To make the 
nonrelativistic treatment consistent we assume that 

rrif,, rric ^ Sm = nib — ^ A, (7-1) 

where A is the typical scale of the quark binding effects in the heavy meson. In the nonrelativistic theory the general 
expression for the hadronic tensor 

Wiqo,q) = ilm J (B|T( J(x) J+(0))|S)e-'«^da; (7.2) 

is reduced to the form 

W{qo, q) = -lm{B\G,s{MB - q° - tO, q)\B). (7.3) 
n 

Here 

G,^iE,q)^iH,M~E)-' (7.4) 
is the full Green function corresponding to the full Hamiltonian operator of the cd system with the total momentum 

q 

H^M = m, + nu + ^^-t^ + ^ + V,,{f). (7.5) 

Thus, the hadronic tensor is the average of the full cd Green function over the ground state of the full bd Hamiltonian 

H,^\B) = Mb\B) = (mb + mrf + eB)\B). (7.6) 
In the rest frame of the i3-meson one has 

H,, = m, + m, + ^ + ^ + 14,(f) 

= TOb + nid + hf^d (7-7) 



The following relation provides a basis for performing the OPE in the nonrelativisitc potential model [ 144 
Hcd - {Mb - 9°) = -{5m - ^ - 9°) + ij^bd - 

-H^r^-^V^i + of^) (7.8) 



2 \mc rrib 

where 7 ~ Age d and we have assumed the following expansion of the potential 



Vq,^Vo + -^ + -\ + ... (7.9) 

2mQ 2mg 



Starting with (7.8) one constructs an OPE series using the amplitude of the free 6 — > c quark transition as a zero-order 
approximation (hereafter referred to as the standard OPE). By virtue of the equations of motion, 

(A^j-es) |S) =0, (7.10) 

one observes the absence of l/mq corrections to the leading order (LO) b ^ c amplitude, so that the corrections 
emerge only at the I/tiq order. A detailed analysis of the duality in the nonrelativistic potential model can be found 
in 1 144 1 . Being completely reliable for the calculation of the integrated decay rate, the choice of the free-quark decay 
as the zero-order approximation turns out to be inconvenient however for calculating differential distributions. In 
particular, the distribution in the invariant mass of the produced hadronic system, AIx, becomes very singular and 
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is represented via 5{Mx — rric) and its derivatives, such that the l/mq corrections are even more singular than the 
leading-order result. This is the price one pays for the choice of the zero-order term. 

It is clear that the free-quark transition amplitude is not the unique choice of the zero-order approximation, at least 
in quantum mechanics. For instance, another structure of the expansion can be obtained if the free cd Green function 
is used as the zero-order approximation. 

In the nonrelativisitc quantum mechanics the relation between the full G{E) and the free Gq{E) Green functions 
is well known and reads 



G-^{E)=H-E, G^^{E) = Ho~E, G'^ (E) ~ G^^ (E) = V, 



(7.11) 



or, equivalently, 



G{E) 



Go{E) 
GoiE) 



Go{E)VG{E) 
Go{E)VGo{E) 



-Go{E)VGo{E)VGo{E) 



(7.12) 



For the heavy quark decay, in most of the kinematical q^-region except for a vicinity of the zero recoil point, the 
Green function Go behaves as l/mg, and, since the matrix elements of the operator V remain finite as mg — > cx), 
t he se ries ( 7.12 ) is an expansion in powers of l/mq. Notic e th at in the nonrelativistic p otenti al model the expansion 

Inserting the expansion (7.12) into the expression for 



(7.12) is fully equivalent to the OPE serie s obtained from (7. 



the hadronic tensor W given by eq. (7.3), we come to the series shown in Fig mJ. 
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Fig. 31. Expansion of the hadronic tensor in the quark model: (a) - the box diagram which provides the LO free-quark 
contribution, (b,c) - diagrams contributing in subleading l/mg orders and containing the final state interaction V = V^^. The 
short-dashed lines represent W bosons. 



The leading contribution in the heavy quark limit is given by the box diagram of Fig. |3l| a with the free c and d 
quarks in the intermediate state. The corresponding analytical expression reads 



t^QM ^ Um{B\Gl^{MB - q\q)\B)- 



(7.13) 



This is the quantity usually taken into account in quark-model calculations. It is easy to see, that the semileptonic 
decay rate calculation based on the box diagram of Fig. |3^ only, reproduces the free quark semileptonic decay rate in 
the heavy-quar k lim it, but should contain also I/toq corrections (see also the general s truct ure of the quark-model 
results of Ref. |145|). Namely, the next-to-leading order (NLO) term in the expansion (7T2) is the diagram with a 
single insertion of the potential between the free c and d quarks (the diagram of Fig. [31|3). It has the order 1 /mg 
and precisely cancels the l/mg contribution of the quark-model box diagram, yielding in this way the absence of the 
l/mg correction in the difference between the decay rates of bound and free heavy quarks. 

In other words, the quark-model box-diagram calculation is just the first term in an alternative expansion of the 
full Green function: unlike the standard OPE series which starts with a single c-quark in the intermediate state, the 
quark-model starts with the free cd pair which is the eigenstate of the Hamiltonian 



rrid 



2mr 



2md 



(7.14) 



Hereafter we refer to the expansion of the decay rate based on the expansion ( 7.12 ) of the Green function as the 
quark-model (QM) expansion. 
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Summing up, the quark model provides an alternative Xjmq expansion with the following properties: 

1. the box diagram of Fig. ^l| a provides the LO l/mq term and reproduces the free-quark decay in the limit 
mg oo. All other terms contribute only in subleading l/mg orders; 

2. the differential distributions in any l/mg order are convergent for almost all allowed q^, except for the region 
close to zero recoil; 



3. before comparing the calculated differential distributions based on the expansion (7.12) with the true distribu- 



tions in the resonance region a proper smearing over some duality interval is required; 
4. the 1/toq correction to the LO term is nonvanishing. 

Clearly, the properties 1-3 are completely equivalent to the standard OPE, while the property 4 makes the quark- 
model expansion much less convenient than the standard OPE, at least for the calculation of the semileptonic decay 
rate. However: 



5. the expansion (7.12) turns out to be more suitable for the calculation of the differential distributions, e.g. for 
the calculation of dT/dMx' in this case the LO result (the box diagram of Fig. |3l|a) is well-defined in the whole 
kinematical region as well as the higher order corrections to it. Thus already the box diagram is appropriate for 
comparison with the experimental dT/dMx at all Mx apart from the resonance region. Beyond the resonance 
region no additional smearing of the calculated dT/dMx is required. 

In full QCD the situation is of course much more complicated. Hadrons are coherent states of infinite number of 
quarks and gluons. Nevertheless, many applications of the constituent quark model have proved the treatment of 
mesons as bound states of two constituent quarks to provide a reasonable description of their properties. From this 
viewpoint the arguments given above remain valid. Namely, the box diagram represents the main contribution to the 
hadronic tensor which reproduces the free-quark decay in the infinite quark mass limit. However the hadronic tensor 
calculated from the box diagram contains l/mg term compared with the free-quark decay tensor. This linear l/mg 
term is known to be cancelled by the l/mg contributions of higher order diagrams. In practice, however the l/mg 
term of the box diagram is not so dangerous: namely, the hadronic tensor calculated from the box diagram, as well 
as all corrections given by the other graphs, are regular in the whole kinematical region. Thus the box diagram should 
provide a reasonable description appropriate for comparison with experiment. Moreover, the box-diagram result can 
be further improved by effectively taking into account the higher order term which kills the 1 /mg correction contained 
in the box diagram. We follow this strategy in our analysis and perform a relativistic treatment of quark-binding 
effects within a constituent quark picture. 

Our consideration of the quark binding effects in inclusive semileptonic decays is based on the relativistic dispersion 
formulation of the quark model discussed in detail in the previous chapters. 

The inclusive decay rates as well as the exclusive hadron transition form factors are given by double spectral 
representations in terms of the soft meson wave functions. Recall, that the double spectral densities of these spectral 
representations are obtained from the corresponding Feynman graphs, whereas subtraction terms should be fixed 
independently. Such subtraction terms in the spectral representations for exclusive transition form factors were fixed 
by requiring the structure of the heavy quark expansion in our approach to match the structure of the heavy quark 
expansion in QCD. In this Chapter we proceed along the same lines in inclusive processes. 

• We construct the double spectral representation of the hadronic tensor within the constituent quark model 
starting with < 0. The hadronic tensor is represented in terms of the soft wave function of the B meson and 
the double spectral density of the box diagram. The hadronic tensor at g~ > is obtained by the analytical 
continuation. Then the l/mg expansion of the spectral representation of the decay rate is performed and the 
LO term is shown to reproduce the free quark decay rate. The subtraction is defined in such a way that the 
l/mg correction to the semileptonic decay rate is absent. This corresponds to effectively taking into account 
other terms beyond the box-diagram approximation which contribute in subleading 1 /mg orders. Moreover an 
account of the l/mg effects of the whole series within the box-diagram expression is possible. This is done by 
introducing a phenomenological cut in the double spectral representation of the box-diagram which affects only 
the differential distributions at large q^. This cut brings the size of the l/mg corrections in the T{B — > Xcive) 
in full agreement with the OPE result and keeps the LO and l/mg correction unchanged. The cut yields 
differential distributions which are finite also in the endpoint g^-region where the HQ expansion series is not 
properly convergent; 

• We calculate various differential distributions in terms of the B-meson soft wave function. These distributions 
are regular in the whole kinematically accessible region and, apart from the resonance region (where the exact 
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distributions are dominated by single resonances and a proper smearing over the duality intervals is necessary) , 
can be directly compared with the observable values. The main effect of quark-binding upon these distributions is 
determined unambiguously through the soft wave function, while the 1 / rrig corrections depend on the particular 
details of an account of the higher-order terms in the series ( 7.12| ). However, in practice these details are not 
essential due to the following two reasons: first they are numerically small, and second, the size of the I/jtiq 
corrections in the integrated semileptonic rate is close to the OPE result. So we expect that the size of the 
I/ttiq corrections is reasonably reproduced also in other quantities; 



• We perform numerical estimates of various differential distributions in inclusive decays in terms of the _B-meson 
wave function determined from the description of the exclusive processes. 

In the next section we present necessary form ulas fo r the free-quark decay and also the OPE prediction for the total 
integrated rate to accuracy. In Section VII B we construct the dispersion representation for the box diagram 

at < and discuss its analytical continuation to > 0. Section VII C gives the l/mg expansion of the hadronic 
tensor in the quark model, and Sections VII D and VII D present numerical results for the differential distributions. 



A. Free quark decay and OPE 

The semileptonic weak decay Q2 ^ Qi is governed by the effective Hamiltonian 

iJcff = ^1^21 • Qi(x)7,,(l - 75)Q2 • ^17^(1 - 75)^2, (7.15) 

where V21 is the corresponding matrix element of the Cabibbo-Kabayashi-Maskawa quark mixing matrix. In what 
follows we assume both leptons to be massless. We imply Q2 = b and Qi = c, and therefore use both notations 
throughout this section. 



I1 
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Fig. 32. The cut Feynman diagram corresponding to the free quark semileptonic decay (32(^2) — > Qi(fci)?i (91)^2(92). Crossed 
lines mean that the corresponding particles are on the mass shell. 



A tree-level decay rate of the free-quark semileptonic decay averaged over the polarizations of the initial quark Q2 
and summed over the polarizations of the final quark Qi is given the cut Feynman diagram of Fig. ^ and has the 
form 

r rr> r> 1 T \ l^r 12 (27r)'* f dkidqidq2 , , 

ro(Q2 Qihh) = —\V2i\ J (27i-)9 - fci - gi - q2) 

x^Sp (7,,(1 - 75)(mi + fci)7,(l - 75)(to2 + ^2)) ^^(fc?)^(fc? - ml) 
xSp (7^(1 - 75)927.(1 - 75)91) 9iq"i)S{ql)eiq°2)Hql) (7-16) 
It is convenient to introduce additional integration by inserting the unity factor 

J dqS{q - 91 - q2)dfl^Siq^ " m') = 1 (7-17) 
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in the integrand. Then we can rewrite the decay rate as follows 

fr, r, 1 T \ Gf\V2i\'^ (271)^ f dkidq ^ q ^ 2n/j/,0nc/,2 2x 
ro((92 -> Qihh) = 2 2^ J "(2^^ ^'^ ^^"^ ^ - "^i) 

x,5(fc2 - fci - g)wO,(fci, fc2)V(g). (7.18) 

Here we have introduced the following quantites: 
The leptonic tensor l^v{<l) 

I^Aq) = / dqidq25iq - qi - q2)B{ql)5{ql)e{ql)5{ql) Sp (7^(1 - 75)'?l7.(l " 75)92) 



47r 

= -j{(it^<iu- gt,,^q )■ (7.19) 

The free-quark tensor 

w%{ki,k2) = ^Sp (7^(1 - 75)(mi + ki)lu{l - l^){m2 + hfj 

= 4 {ki^k2v + k2^kii, - g^i,kik2 + ie^ivai3k2akij3) ■ (7.20) 



We also denote 



CsL{ml,ml,q^) ^ ]^wlAki,k2){q^,qu - g^iuq^) 



= 2kiq ■ 2k2q + q^ ■ 2kik2 

= {ml~mlf +q^{ml+ml)-2q'^. (7.21) 



By virtue of these formulas we come to the representation 

G||Fb,|2 1 

Switching to mf, = TO2 and rric — mi, the differential decay rate takes the familiar form 

:'2 n/; 12 



ro(g2 ^ Qi;ir2) = 3 / d^x''\^'\mlml^i^)CsL{mlml^,''), (7.22) 



dVaib^clv) Gp\Vbc\ 1 1/2. 2 2 2x^ / 2 2 2n oQ^ 
rf^^ " 96^3 {mb,m^,q )CsLim„m^,q ). (7.23) 

The integrated semileptonic decay rate is then given by the expression 

r„(6^cM= / (7.24) 



with 

/o(r) 1 - 8r + 8r^ - r"* - 12r2ln(r), r = {mc/rribf. (7.25) 

The operator product expansion in QCD predicts the absence of the l/mh corrections to the ratio of the bound to 
free quark decay rates |^,|| . With the 1 /ml accuracy the integrated rate of the inclusive semileptonic B ^ Xc decay 
is given by the expression 



A 1+3A2 gA2 (l-r)^ 
2ml ™fc Io{r) 



T{B ^ XM - ro(6 ^ civ) • ( 1 + , - Q-±^——L. ) . (7.26) 



In this expression Ai and A2 are the hadronic matrix elements of the operators of dimension 5 appearing in the 
OPE of the product of the two weak currents. The value A2 = 0.12 GeV^ is well known from the B — B* mass 
splitting, whereas the knowledge of Ai is loose and present estimates range from —0.6 GeV^ to 0. Note that in the 
nonrelativistic quark potential model one has Ai = — < fc^ >, where k is the relative momentum of the constituent 
Qq pair. Typically, the nonrelativistic quark model estimates of Ai range from —0.6 to —0.4 GeV^. 



98 



B. Inclusive meson decay in the dispersion approach 



We now proceed to the inclusive rate for the decay of a pseudoscalar meson with mass Mi containing a heavy quark 
Q2, which we refer to as Pqj, induced by the quark semileptonic transition Q2 Qihh- 

As discussed above, within the quark model the dominant contribution to the inclusive decay rate of any process 
for tuq 00 isgiven by the box diagram of Fig. ^l|a. For semileptonic inclusive decay the corresponding diagram is 
shown in Fig. ^ We start with this diagram and discuss the relevant modifications due to the other terms presented 
in Fig. |3l|b,c later. 



Pi ^ 




Fig. 33. The cut Feynman diagram corresponding to the decay Pq^ XQ^/i(gi)/2(92). Crossed lines mean that the 
corresponding particles are on the mass shell. We denote p2 — fci + fcg. Our notations follow those of Chapter 3 where 
transition form factors have been considered. 



Rec all that the quark structure of the pseudoscalar meson transition into two real quarks is described by the vertex 

(Hi 



G(5i) 
^/N: 



(7.27) 



where pi = si, /c| — rn^, fc| — m§. 

Notice that the upper leptonic block of the diagra m in Fig. |3^ is exactly the same block Z^i, as for the free quark 
decay. The leptonic tensor lf^„ is given by Eq. ( 7.f9 ). 

Within the dispersion approach the differential decay rate reads |^ 



dViPg, ^ XQ.hh) _ Gl 1^211' (2^)4 1 



2 2Mi (27r)i2 

dsi G(si) ds'i G(s;) TT Al/2(s^, g2) 



(7.28) 



2si 



lf,y{q)w^,i,{pi,p2,n)\K^^o. 



The hadronic tensor Wfj^^(pi^p2T n) represents the full double spectral density of the Feynman box diagram 
w^j^(pi,P2) 1^) of Fig. 34 a. This quantity is constructed from the double discontinuity wj^,^{pi,p2, k) 



Wf_„^{pi,P2,K) = discsi disCs; w^^{pi,p2, k), 



(7.29) 



by applying the appropriate subtraction prescription. Fig. |3^ shows the momenta in the integrand of Eq. ( [7.28 ) 
They satisfy the following relations 



"This representation is appropriate for the calculation of the g^-distribution. A more detailed formula is necessary for the 
lepton energy spectrum. It will be considered later in this Chapter. 
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(7.30) 



Pi 




Fig. 34. a: The box diagram for Wfi^(pi,p2,9 — g')', b: its double discontinuity in the variables Si 
'Wfiv{pi,P2, K — q — q'). Crossed lines correspond to particles on the mass shell. 
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1. The spacelike region 



The double discontinuity of the box diagram can be calculated in the region < and q'^ < by placing all 
intermediate particles in the loop on their mass shell as shown in Fig. To obtain the corresponding expression for 
positive values q^ = q'^ > corresponding to the real decay kinematics we shall perform the analytical continuation 



m q"^. 

At q^, q'^ < the explicit expression for wj^^ reads 

wj^^{pi,'P2, k) = J dkidk2dk2dk35{k'l ~ m'l)S{k2 ^ m2)S{k2 ^ m2)S{k'^ — m^) 
xS{pi -k2~ k3)S{p\ - k'2~ k3)5{p2 - ki - fcg) 

x(-l) Sp (ma - A;3)«75 ("^-2 + ^2)7^(1 ^ 75)(nT-i + ^i)7i>(l - 75)(™2 + ^2)) ■ (7-31) 

The function w^^(j)i,p2, n) contains several Lorentz structures multipled by the invariant amplitudes. The latter 
depend on the six independent scalar variables 

si,Si,S2,g ,q ,K. (7.32) 

We are interested in the case q^ — q'^ and — > 0. 

Let us consider the trace in the integrand of w^^. Hereafter we neglect terms proportional to k. Since all particles 
are on the mass shell, we find 

- (?7i2 + /c2)«75 ("I3 - ^3)^75 ("I2 + ^2) = (2A;2^3 + 2TO2m3)(TO2 + k2). (7.33) 

Therefore, 

-Sp (^Z75 (m3 - k3)il5 {^2 + ^2)7^(1 ~ 75)(wi + fci)7i.(l - 75)(w2 + k2)^ 

= 2(si - (m2 - m3)2)u.o (fci, fc2). (7.34) 



From the free quark decay we know that 



w%{k,,k2)l^M) = ^CsLimlmlq^). (7.35) 
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Clearly, the contribution of the traces in the spectral representation for dT / dq^ turns out to be independent of the 
internal integration variables. Therefore, we only have to consider the scalar part of defined according to the 
relation 

AD{si,s\,plq'^,q'^,K'^) = J dkidk2dk'2dk36{kl - m:i)5{kl - ml)6{k'i ~ ml)5{kl ~ ml) 

x<5(pi -k2- k3)S{p[ -k'2- k^)5{P2 - fci - fca). (7.36) 
Explicit calculations give for Ad the following expression 



AD{si,S2,q) 



(7.37) 



Notice that the argument of the 6 function in eq ( 7.37 ) is just the same as for the spectral density of the triangle 
graph f or th e form factor given by Eq. (2.79). 

Eq. (7.36) provides the double discontinuity but, as usual, we need to define in addition a subtracion prescription. 
The form of this prescription cannot be determined within the dispersion approach and should be fixed from some 
other arguments. We shall determine the subtraction prescription by matching to the known QCD result and requiring 
the absence of the I/ttiq corrections in the ratio of the bound and free quark decay rates. As we have discussed, in 
the quark model the absence of the l/mq corrections is provided by the other diagrams shown in Fig. ^l|b,c which 
cancel the 1 /mq term in the box diagram. Therefore requiring the absence of the 1 /mq corrections corresponds to 
an effective account for these terms. 

Our subtraction prescription explicitly reads 



Ml 

Aisi,S2,q^) = —=ADisi,S2,q'^). 



(7.38) 



As we shall demonstrate in the next section, this prescription indeed provides the required absence of the l/mq 
corrections in agreement with the OPE. 

For the analysis of the limit mq ^ 00 it is convenient to have the 'external' integration over si and the 'internal' 
integration over S2 in the spectral representation for the decay rate. The solution of the 9 function then takes the 
form 



(m2 + ms)^ < Si < 00 
s^(si,9^) < S2 < s+(si,g2) 



where the limits are obtained by setting 7] — ±1 in the relation 



S2(si,g^) = (mi +7713)^ + — (si - {m2 +m^Y) + — (cj - l)(si - 



mi 



m2 



m2 



ml) 



I '^1 X 1/2/ 2 2\ / 9~ 

+ 7/ — A' (si, m,, m^) V 
m2 



1. 



The quark recoil uj is defined according to the relation 

q^ = (m2 — mi)^ — 2mim2(cj — 1). 



(7.39) 



(7.40) 



Isolating the free-quark decay amplitude we come to the following dispersion representation for the differential inclusive 
semileptonic decay rate at < 



where 



dT 
d^ 



Koiq^) / dsiip^is) 



si — (m2 — m^)^ m2 
87r2si 



Koiq' 



dsy 



\"\s,,S2.q^) 
\^/'^{m\,m\,q'^y 



^2 



A^/^(m2, m^, (7^) dq'^ 



(7.41) 



(7.42) 



Note that in Eq. (7.41) the free-quark differential rate dTg/dq^ factorizes out, so that the differential rate for a bound 
quark is a product of the free-quark differential rate and a bound state factor, as already noted in |145]. As in the 
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previo us chapters, we use the notation 1^9(5) = G(s)/(s — Mf); the normahzation condition for (p{s) is given by Eq. 

(Hi)- 



It is convenient to rear range Eq. (7.41) by isolating under the integral the structure similar to the structure of the 
normalization condition (3.62) 



dr _ dVo 
dq^ dq^ 



dsiip (si) [si - (m2 - ma) J - - ^ 



87r2si 



-r{si,q^), 



2, _ , drj Ai/2(si,S2,g^) 



2 Ai/2(TO2,m2,92)' 



(7.43) 



and S2 is related to 77 by Eq. ( 7.39 ). As we shall see later, 

r(si,g2)^l 



(7.44) 



in the heavy-quark limit, such that thanks to the normalization condition of the soft wave function Eq. (3.62) one 
gets 



dr dTo 
dq^ dq^ 



(7.45) 



IQ 00. 



2. The timelike region and the anomalous contribution 



To obtain the spectral representation at q^ > we perform the analytical continuation in q^. This procedure is 
done along the same lines as in the case of the transition form factor which has been discuss ed in detail in Chapter 
II. As a result of this procedure in addition to the normal part which is just the expression ( 7.43 ) taken at q^ > 0, 
the anomalous part emerges due to the non-Landau type singularities of the Feynman graph. The location of the 
singularities in the complex si and S2 planes corresponding to the situation of the 'external' si integration is shown 
in Fig. Bq. 



(m2+m3) ^1^1 



S?(q2) 



(m 2+1113) 



S?(q2) 



(mi+mj) 



(m 1+1113) 



Fig. 35. The integration region for m2 > mi and > for the order of integration J dsi j ds2'- a. si < s?, b. si > s\. 
Solid lines - cuts on the physical sheet, dashed lines - cuts on the second sheet. 



The definition of is given by Eq. (7.39) and 



^*f(sl,<z') = (^/i^+v^)^ 
s^(sl,<z') = (^/i^-v^)^ 



q^ + m\ — m\ 



\ 



q^ + m\ — m\ 



(7.46) 
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As one can see, the anomalous cut appears on the physical S2 sheet between sj(si,(7^) and 32(31, q^) in the region 
< (m2 — mi)^ for the values Si > S]'(g^). 

Respectively, the spectral density r(si, q^) acquires the anomalous piece and takes at > the following form 



2 m2 m. 
r{si,q ) ^ 



ds2Ai/2(si, S2, g") + 29iq^)9isi > s°) / ds2X'^H3i, 32,q^) 

S2 



(7.47) 



The q^-behavior of the anomalous term is determined by the lower limit of the si integration, s5(q^). Namely, its 
contribution to the semileptonic rate reads 

1 Wp'i^om A3 

-— — 3 R~{cu), (7.48) 

where 

00 

i?'^"°™(tj) J dkk^w^ik). (7.49) 

mg — 1 

Here k = X^^'^{si, m|, m|) jl^fsl is the r elativ e momentum of the quarks inside the i3-meson. The wave function w(fc) 
is normalized according to the relation ( ^.94| ) leading to 

i?™°'"(l) 1. (7.50) 

Since the soft wave function is steeply falling beyond the confinement region where k'^ ^ A^, the anomalous con- 
tribution becomes inessential already at a; — 1 ~ A/m;,. Only in the endpoint region w — 1 ^ the anomalous 
contribution to the differential distribution becomes strong and diverges like 

1 dr"'""" A^ 1 , ^ 

i aijj rriQ ^uj — 1 

The contribution of the anomalous term to the integrated semileptonic rate of the order 

panom/p ~ / m% (7.52) 

comes from the endpoint region, whereas the rest of the decay phase space provides only the relative A'^/mg anomalous 
contribution to the semileptonic decay rate. 

A 2 

Therefore, the anomalous contribution is negligible at all u except for the endpoint region ui — 1 < which IS m 

fact a very narrow region near zero recoil. As we have discussed, the HQ expansion for the differential distributions 
is anyway ill-defined in this k inematical region. Contributions of the same order of magnitude come also from other 



terms in the expansion (7.12), and keeping this anomalous contribution is beyond the accuracy of our considerations. 



Thus we shall systematically omit the anomalous contribution in numerical calculations. 
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C. The l/mg expansion of the semileptonic decay rate and dV/dq^ 



In this section we perform the l/mg expansion of the meson inclusive decay rate. We show that: 
a. in the leading l/mg order the heavy meson inclusive decay rate is equal to the free quark decay rate; 



b. our subtraction prescription (7.38) leads to the differential distribution dV/dq^ given by eq (7.43) which satisfies 
the relation 



dT{B X^lDi)/dq^ 
dTo{b^ cePi,)/dq^ 



I + Oil/ml) 



(7.53) 



in most of the q^ region except for a close vicinity of zero recoil point. This property guarantees the absence of 
the l/mg corrections in the ratio of the integrated rates, i.e. 



l + 0(l/m2); 



(7.54) 



c. the size of the l/?7ig corrections can be tuned such that they become numerically close to the OPE prediction. 
This is done by introducing the cut in the spectral representation of the decay rate of the B meson. This 
cut affects only the differential distribution dT{B — > Xcivt) / dq^ at large q^ near zero recoil, i.e. in the region 
Lo < 1 + 0(l/mg). 

An important feature of the whole approach is that already the zero order expression provides a realistic Mx- 
distribution. As we shall see later, the modifications b) and c) while affecting the total rate and the g^-distributions 
at large q^, only moderately affect the Mx-distribution, so that the latter is completely determined by the soft 
i?-meson wave function. 

Let us consider the l/mg expansion of the spectral density rjj. The normal part of r£)(si, q'^) has a simple form 



r_D(si,9^) 



ma dr] X^/^{si,S2,q^) 
Ml i_i 2 Xy^mlmlq^)- 



(7.55) 



This representation is a convenient starting point for performing the l/mg expansion. Let us point out that although 
the integration in rj can be easily performed, it is more convenient to work out the HQ expansion before the integration. 

Assuming that m2 is large and that the meson wave function is localized in the region zi ~ A we obtain the 
following expression for A(si, S2, g^) valid at all q^ 



A(si, S2, q ) m^ll-i 

V m2 J 



\{l,f,f') + ^x(^i)v/^i(^i+2m3)(l + e - f2)Ai/2(l, f, f2) 
m2 



zi(zi + 2m3) 



[l + f- + -^ziizi + 2m3)A(l, f, f2) 



(7.56) 



where q^ = q^ /rn^, f = ini/m2 and xi^i) — 1 ^ (^i + '7^3)/2m2. In the limit m2 ^ oo we can expand the A(si, S2, q^) 
in powers of l/m2. Notice however that an actual expansion parameter is not l/m2 but rather 



^Jzl{zl + 2m3) 
m2X}/'^{l, q^, r^) ' 



(7.57) 



and the averaging over the B meson state implies z\{z\ + 2m3) ~ A. Hence the region where the expansion is lastly 
converging is 



This relation can be written as 



m2Al/2(l,g2 f2)^ A. 

I^il = \^l'^(ml,m\,q^)l2m2 > A, 



(7.58) 
(7.59) 
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which means that in the rest frame of the b quark the daughter quark has a 3-momentum much bigger than A. 
The final expression reads 



/ 2n "12 / Zi+m3 

rD{si,q 77- 1 + 

Ml \ 1712 



^ ^ zi{zi + 21113) I ^ 



2m?. \ 3A(l,g^f2) 



(7.60) 



The l/mg term in the ratio of the bound to free quark distributions is generated by the (to2 + zi + m^j/Mi term in 

As we have discussed this linear 1 /mn term contained in the box diagram cancels against the 1 /mg terms coming 
from other terms in the expansion ( [7.12| ). Thus the main contribution of these terms can be taken into account by 
performing the subtraction in the spectral representation for the box diagram which kills the l/mg term. This is 
implemented by defining the subtraction prescription as follows: 

r{suq^)^^rD{si,q^). (7.61) 



After performing the subtraction and taking into account that in the heavy quark limit z{z + 27713) = come to 

the following relation 



This expansion is valid in the region of q^ such that X{\,(f'^r^) ~ 0(1), i.e. in most of the q'^ phase space except for 
the region near ze ro rec oil where A(l, (f' , f^) ~ 0. 



The expression (|7.62| ) has the following features: 

1. In the leading l/mq order the ratio R{q^) is equal to one. Thus the decay rate of the free and the bound quark 
coincide in the heavy quark limit for all q^. Beyond the leading order, the differential distribution coincide also 
within the l/mq accuracy in most of the q^ phase space, except for the region near zero recoil. This guarantees 
the absence of the l/ruQ corrections in the ratio of the integrated rates as well. Thus, our description is in full 
agreement with the OPE results within the l/mq order; 

2. Since the box diagram represents only a part of the l/rriq corrections, we cannot expect the box diagram alone 
to reproduced correctly the I/toq term in the ratio of the integrated rates F/Fo. In fact, the sign of the I /ttiq 
correction in eq. (7.62) turns out to be opposite to the OPE result (cf., e.g., with the results of refs. [138,137 
at q^ — Q). Moreover, the l/jTig term generated by the box diagram is expressed merely in terms of {B\k'^\B) , 
whereas the I/ttiq corrections of the OPE series include also {B\Vi\B), where Vi is the l/mq term appearing 



in the expansion (7.9) of the effective potential (e.g., the chromomagnetic operator in QCD). 



We argue however that it is possible to further modify the spectral representation of the box diagram to bring the 
size of the I/ttiq term developed by this modified representation in agreement with the OPE result. This procedure 
corresponds to phcnomenologically taking into account the contribution of other 1/?77q terms of the expansion ( [7.12 ). 



Omitting the anomalous contribution for the reasons explained above, the differential decay rate takes the form 

^=Ko{q^) r dsi^^is y^-^^l-^^^\ ^f^isiM,ml)^ f ^AV2(.„.2,.^). (7.63) 

where S2 depends on rj through Eq. ( [7.391 ). 

Now, our goal is to modify the differential distribution of Eq. ( 7.63| ) in such a way that the leading order result 



and the l/mq correction in the integrated rate remain intact whereas the 1/mQ term numerically reproduces the 
OPE-based estimate. 

Obviously, we may allow a strong deformation of the differential ^^-distribution at large q^ near zero recoil, where 
the heavy-quark expansion is anyway ill-defined; we may also require dT/dq^{q^ — 0) to exactly reproduce the OPE 
result within the l/mg accuracy. 

Most easily this program may be implemented through the following two steps: 

First, by introducing the factor F{si) ~ 1/(1 + k'^/m'q) which sets the l/mq term in the differential rate at q'^ = 0; 
Second, by cutting the S2-iiitegration in ( 7.63| ) at some upper limit s'!^°'^{q'^) which tunes the size of the l/mq 
effects in the integrated rate. 
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The resulting differential differential distribution takes the form 



I 2^ f , 2, N^l - (™2 " TO3)\l/2/ 2 2n 

— =Ko{q) J dsiifi {si) X ' {si,m^,m^) 



X^n^i) r ^X'/'{s,,s„q^)9{s,<sr^{q')). (7.64) 



In order not to affect the integrated rate in the leading and the subleading l/mg orders, the function S2"°^(q^) should 
satisfy certain properties. Let us look at them more closely. 

Recall, that the soft wave function <^(si) is localized in the region 

Si < ST"^ ~ (m2 +7713+ if 

where 7 is a constant of order A which docs not scale with 7nQ. Let us determine q^ through the equation 

4(sr%9o)=sr^(go) (7.65) 



where is the maximal value of S2 corresponding to 77 = 1 in ( [7.39 ). Furthermore, assume that s™°^(q^) decreases 



with q , and take into account that (■^ij (? ) is a monotonous rising function of both si and q . Then, at q < q^ for 
all Si < s™"^ one finds the relation s^(s 1,(7^) < S2^°-^{q'^), and thus the ^^-distribution does not feel the presence of 
the cut at all. For q^ > q^ the cut becomes really effective and strongly influences the g^-distribution. In order these 
changes in the cut g^-differential distribution not to change the integrated rate in the LO and l/mg order, we need 
the (7q to be not far from zero recoil such that the corresponding luq = 1 + 0{1/771q). Choosing the cut in the form 



s^°'^{q'^) — TOi + 7773 + a (^77i2 — rrii — e — j (7.66) 

where e ~ A and a is a rising function of mg, satisfies these requirements. 

The parameter e accounts for a mismatch between the quark and the hadron threshold, and the form of a{mQ) can 
be found from fitting the size of the I/tuq corrections in the integrated rate to the OPE prediction. Notice also that 
the q^ distributions obtained through the cut expression are even more realistic than those obtained from the uncut 
spectral representations. 

We could of course choose a more sophisticated parameterization of s™"^((7^) to reproduce a correct q^-behaviour of 
dT /dq^ near zero recoil point. For instance, taking into account that the lightest final meson is pseudoscalar, we can 
write ^s™aa:(g2-) _ mi+77i3 + ap{777Q){MpQ —Mp^, — v^)'^^^ yielding the correct behavior near — Mp^ — Mp^, 
where only one P— wave decay channel Pq — s- PqiIDi is opened. In addition, in the heavy quark limit the 5— wave 
transition Pq Yg'tDg requiring another functional dependence y/s^^^jl^ ~ ttli + 7713 + av{7rLQ){MpQ — Mp^, — 
•\/g^)^/^ is opened at = ^'^Pq ^ ^'^Vq, with only small delay in q^ of order A^. So the effects of opening this 
channel are even more important and should be also taken into account. However in the region of large q^ with 
few opened channels the inclusive consideration is anyway not working properly, and taking into account such subtle 
effects is beyond the accuracy of the method. So in numerical calculations we proceed with the phenomenological cut 



provided by eq (7.66) 



The parameters of the model, such as the quark masses and the wave functions were already fixed in Chapter III 
and we now use them for inclusive decays. 

We choose the parameters of the cut according to the criteria above in the following form 

e = (™Q - '7^Q') - {Mq - Mq,) 
a= 1.82 + 0.029TOQ, 

where toq' is the mass of the parent heavy quark and Mq> is the final meson lowest mass. Effectively, this means 
that = (Mq - MqO'. 

A remark is in order here: To compare the calcuatcd integrated rate F as a function of mg with the OPE result 
within the l/mg order, it is reasonable to adopt the expansions of the hadron masses Mq and Mqi also up to 
the second order in l/mg. Then one finds e = -i(Ai + 3A2) • (l/mg- - l/mg). Setting Ai = -0.44 GeV^ and 
A2 = 0.12 GeV^, yields e — O.lO/rng at 771Q1 /niQ — mc/mb ~ 0.28. We point out that in calculations for the real B 
decays we use experimental values of hadronic masses (involving all orders in l/mg). 
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The corresponding integ rated semileptonic decay rate T is plotted in Fig. ^ as a function of 1 /mg and compared 
with the OPE predictions ( |7.26| ). 




0.25 



1 / (GeV 



Fig. 36. The ratio of the integrated rates of the bound-to- free quark semileptonic decay Q — > Q'lvi vs the inverse heavy-quark 
mass 1/mQ at fixed value of niQi /niQ — nic/mb = 0.28. Dotted line - the rate calculated from the initial spectral representation 
of the box diagram which contains l/mg correction, dashed - with the proper subtraction killing the 1/mQ term but without 
tuning the size of the I/tjiq effects. Solid - final result which also includes the cut bringi ng th e size of the l/mg effects in 



agreement with the OPE prediction. Upper and lower dot-dashed lines are the OPE results (7.26) corresponding to Ai = and 
Ai = -0.6 GeV^ (with A2 = 0.f2 GeV^), respectively. 



One can clearly see that the calculations based on both the spectral density td of Eq. (7.55) and the spectral 



density r of Eq. (7.61) which includes a subtraction predict a larger rate for a bound heavy quark compared to the 



free one, that contradicts OPE. 



The introduction of the cut s^°'^{q ) (7.66) brings our dispersion approach results in perfect agreement with the 



standard OPE framework for the whole range of the considered values of mg . 
Figure 37 shows the influence of the cut upon the g^-distribution for the B 



Xctvi decay. 
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Fig. 37. Distribution dT/dq^ in _B — > Xclvi decays vs the squared four- momentum transfer q^: (a) ratio of the bound-to- free 
quark decay R(q^) = {dT / dq^) / (dVo / dq^) , notation of lines same as in Fi g. Is^ . (b) Differential distribution in the bound (solid) 
and free (dot-dashed) semileptonic quark decay. Parameters of the cut ( 7.66 ) are a = 1.96 and e — 0.091 GeV. 



As already discussed, the introduction of the g^-dependent cut in the spectral representation does not change the 
differential distributions at small but it strongly affects the region of large q^. In particular, the cut provides 
the vanishing of dT/dq^ at the physical threshold q^ = {Mg — Af/j)^ = 11.6 GeV^. Let us point out again that this 
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cutting procedure does not affect the integrated rate at the leading and subleading l/mg orders. Note also that the 
differential distributions dV/dq^ given by the dashed and solid lines in Fig ^(a), are equal to each other at = 
and match the OPE result for dT/dq'^{q^ = 0). 

As we shall see later, the improvements on the (/^-differential distributions by approximate account of higher order 
graphs affect only moderately the A/j^-distribution in i? Xcii^e (as well as the photon lineshape in the rare B —>■ Xs"f 
decay). The latter are thus determined by the i?-meson wave function. 

This property allows us to obtain a realistic energy distribution and other observables through the soft wave function 
of the heavy meson. Thus we do not need to introduce any unknown 'smearing function' describing the motion of the 
b quark inside the B meson, but rather directly calculate the effects of the b quark motion with the soft meson wave 
function. 



D. dV/dMx 

We now proceed to the calculation of the Mjf -distribution. 

Recall that in the free-quark decay, which is the leading-order process within the OPE framework, one finds the 
Mx distribution in the form 5{Mx — rric) (neglecting the radiative corrections). Inclusion of the l/mg corrections 
yields a singular series containing derivatives of the S function. For the interpretation of these results one needs an 
introduction of a smearing function. In the quark model the Mx-spectrum is obtained already smeared because of 
the Fermi-motion of the b quark in the B meson. Therefore the shape of the Mjf-spectrum is calculable through the 
B meson wave function determined in the Chapter II. 

Let us calculate dT/dMx fro m our spectral representation for the decay rate. To this end it is necessary to change 



the order of integration in Eq. ( 7.43 ) as follows (M^ = S2 in our notations) 



(mi+ma)^ s^{s2,q^) 

This representation immediately gives d^T /dq^dM^ ■ Integrating the latter over q^ one obtains dT/dM^. The 
calculated distribution is reported in Fig.^ Our result should be compared with the leading order OPE result 
6{mx — rric)- One can see that already the box diagram of the quark model provides a smooth and reasonable 
distribution (beyond the resonance region), which is only moderately affected by a proper account of the subleading 
l/ruQ effects. Most important is that at large Mx the calculated distribution does not require any additional smearing 
and can be directly compared with the experimental results. 
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E. The electron energy spectrum 



In distinction to the Lorentz-invariant distributions dV/dq^ and dV/dMx considered above, the lepton energy 
spectrum should be referred to a specific reference frame. 

The starting point of the electron energy spectrum calculation is the double spectral representation (7.28) for the 
inclusive semileptonic decay rate. We define the variable E to be the energy of the lepton li in the reference frame 
specified by the condition 

i5i = (Vir,0). (7.67) 
This is the rest frame of the (5293 quark pair with the invariant mass 



As the first step, let us consider the hadronic tensor w^^, given by Eq. (7.31). Using the expression (7.34) for the 
trace and performing the integration over dkidk2 leads to the following result 



ne{. . .) (si - (ma - ms)^) 



V*''/\ J- \ ' ' '' Z ''''oil 
W^.{puq) = ^2 „2N ^MPl^^) (7-68) 



where is represented in terms of the 'dispersion momenta' pi and q = pi — P2 (Pi = si, P2 = S2, ~ Q^) as 
follows 

,~ , Plf^Pl,^ , . ~ , Plt^q^ + Pll^Qt, , . ~ 

Wf_i„{pi,q) = -.g^iyiui H W2 + ie^,yai3—=qi3W3 H — W4 + qf_iqi,W5. (7.69) 

-^1 V'^l V -^1 



Lorentz structures containing the momentum k are omitted in this expression. 
The functions wi-a necessary for the energy distribution have the form 

wi = 4(toi +ml~ q^) - 16/3, 
q wi + —W2 = -CsL{m^,m2,q ), 



3 3 

«^3=8Vir(l-ai-a2), (7.70) 
with ai, Q!2, and (3 given by Eqs. (3.52- ^754 ). We have introduced 

k1 = ^^^^1^, 9" = ^4^, (7.71) 

the components of the momentum q in the reference frame (7.67) 

It might be interesting to note that at q^ < 0, and using the reference frame q^ — 0, pi± = (q^ = — f?^) one 
obtains 

P=-{kl- , l-ai-a2 = l-x,, (7.72) 

where 2:3 and k± arc the (+) and (_L) components of the spectator quark momentum, respectively. In the heavy 
meson one finds 

(3~A^, X3~A/mQ. (7.73) 

Altho ugh at > the interpretation of (3 and ai + a2 in terms of k± and 0:3 is not straightforward, the estimates 
( 7.73 ) remain valid also at q^ > 0. This means that in B ^ Xcii^i decay our Wi's differ only slightly from the 
corresponding free-quark expressions. 

Now everything is prepared for the calculation of the lepton energy spectrum. The following steps are standard: 
We must take the convolution of the trace over the lepton loop with the hadronic tensor w^i,, and integrate the result 
over the part of the lepton pair phase space leaving out the integration over dE. 

As a result, we obtain the following double differential distribution 

00 

Gl\V21? [ , 2t , gi-(m2-m3)\ i/2, 2 2^^(Sl) 

1 

J ^qo > E + {S2 < sr%q')) 

-1 

X {2q^wiisi,S2,q^) + [4£;(g° - E) - q^]w2isi,S2, q^) + 2q^{2E - g°)u;3(si, S2, ?')} • (7.74) 



109 



I n thi s formula S2 is connected with 77 through Eq. (7.39), and qq is given in terms of the integration variables by Eq. 



Eq. (7.74) includes modifications which tune the size of the l/mg corrections as explained above. 
Let us check that the integration of the double differential distribution ( [7.74 ) over E brings us to the expression 
(7.64) for dr/dq^ obtained earlier. 
By virtue of the relations 



dEe{q° > E + q^/4E) = \q\, 
EdE9{q° > E + q^/AE) = q^\q\/2, 



E^dE6{q'^ >E + q^/AE) = -\m<ll + f)- 



one finds 



dE{2q^wi + [AE{q° - E) - q^]w2 + 2q^{2E - q^)w:i}e{q^ > E + q^/4E) 



2\q\iq'wi + -(fw2). 



(7.75) 
(7.76) 
(7.77) 

(7.78) 



Thus, integrating the double differential distribution (7.74) over E and using ( 7.7^ ) gives dT/dq^ in th e for m ( [7.64 ). 
Clearly, for the calculation of dT/dq^ we do not need to know all Wi, but only their linear combination (7/7S). On the 
other hand, for calculating the electron energy spectrum d T/dE all functions wi_3 are needed. 

The electron spectrum dT/dE is obtained by integrating ( 7.74 ) over q^. Fig. ^ plots the results of the calculations. 
Here the effects of the subleading 1 /mg orders are more pronounced but nevertheless lead only to a moderate change of 
the quark model box-diagram result. We also compare the dispersion approach prediction with the electron spectrum 
in the free-quark decay process. At low values of the electron energy the spectrum in free quark decay and in the 
inclusive B decay are practically indistinguishable. As desired the spectrum is nonzero at large E above the free-quark 
threshold. At the intermediate E the spectrum in B decay is substantially than the spectrum in the free-quark decay, 
providing the known suppression of the integrated T{B Xclv) compared to ro(6 — *■ civ). 




> 
o 




Fig. 39. Distribution dT/dE in B ^ Xcive decays vs the lepton energy E: (a) ratio of the bound-to- free quark decay 
R{E) = {dT / dE) / [dTo/ dE), lines same as in Fig. (b) Our quark-model result (solid) vs free-quark result (dot-dashed). 
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F. Discussion 



In this Chapter we discussed the dispersion approach to the description of quark-binding effects in the inclusive 

decays of heavy mesons. Our approach allows to express kincmatical distributions in inclusive decays of heavy mesons 
in terms of the heavy meson soft wave function. This soft wave function describes long-distance effects both in 
exclusive and inclusive processes. 

Summarising the main results of this Chapter: 

1. Wo analysed the hadronic tensor in the quark model. We have shown that the diagrammatic representation of 
the hadronic tensor in the quark model yields an expansion in inverse powers of the heavy quark mass as well 
as the standard OPE. However in distinction to the standard OPE which is based on the free quark decay as 
the LO process, the LO process in the quark model is described by the box diagram with a free Qq in the final 
state. This yields some specific features of the hadronic tensor calculated within the quark model, both negative 
and positive. 

On the one hand, a consideration based on the box diagram alone reproduces the correct LO result but contains 
also l/mq corrections. These l/mq terms are cancelled against contributions of other graphs thus leading to 
the agreement with the standard OPE result. Hence the effects of the subleading diagrams should be taken into 
account for obtaining a consistent approach. 

On the other hand, the hadronic tensor calculated from the box diagram is a regular function in the kinematically 
allowed region of all variables, as well as all the subleading terms of order l/mg. This feature makes the quark 
model calculation of the hadronic tensor very suitable for describing the differential distributions; 

2. We have constructed the double spectral representation of the hadronic tensor for the B Xc£Pe decay in 
terms of the B meson soft wave function and the double spectral density of the box diagram, and analysed its 
l/mg expansion in the case of the heavy-to-heavy inclusive transition. The spectral representation is further 
modified in order to take into account essential effects of the other diagrams contributing in subleading orders. 
Namely, the subtraction term in this dispersion representation is determined such that the l/mg correction 
in the integrated semileptonic rate is absent in agreement with OPE. Furthermore, a phenomenological cut 
function is introduced into the spectral representation to bring the size of the l/jTig terms in the differential 
g^-distribution at = and in the integrated rate into full agreement with OPE. Thus our representation of 
the hadronic tensor obeys the OPE predictions in the regions where the latter are expected to be valid; 

3. We have obtained numerical results on differential distributions in inclusive B XciPe decays using the B- 
meson wave function and other quark model parameters previously determined from the description of exclusive 
meson transition form factors within the dispersion approach. So, basically our predictions are parameter-free. 
We notice that modifications of the spectral representations which take into account the subleading 1 /mg effects 
within the box-diagram representation, introduce some uncertainties in our results. However, they do not affect 
our predictions strongly, and the main features of the inclusive distributions are determined by the soft meson 
wave function. Moreover, the size of the subleading corrections is in perfect agreement with the OPE result 
for the integrated rate, and we expect to describe also these subleading effects in differential distributions in a 
proper quantitative way. 

The proposed approach can be applied to the inclusive B X^^s transitions. In particular it is especially suitable for 
the description of the photon line shape in B — > Xgj decays. However, certain subtleties in heavy-to-light transitions 
compared with the heavy-to-heavy decays emerge. They are mostly connected with the fact that in heavy-to-light 
transition the kinematically allowed g^-interval of the hadron semilaptonic decay is larger than that of the quark 
decay, while in case of the heavy-to-heavy transitions the situation is just opposite and the g^-region of the quark 
decay is larger. This feature requires a detailed analysis of the g^-region near zero recoil in heavy-to-light inclusive 
decays. 

It is also worth noting that our approach takes into account only non-perturbative effects in inclusive decays of 
heavy mesons. Perturbative corrections have been ignored. So, for comparing our results with the experimental 
differential distributions perturbative corrections should be also included into consideration. 
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VIII. CONCLUSION 



We presented a review of the relativistic dispersion approach to meson transitions induced by the weak currents. 

This approach is based on the constituent quark picture and takes into account the leading two-particle singularities 
of the Feynman diagrams making use of spectral representations over the mass variables. 
We reported the following main results: 

1. Relativistic meson Qq wave function has been defined which describes in a universal way strong-QCD effects 
related to the meson structure in various meson interactions. Due to confinement, the wave function in terms of 
the relative quark momentum is localized in the region of the order of the confinement scale. The normalization 
condition for the wave function corresponds to the electric charge conservation and is thus related to the meson elastic 
electromagnetic form factor at zero momentum transfer. 

2. We considered the elastic electromagnetic form factor of a pseudoscalar meson and form factors describing weak 
transition of a pseudoscalar meson to pseudoscalar and vector mesons at spacelike momentum transfers. Relativistic 
double spectral representations for these form factors have been constructed in terms of the double spectral densities 
of the corresponding triangle Feynman graphs and the wave functions of the participating mesons. 

The subtraction prescription has been fixed by requiring the dispersion approach form factors to respect all prop- 
erties of the transition form factors known from long-distance QCD in the limit of heavy quark mass. 

3. The dispersion approach form factors were studied in the two limiting cases of the heavy-to- heavy and heavy- 
to-light weak transitions. 

• Heavy-to-heavy transition 

In this case the masses of the initial quark mq and the final quark tuqi participating in the weak transition 
Q ^ Q' are taken to be much larger than the confinement scale A 

rriQ ~ mgi » A. 

The form factors from the dispersion approach have been shown to fully respect the structure of the l/mg 
expansion known from QCD in the leading and subleading l/mq orders. The Isgur-Wise function and the 
subleading-order universal form factors have been calculated in terms of the wave function of a meson containing 
an infinitely heavy quark. 

• Heavy-to-light transition 

In this case the quark masses satisfy the relation 

mg » ruQ' ~ A. 

The dispersion approach form factors have been shown to obey the QCD-based relations between the form 
factors of the vector, axial vector and tensor currents valid within the l/mg accuracy in the region of large g^. 

4. We performed a detailed comparison of the dispersion approach with the light-front quark model for the descrip- 
tion of weak transitions. Both approaches describe the dominance of the qq components of the meson wave function 
in meson interactions and therefore should lead to the same results. 

Indeed, we have shown that for the transition between pseudoscalar mesons the results from both approaches in 
the region of spacelike momentum transfers q"^ < are identical. On the other hand, for some of the form factors 
describing the pseudoscalar to vector meson transition the results are difi^erent. This is related to a different treatment 
of a vector meson, a composite system with spin-1, within both approaches. In the language of spectral representations, 
the difference between the two approaches can be traced back to different subtraction prescriptions. We have shown 
that the subtraction prescription implemented in the light-cone quark model does not fully respect the next-to-leading 
order constraints from QCD for the heavy-to-heavy transition. 

Another advantage of the dispersion approach is the possibility for a consistent description of the region of timelike 
momentum transfers > relevant for the decay kinematics. In the dispersion approach this is achieved by 
performing the analytical continuation in the variable q^ . The light-front quark model in general cannot consistently 
treat the form factor at q^ > 0: it describes only a reference-frame dependent part of the form factor, the partonic 
contribution, leaving out the nonpartonic contribution. 
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5. Wc obtained the form factors in the region > by the analytical continuation. We observed the appearance 
of the anomalous cut in the complex s-plane, where s is the invariant mass of the qq pair. Accordingly, the spectral 
representation for the transition form factor acquires the anomalous contribution in the decay region. The anomalous 
contribution is small for small but becomes increasingly important in the region of large q"^ near zero recoil. 

6. A model for pseudoscalar mesons was considered, which allowed us to study the transition regime to the heavy 
quark limit. Meson decay constants and transition form factors as functions of the quark masses have been analysed. 
It was found that the corrections to the leading order l/niQ relations can be as big as 10-20% for the quark masses 
in the region of the realistic c and b quarks. 

7. We presented the results for the form factors for weak decays of and mesons to light pseudoscalar and 
vector mesons. 

• The effective quark masses and meson wave functions were determined by fitting the quark model parameters 
to the available lattice QCD results for the B ^ p transition form factors at large momentum transfers and to 
the measured D — > {K, K*)lv decay rates. The knowledge of the quark model parameters allowed us to predict 
numerous form factors for many decay channels and for all kinematically accessible q^ values. 

In spite of the rather different masses and properties of mesons involved in weak transitions, all existing data 
on the form factors can be understood in the constituent quark picture, i.e. all form factors can be described 
by the few degrees of freedom of constituent quarks. Details of the soft wave functions are not crucial; only the 
spatial extention of these wave functions of order of the confinement scale is important. In other words, only 
the meson radii are essential. 

• The calculated transition form factors are in good agreement with the results from lattice QCD and from sum 
rules in their regions of validity. 

• Our predictions agree well with all available experimental data. An important example is the ratio of the 
semileptonic branching fractions of the B meson to light tt and p mesons, B{B pli/)/B{B ttIv). A good 
agreement of the predicted ratio with the results of the experiment leads to close values of the matrix element 
Vub extrtacted from the B ^ p and B ^ n channels. 

• We estimated the products of the meson weak and strong coupling constants by extrapolating the form factors to 
the meson pole. The value of each coupling constant can be obtained independently from the residues of several 
form factors. In all cases the values extracted from the different form factors agree with each other within the 
5-10% accuracy. They also agree with the results of the independent direct calculations of these decay constants 
within the dispersion approach. First direct measurement of the coupling constant go'D-w which has become 
available recently is in perfect agreement with our estimate. This gives additional argument in support of the 
reliability of our results for the form factors. 

8. Weak annihilation in a rare radiative decay was analysed. The relevant B jlv form factors were calculated in 
terms of the B meson wave function. Parameter-free estimates were obtained and found to agree well with estimated 
from pcrturbatibc QCD and QCD sum rules. The situation with contact terms was clarified. 

A new contribution to the weak annihilation amplitude was reported, related to the photon emission from the 
light-quark loop. This contribution was believed to vanish in the limit of zero quark masses, and therefore neglected 
in previous analyses. A detailed study shows that this contribution stays finite and should be taken into account. 

9. Non-factorizable effects in the B° — B^ mixing amplitude due to the soft gluon exchanges were studied assuming 
the local gluon condensate dominance. The corresponding correction to the factorizable amplitude was demonstrated 
to be negative. This correction contains specific B meson transition form factors which were calculated within the 
dispersion approach. Numerical estimate for the nonfactorizable effect was obtained. 

10. The dispersion approach was applied to inclusive semileptonic decays. Spectral representations for the differen- 
tial and integrated semileptonic decay rate T{B XJi/) were obtained. The decay rates calculated within the disper- 
sion approach satisfy an important property known from the operator product expansion in QCD in the heavy quark 
limit: namely, the order l/mg correction in the ratio of the bound to free quark decay rate r{B XJi')/To{b dv) 
is absent. 

Different kinematical distributions such as the dilepton q^-spectrum dT/dq^, dT/dMx and the electron energy 
spectrum dT/dEi were calculated in terms of the B meson wave function. This wave function determines the B 
meson form factors and has been tested in exclusive B decays. 
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